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PHOTOELECTRIC MEASUREMENTS OF THE A4200A CN BAND 
AND THE G BAND IN G8Ks SPECTRA 


R. F. Griffin and R. O. Redman 
(Received 1959 October 2) 


Summary 


A description is given of a rapid method of using a spectrometer to give 
reproducible photoelectric measures of integrated intensities over significant 
features of stellar spectra, without the necessity for a sky of photometric 
quality. The accuracy achieved is 1 per cent. The A4z00A CN band 
has been measured in 712 late-type stars and the G band in 212 of these. 

Absolute magnitudes derived from CN intensities have mean errors of 
at least 1-5. There is only a weak correlation of weak CN with high space 
velocity. ‘‘ Weak-line’’ stars have no anomaly in their CN absorption; 
‘* strong-line ’’ stars have slightly enhanced CN. 

For G to K1 stars the spectral type deduced from G band intensities 
has a mean deviation of 0-15 class from the MK type. There is a small 
absolute magnitude effect, the G band strengthening with intrinsic luminosity. 
High-velocity stars have completely normal G band intensities. ‘‘ Weak- 
line ’’ stars show the G band very slightly stronger, “‘ strong-line’’ stars show 
it slightly weaker than the mean, but the individual scatter is fairly large 
and the difference between the two groups is scarcely significant statistically. 





1. Introduction.—Ordinary stellar photometry follows an empirical procedure 
in which the light from a star is measured successively in three or more broad 
bands of wave-length, of the order of 1000A wide. These bands are isolated 
by the combined transmission characteristics of interstellar space, the Earth’s 
atmosphere, the telescope, and a colour filter, and the photometric response 
characteristic of a particular light receiver. Some of the properties of the Earth’s 
atmosphere, and usually of the apparatus also, are generally unknown and 
irreproducible, so that each observer measures in his own ‘‘ system”? (and indeed 
in a succession of different ‘‘systems’’ as he changes his apparatus). For the 
same reasons, the stars are measured relative to one another and not to a terrestrial 
source, and the results are not readily linked to the ordinary physical units. 
A measure of the star’s apparent brightness is obtained, but the procedure is not 
such as to provide the maximum amount of other information about the stars; at 
best one can sometimes get a rough idea of absolute magnitude, and sometimes 
an estimate of interstellar absorption. Fundamentally,the method isamodernized 
version of one devised by the ancient Greeks and it has a number of disadvantages ; 
an important one for us is that the accuracy of modern photometers cannot be 
exploited without a very good sky, which in Cambridge is experienced scarcely 
ten times per year. 

At the other extreme, one may make a detailed spectrophotometry, usually 
photographically. In its most elaborate form this requires a large reflector, a 
good spectrograph, photographic exposures sometimes running into hours, the 
production of microdensitometer records, and a reduction time which can run 
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into months. In a less ambitious form, the telescope can be smaller, the spectro- 
graph more modest and examination of the spectrogram confined to a quick 
visual inspection for estimates of spectral type and luminosity class. The more 
sophisticated procedures of this kind are expensive in telescope time and man- 
hours; and although the resulting information can be very detailed, a good deal 
of care is required to keep errors down to a few per cent. However, the method 
has the outstanding advantage that the sky quality need not be particularly 
good. 

As a kind of middle course between these two extremes, Strémgren and his 
pupils have measured stars photoelectrically through interference filters trans- 
mitting fairly narrow bands centred near selected features of the spectrum. ‘The 
results of measurements through two or more different filters can give estimates 
of spectral type and intrinsic luminosity with a high degree of internal consistency. 

2. Method.—We seek a quick way of obtaining quantitative information 
about stellar spectra which can be related, as an ultimate goal, to such parameters 
as absolute luminosities, chemical abundances, and stellar population character- 
istics, and hence to stellar evolution. At Cambridge during the past two years 
a method has been tried which combines some of the virtues of the other methods 
outlined above. ‘The principle is to isolate some significant feature of the spectrum 
by means of a spectrometer, using a wave-length band usually 30-50A wide, 
and to measure the total light within this band photoelectrically. Measurements 
are also made simultaneously upon two comparison bands, approximately sym- 
metrically placed on either side of the principal band. If the light intensities 
in the three bands are represented by a, 6, and c, the ratio (a+c)/b is a measure 
of the intensity of the spectral feature situated within the middle band, and varies 
from 2 upwards as the intensity of the (absorption) feature varies from zero 
upwards. In principle, a considerable number of features could be measured 
at one time in the same way, as indeed is done in laboratory spectra for rapid 
analysis of industrial products, where up to about thirty channels have been used. 

The transmissionof a spectrometer, and (usually) that of the Earth’s atmosphere, 
does not vary greatly over a 50 A bandwidth, nor does the response of a photo- 
multiplier cathode. Using a spectrometer, therefore, one measures over a fairly 
flat-topped band with a fairly sharp cut-off at either end. The gradient at the 
cut-off depends upon the spectrometer slit width and the seeing, or if the latter 
is very good upon the telescope focusing and guiding as well, but in every case 
it covers only a small fraction of the whole band. The fraction of light which is 
transmitted by the spectrometer may be less than that given by a very good inter- 
ference filter at its peak, but the spectrometer has many advantages; amongst 
these we may note that the transmission is practically constant over the whole 
width of the chosen band, the cut-off at each end can be positioned at will, and 
the wave-lengths accurately specified so that the measurements can, if desired, 
be repeated by others. Most interference filters have a transmission curve whose 
peak varies appreciably in wave-length from point to point (and their properties 
sometimes vary with the passage of time). It is difficult to make an interference 
filter to a tight specification, so that exact repetition of filter measurements by 
others is impracticable. A further important advantage of the spectrometer is 
that the comparison bands can be measured simultaneously with the central 
band, so that changing sky transparency has no more effect than it does in photo- 
graphic spectrophotometry. The chief disadvantage in comparison with the 
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photographic method is that the measured feature is to some extent diluted with 
other unwanted features. 

3. Apparatus and procedure.—The spectrometer used in the work described 
here is vodification of one designed by R. V. Willstrop for three-colour 
photomc: , in bands a few hundred angstréms wide; however, it has served 
for the experimental work on relatively narrow bands. It uses the first order 
of a Wood grating in an Ebert mounting of focal length 36 inches and gives a 
dispersion of 18-4A/mm. ‘The spectrometer is used at the f/18 coudé focus 
of the 36-inch reflector; the image scale is 12”-5/mm, and the image of the slit 
projected into the spectrum is therefore 1-5 A wide per second of arc in the 
star field. ‘The usual width of the slit has been about 2”; in average seeing the 
diameter of the star image is of the same order, and a large fraction of the light 
in it is passed into the spectrometer. This fraction is sufficiently nearly the 
same for each of the three adjacent bands measured simultaneously. The slit 
length is 12 seconds of arc. 

The spectrometer has no thermostat, and there are appreciable wave-length 
shifts with changing dome temperature. It is sufficient to correct these at 
intervals of 20 minutes upwards, according to external conditions. The 
wave-length adjustment is monitored by bringing a line in another part of the 
spectrum into coincidence with a crosswire in the focal plane of an eyepiece, 
using a gas discharge tube as source. The tolerated error is about 1 A. 

The light from each of the three bands of the spectrum transmitted by the 
diaphragm in the focal plane of the spectrometer is fed to an unrefrigerated 
EMI photomultiplier, the output being recorded by pulse-counting (1). After 
amplification, the output pulses from the photomultiplier pass to a discriminator 
which rejects those of less than a given amplitude and permits the acceptably 
large pulses to be counted, the result being displayed on a panel of dekatrons, 
Potentially, pulse-counting gives the maximum useful response from any 
photomultiplier, but in fact (since programmes have been restricted to stars 
brighter than 7™-5) there has usually been more than enough light, and no 
attempt has been made to get maximum electronic efficiency; indeed, neutral 
filters have been used to reduce the light from the brighter stars in order to 
avoid the effects of possible non-linearity of the count at high counting rates. 

The light response of a photomultiplier varies, sometimes grossly, with 
position on the photocathode. We have used Fabry lenses imaging the 
spectrometer collimator aperture, and therefore also the telescope aperture, 
as 4 mm diameter circular patches on the cathodes. Each patch naturally has 
a central shadow due to the Cassegrain secondary mirror of the telescope. 
Moreover, scintillation causes a variation in the distribution of light intensity 
over the illuminated area of each cathode, as part of the light in the star image 
formed by the telescope spills onto the jaws of the spectrometer entrance slit; 
and in addition, poor focusing of the star image upon the slit may cause further 
irregularity of illumination of the patch. 

Despite the fact that the instantaneous pattern of illumination is practically 
identical upon the several cathodes, it seems possible that in unfavourable 
circumstances—in particular, if there are extreme variations in sensitivity within 
the illuminated area of one or more cathodes—errors of a few per cent might 
arise from these causes. Difficulties recently encountered after further 
modifications to the spectrometer, involving the use of larger Fabry images 
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upon the cathodes, have been tentatively ascribed to such origins. (This work 
is still in progress and is not reported here.) Until these modifications were made, 
however, no evidence was found, in spite of various tests, that the effects described 
were of any significance. Although for various reasons the photomultipliers 
have been moved and readjusted quite a number of times, and the telescope 
mirrors (irregularity in whose reflectivities would imply corresponding irregu- 
larities in illumination of the Fabry images) have been twice realuminized, no 
suspicion of systematic change has been found in measurements made over 
periods of up to one year. 

In early trials the three wavebands were measured with one photomultiplier ; 
a chopping device gave a 4}-second ‘‘exposure’’ on each in turn. Thirty 
exposures (ten on each channel) constituted one observation, the total count 
being distributed appropriately among the three outputs. This arrangement 
had the advantage that any change in the photomultiplier sensitivity affected 
all three channels equally. Nevertheless it was inefficient because, quite apart 
from observations requiring three times as long as if the counts ran simultaneously, 
the frequent switching occasioned some additional noise and errors were 
introduced by changing sky transparency and, worse, by movements of the star 
image on the slit. By altering the position and shape of the image, seeing causes 
the fraction of star light which passes the slit to vary, and in the circumstances 
described this variation cannot be relied on to average out sufficiently accurately 
between the channels. 

When a change was made to three separate photomultipliers, working 
simultaneously, both speed and accuracy improved; but the relative sensitivities 
of the three channels may now alter with, say, temperature or mains voltage. 
These changes might be minimized by using exactly similar photomultipliers 
and electronic circuits in each channel; but owing to the manner in which the 
work developed this was not done. In any event, the relative sensitivities of the 
channels must be checked sufficiently often; in our case checks are made at 
intervals of 10-15 minutes. 

For this purpose a tungsten-filament lamp is used; it is adjusted to have 
about the same colour temperature as the stars under observation. Our practice 
has been to reduce observations of stars using lamp counts interpolated to the 
time at which any given star observation was made. If the star counts on the 
three channels (after subtraction of dark counts) are A, B, and C, and the 
interpolated lamp counts A,, B,, C,, then the intensity ratio we require is 


A C B 
(4. +7) /s,- 


In effect one is forming the ratio (a+c)/b, referred to previously, using not 
the absolute light intensity of the star in each waveband but its intensity relative 
to the lamp. As well as taking care of drifts in the relative sensitivities of the 
three channels, this procedure has the advantage that the two comparison bands 
are equally weighted in the ratio. Suppose, for instance, that the chosen bands 
are in a region of steep spectrophotometric gradient, or that the comparison 
bands are of greatly unequal width: then, although the amount of light in one 
comparison band will be much greater than in the other for the star, the same will 
be true for the lamp and the quotients A/A, and C/C, will be approximately 
equal. This condition is desirable since, when it is fulfilled, any change in the 
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colour of either lamp or star has only a second-order effect upon the intensity 
ratio; for if, owing to a colour change of the lamp or star, A/A, rises relatively 
to B/B,, then to a first approximation C/C, is reduced by a corresponding 
amount and the intensity ratio stands unchanged. (In the special case in which 
the amount of light in the comparison bands is approximately equal, a worthwhile 
simplification can be effected by receiving the light from both comparison bands 
upon a single photomultiplier.) Changes in the apparent colour of the star 
may be caused by differential absorption and also, in the case of our particular 
spectrometer, by differential refraction in the Earth’s atmosphere. The 
spectrometer slit is horizontal so that in observations near the meridian 
atmospheric dispersion lies across it, and movement of the star image on the 
slit is tantamount to altering somewhat the effective colour of the star. For 
the same reason, when observing bands in the violet it is appropriate to set the 
visual image of the star off to one side of the slit by a suitable amount; the effects 
of dispersion do not however become appreciable until zenith distances of 40° 
are reached. 

Changes in the general brightness of the star or lamp do not affect the intensity 
ratio. ‘The lamp therefore does not need to be controlled at all accurately and we 
have in fact used ordinary torch bulbs, behind two spaced diffusers of flashed 
opal, run from a 44-volt dry battery. The original battery is still in use. The 
spectrometer is diaphragmed so that the lamp uses as far as possible the same 
optical path as the star; the patch which the lamp illuminates upon the 
photomultiplier cathodes does not have a central shadow as does the star-lit 
patch, although even this could be duplicated by placing a central stop on the 
collimator mirror if it were deemed worth while. 

Experience has confirmed the most optimistic expectation in that gross changes 
in sky transparency do not affect the intensity ratio detectably. Satisfactory 
results are obtained through passing cloud and on extremely hazy nights with 
a zenith extinction as great as 2™-5. On such occasions there is sometimes 
pronounced reddening, the absorption in regions only 150A apart differing by 
(exceptionally) 20 per cent. Observations were once successfully made through 
a fog which reduced horizontal visibility at the telescope to less than thirty yards. 
Changes in the telescope transmission are equally without effect: realuminizing 
the mirrors resulted in a light gain of 1™-3 but there was no sign of a systematic 
change in the measured ratios. 

Owing to the small area of sky covered by the slit, and the fairly small 
bandwidths employed, the contribution of the sky background to the noise 
count is usually quite negligible; only in strong twilight or near a full moon has 
it been necessary to distinguish between ‘‘sky’’ and ‘‘ dark’’, and observations 
are practicable almost up to the limb of the Moon. This convenient state of 
affairs is chiefly a result of programmes so far having been confined to fairly 
bright stars. If the method is pursued to much fainter stars, as it certainly can 
be (by using selected photomultipliers, improving the counting circuits and 
possibly increasing the counting time) more attention to sky brightness will 
become necessary. 

4. Observations.—The routine finally adopted is as follows. Each count or 
‘‘exposure’’ lasts the same time, about one minute, and is terminated auto- 
matically. ‘Two such exposures are made successively on each star before the 
telescope is moved to the next. After every group of three stars a lamp count 
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is taken; sometimes this is repeated using a neutral filter in front of the lamp 
as a rough and ready check on the linearity of the photometric response. At 
longer intervals, as experience suggests, the spectrometer wave-length setting 
is checked and if necessary adjusted, and dark counts (and, in a bright sky, sky 
counts as well) are taken. Normally 12-16 stars are observed per hour; ona 
number of occasions over 100 stars have been measured in a night. 

Each star is usually observed on at least three different nights, i.e. at least 
six times. Where the results on a particular star have appeared less than usually 
consistent further observations on subsequent nights have been added, although 
from a statistical standpoint such a procedure is perhaps not to be recommended. 
Very few individual observations have been rejected; but whole nights have 
occasionally been rejected following unsatisfactory readjustments of the apparatus. 

The present paper deals with measurements of the A 4200 A CN band in over 
700 stars and of the G band in over 200 stars. Other work will be reported 
elsewhere: it includes measurements of the magnesium b lines in 539 of the 
same stars by ‘I’. J. Deeming, and short programmes on the sodium D lines by 
Deeming and on the A3800A CN band by G. A. H. Walker. 

5. The 44200 A CN band.—-The first programme carried out by the method 
outlined above was the measurement of the intensity of the 44200 A CN band 
in late-type stars, principally those of types G8-Ks5. All the stars lie between 
declinations o° and +73’, the upper limit being imposed by the presence of the 
spectrometer at the coudé focus to the north of the telescope. 

An attempt was made to include all stars, within these ranges of spectral type 
and declination, brighter than 5™-5; most of these have been listed by Roman 
(2). ‘To these were added dwarfs, subgiants, supergiants and high velocity stars 
collected from various sources, to a limiting magnitude of 7-5. Later a number 
of F stars, especially supergiants, were added to the programme after it had 
been found that high-luminosity stars show appreciable CN absorption as 
early as F5, although it has been stated that this absorption is undetectable 
in MK types earlier than Go (3). In the magnitude range 5™5~7™-5, 
preference was given to those stars for which MK spectral types were available. 

The purpose of the programme was to elucidate quantitatively the relationship 
between CN strength on the one hand and spectral type and luminosity on the 
other; and to investigate the divergences of individual stars from the mean for 
their type and luminosity, particularly with regard to the supposed weakening 
of CN in high-velocity stars. 

The wave-lengths of the three observed bands were 


A 4097-4149 A, 
B 4164-4214A, 
C 4230-4283 A. 


These bands were chosen to avoid the strongly luminosity-sensitive Sri lines 
at AA 4077, 4216A, the intense Cat line at 44227A and the G band at about 
AA 4285-4315 A. The A channel certainly includes a good deal of the weaker 
structure of the CN band, and also covers Hd, but this seemed unavoidable. 
It is inevitable, too, that owing to the complexity of a late type spectrum a large 
number of other more or less strong lines, including many of Fe1, are included 
in each of the three bands, and it should be made clear that the measurements 
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must to some extent be affected by this. Our observations simply define the 
degree in which the total absorption in the waveband AA 4164-4214A exceeds 
that in the specified adjacent regions; and we shall call the results CN intensities, 
because, although this term may be open to semantic objection, for most of the 
stars on the present programme the principal cause of the increased absorption 
around A 4200 A is nevertheless CN. 

The slit width used throughout the programme was 0-161 mm = 2:0 seconds 
of arc = 3:0 A width projected in the spectrum. 

A total of 5,032 observations was made on 67 nights in 1958 and one night 
in 1959, the number of stars observed being 712. ‘The mean CN intensity ratio 
derived for each of these stars is given in ‘Table I. 

The number of observations rejected was 65. ‘The deviations of the remainder 
from their respective means form an accurately Gaussian distribution with a 
root-mean-square value of 0°0425. Owing to the derivation of the means for 
all the stars from the same observations and the consequent reduction of the 
number of degrees of freedom by 712, the r.m.s. error of a single observation 
is slightly greater, 0-046. Since the average of the CN intensity ratios for all the 
stars measured is about 2:2, this error corresponds to 2-1 per cent; the theoretical 
value to be expected from the finite number of photoelectrons counted is about 
half this. The error of the mean CN ratio for one star, which in most cases 
depends on six observations, made in pairs so that only three are entirely indepen- 
dent, may be estimated at 1-o per cent. ‘The intrinsic scatter in CN strength among 
stars of the same spectral type and luminosity is found to be much greater than 
this, so that the errors of measurement can be ignored in what follows. 


Taste I 


Type CN G band Notes 
ratio ratio 


dG8 2°08 
G3 ll | 
Ko III ‘21 
dGg ‘In: 
F2 IIt 97: 
K4 Ill “28 
Fs Ib-II ‘98 
G8 Il ‘27 
dKs 13 
F3 II ‘95: 


- 
oc mu OMA WwW N = 
"NN N= BS NN N 


Gs Ill 2°10 
K2 Il 2°30 
K2 III 2°17 
Ko III 2°10 
K2 Ib-II 2°36 
Ks Ill 2°16 
K4 Ill 

K3 Ib-II 

G8 IIIp 2°08 
K3 ITI+ 
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TABLE I (continued) 


H.D. Type CN G band Notes 
ratio ratio 


Ko Vt 2°11 
Ko II-III 2°45 
K2V 2°09 
G8 III 2°19 
Gg III-IV 2°20 
Ks III 2°14 
Go Ibt 2°15 
K1 II 2°25 
sgKo 2'21 
K2V 2°01 


Ks III 2°16 
Go Vt 1°96: 
K3 Ia 

K2 Ill 

sgKr 2°26 
G8 III-IV 

G8 III-IV 

G8 III-IV 

Ko III 

Go Ia 


K2 Ill 
Gs VI 
Kr Il 
Ko III 
Ko III-IV 
Ko III 
Fo Iat 
Ko III-IV 
G8 IV 
Ko III 


K2 pIl 
Fs III 
F3 Ib 
Ko III 
K4 III 
G8 IIIt+ 
K3 III 
Go Ib 
cK1 
G8 III 


K1 Ill 
G8 II-III 
Fs Ib 

K3 IIIt 
G8 III 
G2 Vt 
K3 III 
Ki V 
sgK2 

Fs Iat 
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TABLE I (continued) 


No. Name H.D. Type CN G band Notes 
or H.R. ratio 


o Psc 10761 G8 Ill d 
511 10780 Ko V 
521 10975 Ko III a 
11092 Ks Iab-Ib 
a Tri 11443 F6 IV 
é Psc 11559 Ko III 
11544 Gz Ib 
56 And 11749 Ko III 
t Ari 11909 Kip 
12399 Gs la 


y' And 12533 K3 Il 
a Ari 12929 K2 IIIt 
60 And 13520 K4 Ill 
&' Cet 13611 G8 Il 
645 13530 Ko III 
13686 K3 Ib 
13725 K4 Il 
6 Tri 13974 GoV 
690 14662 F7 Ibt 
64 And 14770 G8 III 


65 And 14872 K4 Ill 

27 Ari 15596 Gs III-IV 
737 15694 K3 III 

14 Tri 15656 Ks III 
v Cet 16161 G8 III 
753 16160 K3 V 
743 15920 G8 Iil 
747 16024 cKs 

14 Per 16901 Go Ib 

39 Ari 17361 K1 III 


825 17378 As lat 
7 Per 17506 K3 Ibt 
17 Per 17709 Ks III 
17971 Fs lat 
861 17958 K3 Ib 
24 Per 18449 K2 Ill 
18391 Go lat 
885 18474 Gp 
918 18970 Ko II-III 
t Per 19373 Go V 


« Per 19476 Ko III 
w Per 19656 K1 III 
8 Ari 19787 K2 III 
949 19735 Ks Ill 
969 20123 Gs Il 
978 20277 sgG8 
991 20468 K2 II 
20630 Gs Vt 
20618 sgG5 
20644 K4 Ill 
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TABLE I (continued) 


Name H.D. Type CN G Band Notes 
or H.R. ratio ratio 


63 Ari 20893 K3 III 2°26 
a Per 20902 Fs Ibt 2°04 
o Tau 21120 G8 IIIt 2°24 
o Per 21552 K3 III 2°19 
5 Tau 21754 Ko II-III "34 

36 Per 21770 F4 IIIt ‘97: 
1085 22072 dG7 ‘O07 

14 Tau 23183 Ko III 20 
v Per 23230 Fs IIt 03 

23841 K2 Ill “16 


25056 Go Ib 
1242 25291 Fo IIt 
37 Tau 25604 Ko III 
1255 25602 Ko III-IV 
25893 dK2 
49 Per 25975 K1 III 
1270 25877 G8 Il 
1286 26311 cKs 
pe Per 26630 Go Ibt 
1327 27022 Gs IIIt 


nO - WO 
Ntnenunw 


NR NNN NDNDNN KN 
) wW Ww 
wu Oo 


Oo 
Iw 


y Tau 27371 Ko IIIt 
dKs5 

54 Per 27348 G8 III 
¢ Tau 27382 K1 Ill 
5 Tau 27697 Ko IIIf 
1390 27971 K1 Ill 
7 Tau 28100 G8 III 

75 Tau 28292 K2 Ill 
e Tau 28305 Ko IIIt 
6 Tau 28307 Ko IIIt 


NNNNNN NNN WN 


a Tau 29139 Ks IIIt 
58 Per 29094 G8 II 
3 Cam 29317 Ko III 
1533 30504 K4 II 
2 Aur 30834 K3 III 
o* Ori 31421 K2 Ill 
« Aur 31398 K3 Il 
m® Ori 31767 K2 II 
31782 Ko IV 
e Aur 31964 Fo lap 


N 
oo 


Ne NNNN NNN WN 


B Cam 31910 Go Ib 
104 Tau 32923 G4V 
33299 K1 Ib 
1684 33554 Ks Ill 
p Ori 33856 K3 III 
16 Aur 34334 K3 III 
A Aur 34411 G2 IV-V 
1720 34255 cK4 
109 Tau 34559 G8 III 
21 Ori 34658 Fs Il 
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TABLE I (continued) 


No. Name H.D. Type CN G band Notes 
or H.R. ratio ratio 
dG6 2°04 
o Aur 35186 K4 Ill 2°28 
¢ Aur K3 p 2°41 
1884 G3 Ib 2°20: 
¢? Ori G8 IIIp 2°03 
1908 Ks III 2°03 
1925 dK1 2°00 
24 Cam sgG9 2°12 
51 Ori Ko III 2°15 
sgK1 2°18 


G8 Iab 2°44: 
dK1 2°06 
G8 III 

G8 III 

Ko III 

K2 IIt 

Go V 

Gog II 

Ko III+ 

K1 Ill 


Ko III 
Ko III 
G8 III 
K2 II-III 
G8 III 
K2 Ill 
Mo Iabt 
Gs Ip 

Cc 


K4 Ill 


dK3 
K1 IIt 
sgG8 
Ko Iab 
Cs 
sgG7 
K3 III 
Gs Ibt 
Ks III 
G8 Ibt 


K1 III 

Fs III 

Ko III-IV 
Go V 

Kz III 
K4 III 
Ko III 

K3 III 
K2 III 
cK4 
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Name H.D. Type CN G band Notes 
or H.R. ratio ratio 


52071 K2 III 2°13 

w Gem 52497 Gs IIt 2°22 
2649 52960 K3 Ill 2°29 

54371 dG6 2°Or 

2692 54563 dG7 2°02 

7 Gem 54719 K2 III 2°38 

63 Aur 54716 K4 II-III 2°19 
18 Lyn 55280 K2 III 2°20 
56224 Kr Ill 2°30 

65 Aur 57264 G8 III 2°16 


57 Gem $7727 G8 III 2°12 

66 Aur 57669 Ko III 2°41 
« Gem 58207 Ko III 2°22 
e CMi 58367 G8 III 2°27 
y CMi 58972 K3 Ill 2°16 

65 Gem 59148 K2 III 2°35 
6 CMi 59294 K2 Ill 2°35 
2896 60318 Ko III 2°27 
v Gem 60522 Mo IIIt 
o Gem 62044 Kr Ill 2°23 


76 Gem 62285 Ks III 
« Gem 62345 G8 IIIt 
B Gem 62509 Ko IIIt 

61994 dGs 

81 Gem 62721 Ks III 

63410 G8 III 
14 CMi 65345 Ko III 

65583 G8 V 
3145 66141 K2 Ill 
x Gem 66216 K2 III 


14 
‘21 
25 
95 
‘08 
10 
‘21 
“95 
18 
24 


NNs NN DH NN N 


pe Cne 67228 G2IV 
55 Cam 67447 G8 II 
BCne 69267 K4 IIIt 
31 Lyn 70272 Ks III 
3306 FIIISs G8 Il 
71597 K2 III 
o UMa 71369 Gs Il 
71952 sgKo 
3360 72184 K2 III 
7 Cnc 72292 K3 III 


‘Ol: 
23 
18 
08 
20 
19 
13: 
II 
31 


NNNNNNNN NN 


v? Cne 72324 Gg III 
o Hya 73471 K2 Ill 
nm UMa 73108 K2 III 
34 Lyn 73593 G8 IV 
5 Cne 74442 Ko III 
«Cne 74739 G8 II 
e Hya 74874 Go III 
35 Lyn 75506 Ko III 
55 Cnc 75732 dKo 
p® Cnc 76219 G8 II-III 
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TABLE I (continued) 


No. Name H.D. Type CN G band Notes 
or H.R. ratio ratio 


¢ Hya 76294 Ko II-III 2°28 b 
3545 76291 Kr IV 2°19 
w Hya 77996 K2 II-III 2°29 
3612 77912 G8 Ib-II 2°27 
o' UMa 77800 Ks III 2°06 
7 Cnc 78235 G8 III 2°13 
78249 Ki IV 
78479 K3 II 
€ Cnc 78515 Ko III 
81 Cnc 79096 dG7 


~~ NN N 
Seven 
COO mH 


17 UMa 79354 Ks III 

3664 79452 G6 Ill 
79969 dK4 

«x Leo 81146 K2 Ill 

81192 G8 III 

A Leo 82308 Ks III 

& Leo 82395 Ko III 

6 Leo 82381 K3 III 
82394 cG7 

8@UMa 82328 F6 IV 


~ NNNNNN | N WN 
SSne5eusee 
eonuwnornn nw 


- 


82443 dGg 

24 UMa 82210 Gs IV 
10 LMi 82635 G8 III 

3809 - $2741 Ko III 
11 LMi 82885 G8 IV-V 
10 Leo 83240 K1 III 

3834 83425 K3 III 
27 UMa 83506 Ko III 
43 Lyn 83805 G8 III 

e Leo 84441 Go IIt 


Saas 


NW 
e+ wr 


ai 
na 
a 
rae 
na 
2°33 
va 
nat 
a 
ra 


- 
~“ 


14 LMi 84453 Ko IV 
84406 sgKo 
3881 84737 G2V 
pe Leo 85503 K2 III 
31 Leo 87837 K4 III 
89269 dGs 
y' Leo 89484 Ko III 
29 LMi go250 Ki Ill 
90572 sgKo 
B LMi 90537 G8 III-IV 


NNNNF NNN NN 


31 Sex gIoll sgKo 

48 Leo 91612 G8 II-III 

37 LMi 92125 G2 Il 
4165 92095 K3 III 

38 UMa 92424 K2 Ill 
4181 92523 K3 Ill 
4242 94084 K2 III 
4243 94132 dGog 

46 LMi 94264 Ko III-IV 

44 UMa 94247 K3 III 


NNNWNNNN DN WN WN 
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TABLE I (continued) 


Name H.D. Type CN G band Notes 
or H.R. ratio ratio 


46 UMa 94600 Kr Ill 2°21 29S: 
94549 dG8 2°13 2°75 
4264 94669 K2 III 2°20 
47 UMa 95128 Go V 1°99: 
58 Leo 95345 Ki III 
a UMa 95689 Ko III+ 
65 Leo 96436 sgG7 
ys UMa 96833 K1 III 
97561 G7 IV 
73 Leo 97907 K3 III 


€ UMa 98230 Go V 
v UMa 98262 K3 III 
98824 Ki III-IV 
56 UMa 98839 G8 IIt 
4404 99196 K4 Ill 
83 Leo 99491 dKo 
t Leo 99648 G8 II-III 
4452 100470 Ko III 
2 Dra 100696 Ko III 
92 Leo 101484 Ko III 


61 UMa IOI5§01 G8 V 
3 Dra 101673 K3 Ill 
x UMa 102224 
4521 102328 
4550 103095 

104556 
o Vir 104979 
4610 105043 


105475 
105631 


105963 
106714 
106760 
107328 
107383 
107469 
107950 
108225 
108381 
108722 


109317 
109358 
110833 
110914 
111028 
111067 
111812 
112033 
112989 
113092 
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TABLE I (continued) 


Name H.D. Type CN G band 
or H.R. ratio ratio 


e Vir 113226 Gg II-III 2°26 2°68 
41 Com 113996 Ks III 2°18 2°83 
114960 Ks III 2°28 2°91 
4997 115004 Ko III 2°34 2°80 
115404 dK3 2°02 aoa: 
59 Vir 115383 Go V 1°98 : 
5013 115478 K3 III 2°86 
115539 G8 III-IV 2°76 
70 Vir 117176 Gs V 2°59 
5102 117876 Ko III 2°69 


118643 cK3 2°89 

84 Vir 119425 K2 Ill ° 2°87 
v Boo 120477 Ks III . 2°80 
6 Boo 120539 K4 III 2°82 
§227 121146 sgK2 - 235 
7 Boo 121370 Go IVt 

9 Boo 121710 K3 III ° 2°84: 
5273 122742 G8 V : 

5302 123977 Ko III 2°83 
15 Boo 124679 Ko III 2°72 


a Boo 124897 K2 IIIp : 2°85 
A Boo 125351 Kr III 2°81 
125560 K3 III 2°88 
p Boo 127665 K3 IIIt 2°90 
128165 K3 Vt : 2°75 
128750 sgK2 2°86 
128902 K4 III 2°75 
129312 G8 III 2°73 
129336 G8 III 2°74 
129580 dKo 2°82 


129972 Ko III 2°73 
IZ11II Ko III-iV 2°81 
131156 G8 Vt 2°59 
131507 K4 III . 2°87 
131511 dKi 2°69 
132132 sgK1 2°83 
133124 K4 III 2°83 
133208 G8 II-III 2°73 
133165 Ko III 2°79 
133582 K2 Ill 2°85 


134190 G8 III 
135722 G8 IIIt 
136726 K4 Ill 
136512 Ko III 
136514 K3 III 
137107 Go V 
137759 K2 IIIt 
137704 K4 Ill 
138481 Ks III 
139195 Ko p 





R. F. Griffin and R. O. Redman, Photoelectric measurements Vol. 120 


TABLE I (continued) 


Name H.D. Type CN G band Notes 
or H.R. ratio 


139341 dK4 2°80 
¢ Boo 139641 G8 IV 2°70 
wb Ser 140538 dGs 2°55: 
a Ser 140573 K2 Ill : 2°89 
A Ser 141004 Go V 
5 CrB 141714 Gs III-IV . 2°60: 
w Ser 141680 G8 III ° 2°81 
p Ser 141992 Ks III ° 2°80 
« CrB 142091 Ko III-IV ° 2°78 
5924 142574 K4 Ill "06 : 2°78 


¢ Ser 142980 Ki IV 
e CrB 143107 K3 II 
5 Her 143666 Ko III 
pCrB 143761 G2V 
v Her 144046 F2 II 
144287 G8 V 
144579 dG8 
«' Her 145001 G8 III 
«* Her 145000 
6014 145148 Ko IV 


2°86 
2°82 
2°80 
2°50 


2°78 
2°71 
2°62 
280 
2°80 


N NNR = me = NN N 


7 CrB 145328 Ko IIIt 
14 Her 145675 dKr 
49 SerA 145958 dKo 
49 Ser B dKo 
€ CrB 147677 Ko III 
v? CrB 147767 Ks III 
6126 148293 K2 III 
7 Dra 148387 G8 IIIt 
6136 148513 K4 IIIp 
148653 dK2 


2°82 


2°67 
2°63 


NNNNN NNN DN WN 


Bb Her 148856 G8 IIIt+t 
6152 148897 G8 p 

29 Her 149161 Ks III 
6199 150449 K1 III 

¢ Her 150680 Go IV 

18 Dra I51I01 Ki p 

n Her 150997 G8 III-IV 
43 Her 151217 Ks Ill 
48 Her 151937 Kr II-III 
51 Her 152326 K2 II-III 


NON 


NN wHNNNN NN 


152391 G6 V 
6286 152812 K2 III 
6287 152815 G8 Ill 
54 Her 152879 K4 III 
6301 153226 Ko IV 
«x Oph 153210 K2 IIIt+ 

154345 dKo 
6342 154278 K1 Ill 
6364 154733 K4 III 
6388 155410 K3 III 


NNNNNNKNK NN WN 
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TABLE I (continued) 


No. Name H.D. Type CN G band Notes 
or H.R. ratio ratio 


am Her 156283 K3 IIt 2°31 2°92 
6433 156681 K4 II-III 2°18 2°80 
72 Her 157214 G2V 2°03 ; 

o Oph 157999 K3 II 2°33 2°90 
6518 158633 dKi 2°00 

A Her 158899 K4 III 2°28 2°88 
6 Dra 159181 G2 Ib 2°16 ' 
27 Dra 159966 Ko III 2°22 

160346 dK3 2°00 2°74 

83 Her 161074 K4 III 2°14 2°84 


B Oph 161096 K2 IIIt 2°36 2-91 
161198 dG8 268 
161796 F3 Ib 2°08; 
p Her 161797 Gs IVt 2°65 
6638 162076 sgG5 2°76 
87 Her 162211 K2 Ill 2°85 
go Her 163217 K3 Ill , 2°78 
€ Dra 163588 K2 IIIt 2°83 
89 Her 163506 F2 Iat 
6 Her 163770 K1 IIt 2°92 


y Dra 164058 Ks IIIt 2°87 
€ Her 163993 Ko III 2°71 
93 Her 164349 Ko II-III 2°89 
164922 Ko V 2°61 
70 OphA 165341 Ko Vt 2°68 
70 Oph B K7V : 2°69 
71 Oph 165760 G8 III-IV 2°76 
6791 166208 Ko p 2°53 
6793 166229 K2 III 2°88 
6806 166620 K2V 2°75: 


6817 167042 Kr III 2°78: 
6820 167193 K4 III 2°81 
6853 168322 Ko III 2°83 
105 Her 168532 K4 II 2°87 
« Lyr 168775 K2 III 
74 Oph 168656 G8 III 
108 Her 168913 Fo Ib 
6885 169191 K3 Ill 2°83: 
109 Her 169414 K2 Ill 2°85: 
42 Dra 170693 K2 III 


6950 170829 G8 IV 

45 Dra 171635 F7 Ibt 
6983 171779 Ko III 
6995 171994 G8 IV 
7008 172365 F8 Ib-II 
7064 173780 K3 Ill 
o Dra 175306 Ko II-III 
7123 175225 dG8 
7137 175535 G8 III 
v Dra 176524 Ko III 
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TABLE I (continued) 


Name H.D. CN G band Notes 
or H.R. ratio ratio 


e Aql 176411 2°41 2°82: 
718% 176527 2°20 2°85: 
A Lyr 176670 2-37 2°92: 
7275 179094 2°18 aya: 
53 Dra 180006 3 2°33 
179784 2°40 
179785 2°34 
5 Dra 180711 ; 2°20 
180028 J 2°09 
54 Dra 180610 2°26 


7308 180583 2°10 
6 Lyr 180809 Ko IIt 2°38 
«x Cyg 181276 Ko IIIt 2°28 
23 Aql 180972 K2 II-III 2°34 
7345 181655 G8 V 2°02 
182293 K3 III a°at 

182296 G3 Ib 2°33 

7368 182488 Ko V 2°09 
31 Aql 182572 G8 IV 2°14 
4 Vul 182762 Ko III 2°20 


7421 184010 sgG8 

7428 184398 K2 II-III 
184467 K1 V 

p Aql 184406 K3 III 

o Dra 185144 Ko Vt 

e Sge 185194 G8 III 

7468 185351 Ko III 

7475 185622 cKs 
185662 dK4 

¢ Cyg 185734 G8 III-IV 


15 
15 
‘Ol 
29 
‘02 
33 
‘22 
26 
‘19 
24 


NNNN NNN NN WN 


a Sge 185758 Go IIt 
B Sge 185958 G8 IIt 
186120 sgKo 
10 Vul 186486 G8 III 
15 Cyg 186675 G8 III 
y Aql 186791 K3 IIt 
7542 187203 Go Ib 
187299 Gs Iab-Ib 
e Dra 188119 G8 III+ 
20 Cyg 188056 K3 Ill 


15 
35 
‘37 
‘22 
25 
-30 
17 
"39 
7: 


“41 


NNNNNNNN DNDN 


7573 187982 A2 la ‘95: 
& Aql 188310 Ko III ‘23 
188507 K4 II-III 18 
B Aql 188512 G8 IVt 05 
188753 dKo “s2: 
7 Cyg 188947 Ko III 29 
7633 189276 Ks II-III 2°17 
y Sge 189319 Ks Ill 2°15 
189671 G8 II 
26 Cyg 190147 Kr II-III 2°38 
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TABLE I (continued) 


No. Name H.D. Type CN G band Notes 
or H.R. ratio 


16 Vul 190004 Fs II 
331777 F8 Ia 
7670 190360 G6 IV 
190403 Gs Ib-II 
190404 Ki V 
p Dra 190940 K3 Ill 
7 Sge 190608 K2 III 
27 Cyg 191026 Ko IV 
66 Dra 191277 K3 II 
191046 Ko III 


NNNNNH ND D 
~NeNWOW & ww 
ous n Oe OI 


— 
+ 


191785 Ki V 
192004 K3 II-III 
192577 K2 II 
192713 G2 Ibt 
192806 K3 Ill 
192909 K3 Ib-II 
192944 G8 III 
193092 K4 Il 
193217 K4 Il 
193370 Fs Ibt 


by NNNNNNDNDHDN 
~NNNH HK NERO 
CONS HVMROGSH 


193469 Ks Ib 
194069 G2 Il 
194093 F8 Ibt 
194013 G8 III-IV 
194317 K3 Ill 
194526 Ks III 
195100 - Gs Ill 
195338 G7 Il 
195295 Fs IIt 
195506 K2 III 


2°23 
2°29 
2°11 
2°27 
2°25 
2°16 
218 
2°30 
2°01 
2°17 


195593 Fs Iabt 
195987 Go V 
196093 K2 Ib 
196725 K3 Ib 
196755 Gs IV 
196758 K1 III 
' 197345 Az lat 
197752 K2 III 
197912 Ko III 
197913 dGog 


N wre oo 
Ssz oroun 


NNN NNNDNH DN 
w 
w 


ON 
nau 


198149 Ko IVt 
197989 Ko III 
197964 sgK1 
198134 K3 III 
198387 dKo 
198809 G8 III 
199191 G8 III 
199169 K4 III 
199253 Ko III 
199612 G8 II-III 
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TABLE I (continued) 


Name H.D. Type CN G band Notes 
or H.R. ratio ratio 


199580 K1 IV 2°15 

200102 G1 Ib 2°22 

€ Cyg 200905 Ks Ibt 2°19 
8082 201051 sgK1 2°21 
61 CygA 201091 Ks Vt 1°89 
61 Cyg B 201092 K7 Vt 1°81 
63 Cyg 201251 K4 Il 2°28 
C Cyg 202109 G8 II 2°39 
203358 G8 IV-V 2°08 

8165 203344 Ko ITI-IV 2°28 


1 Peg 203504 Kr III 2°30 
71 Cyg 204771 Ko III 2°28 
8248 205349 K1 Ib 2°44 
p Cyg 205435 G8 III 2°27 
72 Cyg 205512 K1 Ill 2°35 
206078 G8 III 2°06 

205941 dG8 2°08 

206121 Gs II 2°10 

25 Aqr 206067 Ko III 2°23 
206312 K1 Il 2°43 


206374 G8 V 2°03 
206731 G8 II 2°30 
206952 Ko III 2°41 
206778 K2 Ibt 2°42 
206936 M2 Iat 2°11: 
206859 Gs Ibt 2°40 
207130 Kr III 2°36 
207119 Ks Ib re 
207089 Ko Ib 2°36 
207260 A2 lat 2°01: 


207134 K3 III 2°26 
207991 Ks Ib 

208606 G8 Ibt 

209747 K4 Ill 

209960 K4 III 

209945 Ks III 

210220 G6 Ill 

210459 Fs II-III 

210807 G8 III 

210745 Kr Ibt 


210889 K2 Ill 
211073 K3 III 
211076 K4 Ill 
211388 K3 II-III 
212466 ? la 
212496 Go III 
212943 Ko III-IV 
213310 Mo Iab 
213893 Ks III 
214868 K3 III 


NNNNNNNN DN 
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TABLE I (continued) 


No. Name H.D. Type CN G band Notes 
or H.R. ratio ratio 


n Peg 215182 G2 II-III ae: 
13 Lac 215373 Ko III 2°30 
215549 K1 III-IV 2°10 

A Peg 215665 G8 II-III 2°39 
p Peg 216131 Ko III 2°21 
« Cep 216228 K1 III 2°25 
8692 216206 G4 Ibt 2°28 
8726 216946 Ks Ibt 2°15 
8752 217476 Go lat 2°17 
8761 217673 K2 Il 2°36 


8779 218029 K3 III 2°40 
3 And 218031 Ko III 24 
8784 218101 G8 IV 14 
56 Peg 218356 Ko IIp es | 
4 And 218452 Ks III "aa: 
60 Peg 218395 Ko IV 14 
8832 219134 K3 Vt 08 
y Psc 219615 G7 Ill "15 
8857 219668 sgKo "32 
o Cep 219916 Ko III 16 


NNNNNNN NWN 


‘22 
: 
34 
"a5 
34 
‘97: 
33 
-_ 
*20 
"10 


11 And 219945 Ko III 
219978 Ks Ib 
8875 219962 K2 III 
7 Psc 220009 K2 Ill 
66 Peg 220363 K3 III 
v Peg 220657 F8 IV 
@ Psc 220954 Kr Ill 
70 Peg 221115 G8 III 
14 And 221345 Ko III 
221354 Ko V 


N NNN FN NNN DN 


221585 G8 IV 
221639 sgGo 
221673 K4 Ill 
221861 Ko Ib 
222107 G8 ITI-IV 
222842 Ko III 
223047 Gs Ib 
223075 C6 

223173 K3 II 
223165 Ki Ill 


05 
17 
24 
42 


NNN N 


25 
+31 
“58: 
°34 
28 


NNN NHN DN WN 


224014 ?la 26: 
224617 F4 III 96: 


Notes to Table I 


(1) The order is that of R.A. for epoch 1960. 

(2) No observations of the following northern Gs5-Ks stars brighter than 5™-5 were 
secured: HR 1205, 47 Tau, 4 UMa, 6 Dra, 24 ComA, e€' Boo, 3 Ser, 95 Her B or 
f} Cyg; or upon the following, which are in declination > +73°: HR 285, 49 Cas, 
HR 1523, HR 2527, HR3751, HR 4126, 4 UMi, 5 UMi, 6 UMi, @ UMi, e UMi, 
HR 7117, 7 Dra, HR 8702, HR 8748, 7 Cep, y Cep. 
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(3) Meaning of symbols: 


Mean of observations made on fewer than three nights. 
MK standard spectrum. 

Dynamical space velocity greater than 80 km/sec. 

“* areo™. 

‘“* Weak-CN ”’. 

** Strong-line ’’. 

“** Weak-line ’’. 

B.D. + 42° 939. 

Compcsite spectrum, early-type companion. 

Double, companion star only partly excluded by slit. 


Treen Op +. 


In the cases of all other double stars, where the angular separation is less than 2” the 
two stars have been measured together, and where the separation is more than 2” 
the brighter component alone has been observed. 

(4) The CN anomaly of a class III star may be obtained by subtracting the following 
quantity from its CN ratio: G8 2-19, Ko 2:23, Ki 2-25, K2 2-28, K3 2-27, K4 2-21, 
Ks 2:13. The G band anomaly may likewise be obtained by subtraction of: G8 2-71, 

Ko 2°76, K1 2-80, K2 2°84, K3 2°85, K4 2°83, K5 2-80. 


The CN ratios given in Table I are compared in Fig. 1 with results published 
for some of the same stars by Lindblad (4), Ohman (5), Iwanowska (6), Keenan 
(7), Plassard (8), and Gyldenkerne (9). 

Although Lindblad used the method of exposure-ratios on objective-prism 
spectra, modestly remarking (4), ‘‘It is, of course, freely admitted that the 
method of exposure-ratios is a fairly crude spectrophotometric method in the 
case of bright stars’’, his observations show the best correlation with ours. 

The earlier measurements just cited and used in constructing Fig. 1 are not 
exactly comparable with each other or with ours. Most of them are on a 
logarithmic or magnitude scale and most are based on slitless or objective prism 
spectra, although Gyldenkerne used interference filters. There are differences 
in the wave-length bands used and for instrumental reasons there is much 
variation in the sharpness with which these bands are defined. The rather large 
scatter seen in Fig. 1 is therefore not unexpected. The tail of the distribution 
in the cases of Plassard and Gyldenkerne is due to stars of spectral types K4 
and Ks. 

5.1. Variation of CN intensity with spectral type and luminosity.—For most 
of the programme stars an MK spectral classification is available. An attempt 
was made to assemble more accurate absolute magnitude data, from the various 
published catalogues of spectroscopic parallaxes or luminosities, but this effort 
only served to convince us that the most reliable and homogeneous collection 
of these data lay in the MK classifications themselves, so that these have been 
used wherever they exist. 

The CN intensity ratios were first plotted against spectral type at each whole 
luminosity class. The mean curves are shown in Fig. 2. It is important to 
notice that these curves simply summarize the average results for all stars; 
individual stars show quite large deviations from the appropriate curve, the 
r.m.s. deviation being about 0:07, i.e. over 3 per cent. The deviations for 
luminosity class V, where there is very little CN absorption at any spectral type, 
are rather less; and rather remarkably those of the supergiants of luminosity 
class Ib are also smaller than the average. Since a scatter of 0°07 is over three 
times the error of measurement, it must be due either to an intrinsic scatter in 
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Fic. 1.—CN ratios (abscissae) compared with intensities derived by previous authors, 
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the CN intensities or to errors in the MK types and luminosities. However, 
our measurements of the G band (see Section 6) suggest that the MK spectral 
classification errors are not normally greater than about one-tenth of a spectral 
class; they would need to be three times as great to explain the CN scatter. 
Alternatively, mean errors of a whole luminosity class would be needed. We 
think that such an explanation is untenable and that the major part of the CN 
deviations represents a real scatter. 


26 














+9 © 


Fic. 2.—Mean curves of CN ratio as a function of spectral type for each luminosity class. 


Owing to this large scatter it is not possible accurately to deduce the absolute 
magnitude of a given star from its CN ratio, even if the spectral type is accurately 
known. If one attempts to derive what we might call ‘‘ CN absolute magnitudes ’’, 
the r.m.s. difference of these from the MK luminosity is of the order of 1™-5 
even in the range of spectral type G8-K1 where the sensitivity of the average 
CN ratio to luminosity is particularly favourable. The CN luminosity is not 
rarely greatly at variance with the value reliably determined by other methods, 
or in direct conflict with the trigonometrical parallax. In quoting examples we 
do not need to fall back upon well-known peculiar stars such as ¢? Orionis. The 
CN absolute magnitude for « Dra (G8III), for instance, is 4-4, which would 
require a parallax of o”-120 compared with the value of 0”-001 + 0-006 observed 
trigonometrically. Again, the consensus of opinion of several authors prevents 
us from accepting the absolute magnitudes of about — 2 predicted for y, 5, « and 
ma Tau. The high accuracies (o™-5 or even o™-3) claimed by previous authors 
(4, 8) for CN absolute magnitudes appear to be fictitious and due in part to 
observational selection of a rather homogeneous group of stars—which we have 
intentionally avoided as far as possible—and in part to the lack of alternative 
luminosity data for comparison purposes. For instance, comparison of Lindblad’s 
(4) derived magnitudes with the MK data reveals a standard error not far short 
of two magnitudes. 

5.2. Correlation of CN intensity with space velocity.—In order to correlate 
CN intensity with space velocity it is desirable first to remove the grosser effects 
of spectral type and luminosity. This is readily performed by subtracting from 
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the measured CN ratio for each star the mean value for all observed stars of the 
same spectral type and luminosity: the difference we shall call the CN anomaly. 
Space velocity data are available (2, 10, 11) for almost all of the 342 luminosity 
class III stars on the programme but are lacking for many of the others, so that 
the present discussion will be restricted to normal giants. It is difficult to say 
just how accurate these data are, but in so far as the velocities depend upon 
proper motions and upon absolute magnitudes estimated spectroscopically, mean 
errors of perhaps 25 per cent or more can be expected. 
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Fic. 3.—Frequency distribution of CN anomalies for luminosity class IIT stars. 
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Fic. 4.—CN anomalies as a function of space velocity for class IIT stars. Crosses distinguish 
““ weak-CN ”’ stars, circles “ 4150"” stars. 


The frequency distribution of CN anomalies for all observed giant stars is 
shown in Fig. 3. High-velocity stars are plotted below the horizontal axis, 
low-velocity stars above; the dividing velocity is taken as 80 km/sec. In Fig. 4, 
CN anomalies are plotted against space velocities directly. 

It is quite clear that there is nothing like a one-to-one correspondence between 
high velocity and weak CN absorption. Some high-velocity stars do have weak 
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CN, but in many others the CN ratio is abnormally high. Many low-velocity 
stars, on the other hand, have weak CN. There is indeed a small statistical 
difference in the CN ratio between the two groups of stars shown in Fig. 3, viz. 
0042 + 0010 (s.e.), which is little more than half the random scatter within either 
group. 

These conclusions are surprisingly at variance with what may have become 
established as a supposed fact. It was early found (12) that certain objects which 
displayed spectral peculiarities, sometimes including a weakening of the CN 
absorption, had high velocities, and it perhaps came to be supposed that such 
peculiarities were. a quite general characteristic of high-velocity stars (11, 13); 
and objects in whose spectra these characteristics were detected have sometimes 
been labelled ‘‘high-velocity’’ irrespective of their space motions (14). To 
apply this term to stars merely on the basis of spectroscopic peculiarities and not 
on a measured velocity appears to us most undesirable. Especially is this the 
case when the peculiarity concerned is so weakly correlated with velocity as 
is CN. We emphasize that not all weak-CN stars have high velocities and not 
all high-velocity stars have weak CN. 

It is worth pointing out that if for any reason the intrinsic luminosity of a star 
is overestimated or its parallax underestimated, its computed space velocity will 
be too high and also the CN intensity will be weaker than expected. One might 
in this way obtain a spurious correlation between high velocities and weak CN. 
It is not, of course, suggested that this is the whole story, since many stars have 
radial velocities which are themselves sufficient to warrant inclusion of these 
bodies in the high-velocity group; but it is an effect that will tend systematically 
towards a weak correlation of the observed sign. 

5.3. Weak- and strong-line, weak- and strong-CN stars.—Miss Roman (2) 
and others (11) have classified stellar spectra on the MK system and have noted 
in particular cases an abnormal weakness or strength of the CN band. Miss Roman 
has in addition detected spectroscopic differences leading to her well-known 
classes of ‘‘ weak-line’’ and ‘‘strong-line’’ stars. ‘These classes have aroused 
interest because they have been associated with the velocity characteristics of 
their members and hence with their membership of Populations I or II, or at 
least of distinguishable population groups. ‘The spectroscopic characteristics 
are supposed to reflect variations of chemical composition and these may be 
related to stellar age or the stage of evolution reached. 

Owing to a certain amount of confusion we must quote briefly what Miss Roman 
(2) means when she classifies G8-Kz1 stars into ‘‘ weak-line’’ and ‘‘ strong-line’”’ 
groups. She says that 85 per cent of all G8 to K1 stars can be assigned to one 
of two groups typified respectively by HR 4126 and 2 Draconis. The HR 4126 
stars have A4227A Cal and the G band both weaker than average, but their 
velocity distribution is similar to that of the ‘‘strong-line’’ stars of late F and 
early G types. She describes the HR 4126 group as ‘‘strong-line’’ by this 
analogy, ‘‘ although the name is no longer descriptive’’. Conversely, the 2 Dra 
stars have strong A44227A Cal and G band but are associated in velocity 
characteristics with the weak-line stars of earlier type and are described as 
‘“‘weak-line’’ stars. The remaining 15 per cent of the G8-Kz1 stars share the 
spectroscopic characteristics of the latter group but also have the A4200A CN 
band either unusually weak, being then called ‘‘ weak-CN ”’ stars, or unusually 
strong, in which case they are designated ‘‘4150’’ stars. This last name has 
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received objection as ‘rather ambiguous, for some objects designated as *‘ 4150 
stars ’’ in the literature (e.g. y Cep) actually do not have very strong CN bands’ 
(14). 

The nomenclature in this subject seems to us to be greatly in need of 
improvement. 

Our own measurements show that ‘‘4150’’ stars in fact have positive CN 
anomalies, and that all but two of the ‘‘ weak CN ”’ stars have negative anomalies, 
the exceptions being 16 Vir and 7 Cep. However, many stars with equally or 
even more anomalous intensities have passed without comment. 


J 
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(mean) (mean) 











Fic. 5. 


Of the giant stars observed by us, 265 are in Miss Roman’s catalogue (2). 
Nine are noted by her as ‘‘weak-CN”’ stars and their average anomaly is 
—o-051 +0018 (s.e.). Only three of the nine have anomalies more negative 
than this, but in our own measures as many as fifty stars altogether, 19 per cent 
of the total, have such anomalies, while in 14 of these the anomaly exceeds twice 
the quantity just cited. The mean anomaly of the eleven “‘4150’” stars is 
+0:087+0°014; five of these have greater anomalies, but so have nine other 
stars whose excessive CN intensity apparently has passed unnoticed. 

The ‘‘ weak-line’’ stars show no systematic anomaly, but the “‘strong-line ”’ 
stars have slightly enhanced CN absorption, their mean CN anomaly being 
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+0°039+0°007. The anomalies of the individual stars are shown in histogram 
form in Fig. 5. 

6. The G band.—During an intermission in the CN programme in 1958 
April-June some measurements were made on the G band, using the same method 


and apparatus and observing some of the same stars. The wave-length regions 
used were 


A 4230-4270A, 
B 4285-4315 A, 
C 4342-4380A. 


The slit width was 0-121 mm=1":5=2:2A. Since the consistency of the early 
measurements was rather lower than expected, all observations made in the 
period 1958 April 8-23 were rejected en bloc. The total number of later 
observations was 1259, made on 23 nights, on 212 stars. Two only of these 
observations were discordant and were rejected. The rest fit a Gaussian error 
curve very well; the r.m.s. width is 0-047, implying a root-mean-square error 
per observation of 0-052, or slightly better than 2-0 per cent. 
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Fic. 6.— Mean curve of G band ratio as a function of spectral type for luminosity class ITT. 


The mean G band ratio for each of the observed stars is shown in Table I. 
125 of these stars are normal giants; the variation of G band ratio with spectral 
type at this luminosity is shown in Fig. 6. Since the ratio reaches a maximum 
in the K types, it is not by itself a good criterion of type there. Upto K1, however, 
the spectral type predicted from the G band has a mean deviation of 0-15 class 
from the MK type; this is satisfactory but not unexpected, as the MK types 
rely heavily on the G band in the range G8-K1. Assuming the accuracy of the 
MK types, the r.m.s. deviation per star in G band ratio from the mean line drawn 
in Fig. 6 is 0-05. It is very unlikely that this scatter is due wholly to errors of 
measurement; nor can it be due to errors in MK type since it persists at the 
same value throughout the range of types studied, including the mid-K’s where 
the G band ratio reaches a rather flat maximum and type errors are unimportant. 
At least part of the scatter must therefore be regarded as intrinsic, as in the case 
of CN; but it is small enough to make the type corresponding to the G band 
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intensity a useful first approximation to the MK type for G stars and, it may well 
prove, in the F’s as well. 

The mean ‘‘G band anomaly’’ for fifteen high-velocity stars of luminosity 
class III is 0-000 + 0°013 (s.e.); seven of these stars have positive anomalies, 
seven negative, and one zero. The G band anomaly is plotted as a function of 
space velocity in Fig. 7. It seems clear that, contrary to what has been supposed 
(11), there is no general systematic dependence of the intensity of this band 
upon stellar velocity. The mean anomaly of the 23 ‘‘weak-line’’ stars is 
+0°'006+0-010, and that of the 22 ‘‘strong-line’’ stars is —0-023 +0°012. 
The difference between the two groups is therefore 0:029’+ 0-016, which is 
hardly of statistical significance. The scatter within the groups is much greater 
than the difference between them. 
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Fic. 7.—G band anomalies as a function of space velocity for class III stars. 


There are too few stars of luminosity classes other than III and whose G bands 
have been measured to permit any detailed conclusions to be drawn from them. 
So far as the data go, they show a weak dependence of G band ratio upon 
luminosity, the higher ratios corresponding to the higher luminosities; the 
dependence upon spectral type at each luminosity is indistinguishable in form 
from the curve shown in Fig. 6 for luminosity class III. The mean vertical 


Tas_e II 
Luminosity class Displacement No. of stars 
II +0°07 +0°015 7 
II-III +0°015 +0°02 9 
III-IV 0700 +0°02 11 
IV —0°035 +0°025 8 
Vv —0'07 +0015 12 
distance through which this curve should be displaced to represent the G band 
intensities at other luminosities is shown in Table II. At the season when the 
observations were made, practically no apparently bright stars of class Ib were 
in the easily accessible part of the sky. 
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PRELIMINARY LIST OF RADIAL VELOCITIES FOR 33 STARS 
NEARER THAN 20 PARSECS 


B. R. Leaton and B. E. J. Pagel 


(Communicated by the Astronomer Royal) 
(Received 1959 September 22) 


Summary 


‘To extend the available data on space motions of stars within 20 pe, 
140 spectra of standard velocity stars and stars from Gliese’s catalogue were 
taken with the 80 A/mm Xf spectrograph camera of the Mt Wilson 60-inch 
telescope by the Astronomer Royal during 1958 February-March, Using 
effective wave-lengths from 12 plates of § stars with quality “ a’’ velocities 
in R. E. Wilson's catalogue, radial velocities have been determined for 33 
stars with no previous determinations and for 17 stars with previous deter- 
minations of quality poorer than “a"’. The scatter among indfvidual plates 
is sometimes as large as 10 km/sec, but the mean from 3 or more plates 
should be accurate to +4 km/sec within the wave-length system employed. 
A further series of observations of this type is planned. 





1. Introduction.—A large amount of observational material exists for the 
small region of the galaxy in the immediate vicinity of the Sun, but there are a 


number of deficiencies. For instance, in the best available collection of data for 
nearby stars, Gliese’s Katalog der Sterne naher als 20 Parsek fiir 1950.0 (Gliese 
1957), of 1094 stars listed whose parallax is supposed to be 0”-o50 or greater, 
over 300 stars have, as yet, had no determination made of their radial velocity, 
even allowing for binary stars which may be presumed to have a common motion. 
Hence, although there are observations of the proper motions of these stars, it 
is impossible to deduce their space motions. It is obvious that any conclusions 
drawn from analyses of such motions (e.g. Woolley 1958) would be greatly 
strengthened if such deficiencies were made good. 

Most of the stars in Gliese’s catalogue for which there are, as yet, no measures 
of radial velocity are late-type dwarfs or sub-dwarfs, whose relative faintness 
has made such observations over-tedious or impracticable with even moderate- 
sized telescopes. Advantage has been taken of an opportunity to use the X-spectro- 
graph on the 60-inch telescope at Mt Wilson to make the data on nearby stars 
more complete*. ‘The comparatively large aperture of the 60-inch coupled with 
the speed of the four-inch X-spectrograph camera (f/1) renders possible the 
production of a usable spectrogram, dispersion 80A/mm, in a reasonably short 
observing time, e.g. a tenth-magnitude Mo in about one hour. 

This paper gives the results of measurements of spectrograms taken by 
Dr R. v. d. R. Woolley in 1958 February and March. It is intended to make a 
further series of observations of this type in the near future. 

*I am very much indebted to the Carnegie Institution of Washington for appointing me a 


Research Associate in 1958 and to Dr I, S. Bowen for making the 60-inch telescope available to 
me for securing these spectra. R. v. d. R. Woolley. 
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2. Summary of measurements.—The plates have been measured by Woolley 
(RW), Leaton (RL), Pagel (BP), and Alexander (JA) as shown. The whole 
series of measurements was analysed by Leaton, with some help from Pagel, 
using adopted wave-lengths (Table 1) derived from measurements of stars of 


‘ ” 


velocities having quality ‘‘a’’ in Wilson’s catalogue (Wilson 1953), which are 
given in Table II. Results for other stars with published radial velocities are 
given in Table III, and new velocities for 33 stars with no previous measurements 
available are shown in Table IV. 


TaBLe I 
Adopted Wave-lengths recommended Original 
wave-lengths, A. by PAM* Victoria systemt 
4455°03 4454°91 
4435 °28 35°23 
4415°26 15°15 
4404°75 04°75 
4383 °82 4383 °56 
4352°08 51°85 
51°74 
4325 °32 25°64 
4308 °05 07°93 
4290°00 
4271°69 4271°56 
4260°52 60°48 
4254°7! 54°35 
4°54 
4250°30 50°47 
4226°81 26°74 
4191°37 
4187°37 
4143 °60 
4101°66 (G—K3) 4101°74 
407183 4071°70 
‘52 (G—K3) 
4063 2, (K4—M) 
4045 °72 4045°76 
4030°67 4030°70 
4005 °35 4005 °30 


* Petrie, Andrews and McDonald (1958) + Harper (1938) 


Taking means from Table I, we find 


adopted minus Victoria = +0-014 A+0-02(p.e.) 
adopted minus PAM = +0001 A+o-014 (p.e.) 


Thus the present system agrees with that of Petrie, Andrews and McDonald 
within o-1 km/sec, which is very satisfactory. 

Wave-lengths derived by averaging all stars of ‘‘known’’ velocity (i.e. 
Wilson ‘‘a’’, “‘b”’ and ‘“‘c’’ and some measured by Dyer (1954)) gave a system 
of wave-lengths that are low on PAM by 


0°036 A +0014 (p.e.), 


i.e. these give velocities that are 2-5 km/sec greater. 
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Of the 140 plates, 
12 are of 5 stars with ‘‘a’’ velocities in Wilson ; 
34 are of 17 stars with some sort of measurements available (‘‘b’’, ‘‘c”’ 
and Dyer); 
73 are of 33 stars with no previous measurements available ; 
1g were not measured because of overexposure, underexposure, lack of 
focus, or double Fe lines; 
I appears to be the wrong star ; 
I remains unidentified. 


‘The measures for the first three categories appear in Tables II, III, and IV 
respectively. 


TABLE II 
Stars with “a” velocities in Wilson 
‘tiie Wilson Wilson Gliese Plate Exp. Measured Mean 
ae No. — velocity No. Xf time velocities * R.G.O 
7 CasA 3° 447 +940 34A 34370 r™ +11 (RL) 
(GoV) 3438 + — 7(RL) +2+3 (p.e.) 
3460 3™ + 3(RL) 


« Cet “{ +19°3a 137 3440 %r™ +15 (RL) +15 
(Gs5V) 


119 Tau 4: 22°5 3332 +21 (RL) +15 (RW) +18+1 (p.e.) 
(cM2) 33360 +21 (RL) +16 (RW) 


36 Com 5: 7 3342 2(RL) —1 (BP) 
(gMo) 3347 4(RL) +2(BP) +3+1 (p.e.) 
3357 4m 7(RL) +6 (BP) 


& Boo A 4 +3°9a 566A 3456 ym (RL) 
(G8V) 3457 408 3 (RL) +5+1 (p.e.) 
3491 ym (RL) 


Mean Wilson “ a’’—RL (km/sec): +0°5 +1°3 (p.e.) 


3. Quality of the results.—The best indication of the accuracy of the results 
is given by an inspection of ‘Table II. Different plates of one star occasionally 
give startlingly large discrepancies, perhaps owing to the increased effect of 
asymmetrical guiding on the wide slits employed in modern spectrographs. 
The effect is worst of all among the stars used as standards, since the exposures 
on these were very short; experience with other spectrographs suggests that 
5 minutes is the minimum exposure time for adequate guiding, while the brightest 
standards used here were exposed for 1 minute only. However, discrepancies of 
the order of 10 km/sec still exist between plates of one star that were exposed for 
longer than § minutes, e.g. for 7 Cas B (Table ITI). 

Thus the average of three plates should be correct with the probable error 
stated in Table IV within the present system of wave-lengths. Further observations 
with the Xf spectrograph are, however, needed to check the system itself. Values 
obtained from single plates may be out by as much as 10 or even 15 km/sec. 
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We are grateful to the Astronomer Royal for securing the observational 
material on which this paper is based and to the Director of the Mt Wilson and 
Palomar Observatories for making the 60-inch telescope available for this purpose. 
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H1, Her AND Her INTENSITIES IN PLANETARY NEBULAE 
M. ¥. Seaton 
(Received 1959 September 15) 


Summary 


Recombination spectra have been calculated for H1, Heit and the 
zp-nd series of Het. Expressions are given for the calculation of He*?/H* 
and He'/H* abundance ratios. It is shown that the ionized region of a 
nebula will contain a negligible amount of neutral helium if 


I (He 11, A 4686)/7 (He 1, A 5876) 
exceeds 1-2. 

For NGC 7662 and 2392 the measured relative intensities within each 
spectrum are in good agreement with theory but for NGC 7027 the agreement 
is poor for the weak lines and for the infra-red. ‘The suggestion that the 
intensity anomalies for NGC 7027 may be due to collisional processes 1s 
found to be difficult to accept. "The possibility is considered that there 
may be systematic errors in the intensity measurements. 

‘The He/H abundance ratios, by numbers of atoms, are found to be 0-20 
for NGC 7662, 0:16 for NGC 2392 and o-19 for NGC 7027. 


1. Introduction.—\t is usually assumed that the H1, Het and Hetr spectra 
in gaseous nebulae are excited by radiative recombination. The aim of the 
present paper is to check the theory by comparing observed and calculated 
relative intensities within each spectrum, and if these checks prove satisfactory, 
to obtain improved estimates of the He/H abundance ratios. We consider three 
planetary nebulae for which detailed spectrophotometric observations have 
been made. For NGC 7662 and 2392 the measurements by Minkowski and 
Aller (1) are found to be in good agreement with theory but for NGC 7027 the 
measurements by Aller, Bowen and Minkowski (2) are found to be in disagreement 
with theory for the weak lines and for the infra-red. We can find no satisfactory 
explanation in terms of physical processes in the nebula and therefore consider the 
possibility of systematic errors in the observations. 

2. Theory 

(a) Ht and Heii.—In earlier calculations (3) it was assumed that the al 
states were populated in proportion to the statistical weights (2/+ 1): 

Now Cs Dy. (1) 

Calculations by Burgess (4) and by Searle (5) show that (1) is not true but that 

calculations made assuming (1) nevertheless give quite good results for the 
observed intensities, 

I(n,n')= ¥ I(nl,n'l’). (2) 


u 
We make the usual distinction between Case A (optically thin in Lyman lines) 
and Case B (optically thick in Lyman lines). In Case B transitions to the ground 
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state are not taken into account, it being assumed that these are counterbalanced 
by absorptions from the ground state. It has been shown (6) that, in calculations 
assuming (1), the errors are greater for Case A than for Case B, the reason being 
that in reality transitions to the ground state take place only from the np states 
and that the importance of such transitions is exaggerated when (1) is assumed. 
For bright planetaries we may expect Case B to apply and calculations based 
on (1) to give good results for the H1 and Het! lines. In the present paper only 
Case B is considered. 

For the nn’ line the absolute quantum emission rate (quanta cm~* sec!) 
is denoted by 9,,,,.. For a recombination line of the ion X*+™ we have 


Yn,n(Xt™) = NN(X*™* Joey, p(X") (3) 


where N, is the electron density, N(X*+™*!) the density of recombining ions and 
x, (X*") is the effective recombination coefficient for the nn’ line. From formulae 
given by Seaton (6), 
: o ; 
Snow = 2% ZN *be7n ale a * re) 
? nn'(n* —n'*) 
where Z is the ion charge before recombination, 
2°77) 2a4ca,? 
QG= — a =5°:197 x 10 'cm*sec 


(« being the fine structure constant, ¢ the speed of light and a, the Bohr radius), 


- -_ 
A= ed = 157890 (7) (T in “K) (6) 


and x, =A/n*. ‘The factors b,(A)e™ are tabulated by Seaton (6) and the g,,_, 
by Baker and Menzel (3). 

Table I gives coefficients for recombination to the hydrogenic tons H® and 
He’; this includes ground-state recombination coefficients z,, coefficients 
z, defined by 


~R= > a, (7) 
n=2 


and effective recombination coefficients x, .(H"), %3¢,2,(H°) and a, ,(He*). 
Values of x, .(H®°) and «,, .,(H°) are obtained from the work of Burgess (4), 
allowing for /-degeneracy; all other coefficients are obtained from the work of 
Seaton (6, 7). 
The relative intensities in each spectrum are insensitive to the electron 
temperature. In subsequent discussions we consider only Case B and 
T,.=1'5 x 10° °K. 


‘The results of Burgess are used for the Ha/Hf ratio and the results of Seaton 
for all other H1 and Heli ratios. 

(6) He1.—For He® orbital degeneracy is removed and it becomes essential 
to take proper account of the individual m/ levels. Since (1) is not true calculations 
made assuming (1) will give incorrect intensities for the n/-+n'l' lines. Such 
calculations will also be over-sensitive to transitions to the ground state, and the 
2s metastable states. ‘This explains the poor agreement, noted by Aller and Menzel 
(8), between observed singlet to triplet intensity ratios and those calculated by 
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TABLE I 


Recombination coefficients in units cm® sec~ 
T,=1 x 108 
10'4a, (H”) 15°5 
10'*a,, (H°) 25°8 
10'*a4,9 (H°) 2°99 
10'4a54,2» (H") 5°27 2°38 


10'*a, (He*) 63°6 44'5 
10a, (Het) 154°0 gor! 
10'a,,, (He*) 20'8 11°8 


Goldberg (9) assuming (1) for m>2. Further evidence that Goldberg’s calcula- 
tions are inexact was obtained by Wurm (10), who noted that the deduced He/H 
chemical abundance ratios depended systematically on the He*/He*? ratios. 

Greatly improved He! calculations have been made by Mathis (11). These 
give singlet to triplet ratios in satisfactory agreement with observation and ratios 
He/H showing no systemmatic dependence on He*/He** but the calculations 
still do not make full allowance for all of the individual n/ levels. 

Most of the observed Het lines are of the type ms>2p and nd->-2p. We 
consider only the nd->2p transitions. ‘The d states, and all states with /> 2, are 
nearly hydrogenic. The rate of population of a level nd is determined almost 
entirely by capture on nd and by capture on states with /> 2 followed by cascade 
to nd. Allowing for statistical weight factors, a good approximation should be 


Xp ip, p(He®) = 4%, d, 2p(H") } (8) 


tn *D,2 (He?) wi 3%, 4,2)(H?). 


Some error will result from the fact that the ratios 
A, d, 2p! > Ana, n’p (9) 
~. 


will be different for He® and H® but this error is negligible for 3d->2p and is 
probably small for the higher states. We use values of ~,,4,,(H°) obtained 
from the calculations of Burgess (4); these will be accurate for small n but less 
accurate for large n. Finally it may be noted that the «,, »,(H°) are insensitive 
to transitions to the ground state (Cases A and B giving results which never 
differ by more than two per cent) and that we may expect ~,,:,, »p to be insensitive 
to transitions to or from 1 1S and 2'S and z,,,,, »p» to be insensitive to transitions 
to or from 2 3S. 

(c) He/H abundance ratios—For abundance determinations we consider 
the lines: HB, 42, 4861; Hel, 3 *7D—2 °P, A5876; Hert, 4->3, 14686. 

The relative quantum emission rates are 


(4686) _ J %2a(He*)N_V(He**) av 


44861) | a4, (H°)N, N(H*) dV 





q(5876) | a3 sp.» ap(He®)N,N(He*) dV 
q(4861) { 





a »(H°)N.N(H*) dV 
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where the integrals are over the nebular volume. ‘The ratios of recombination 
coefficients are insensitive to 7’. For constant temperatures and densities 
we obtain the abundance ratios in terms of intensity ratios: 
N(He**) <oasgeeet (12) 
N(He*) (4861) 
N(He*) = 1-00 {45876) (13) 
N(H*) 1(4861) . 
where the numerical factors are calculated for T,=1°5 x 10 K. ‘The chemical 
abundance ratio is given by 
N(He) _ \(He*?) + (He) (14) 
N(H) N(H*) 
if there is a negligible amount of He" in the ionized region. We now show that 
He® will have a negligible abundance if 
1(4686) eile (1s) 
1(5876) 
Consider the case of a pure He nebula, assume complete absorption of all 
stellar quanta capable of ionizing He® and let .4 (He®) be the number of stellar 
quanta capable of ionizing He® and .4 (He*) the number capable of ionizing 
He*. ‘The number of He* ionizations due to quanta produced in the nebula 
is equal to the number of recombinations to the He* ground state and the number 
of ionizations due to stellar quanta is therefore equal to the number of recombina- 
tions on excited states; hence 


WN (Het) = [ xy(He*)N.N(He**) dV (16) 





with x, defined by (7). In the same way one obtains 
MW (He®) = | 1,,(He®)N,N(He*) dV (17) 


noting that nearly all recombinations to He* produce one quantum capablet of 
ionizing He". ‘The integrals in (16) and (17) are proportional to observed quantum 
emission rates. Using the approximation x, (He°)=,(H°) we obtain* 
1(4686) _ inal (He ) (18) 
1(5876) A (He’*) 
and hence /(4686) /(5876)< 1-2. If the ratio is greater than 1-2 we may conclude 
that absorption by He® is incomplete and for typical dimensions and densities 
of planetary nebulae this implies a negligible amount of He’. 
For the pure He nebula a ratio exceeding 1-2 implies a limit on the amount 
of gas available to be ionized. For a real nebula, containing a mixture of H 
and He, this is no longer true but a ratio greater than 1-2 still implies a negligible 
amount of He® in the ionized region. Fig. 1 illustrates two possible situations. 
In Fig. 1(a) the He* region is limited by He® absorption and the intensity ratio 
must be less than 1:2 while in Fig. 1(6), which shows the position for higher 
star temperatures (13), the extent of the He* region is limited by H® absorption 
and the intensity ratio may exceed 1-2. 
If the ratio 1(4686)//(5876) is less than 1-2 the position may be that of Fig. 1(a) 
and in this case a reliable estimate of the He/H abundance ratio can be made 
only by examining the sizes of monochromatic images of the nebula in H1 and 


* A similar relation for /(4686)//(4861) has been discussed by Wurm and Singer (12.) 
+ But some of these quanta may ionize H”. 
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(a) Fic. 1. 


Het radiations. Mathis (11) tinds that a number of nebulae have similar sizes 
in Hf1 and Het. It would be interesting to know if this is the case for objects 
having very small values of the intensity ratio. 

3. Relative intensities for NGC 7662 and 2392.—We consider the intensities 
measured by Minkowski and Aller (1) and ignore all lines for which the intensities 
may be uncertain due to blending. Let /,(A) denote the observed intensities 
and /,(A) the intensities which would be observed if there were no interstellar 
reddening; then 


log I,(A) = log I,(A) + ef(A) (19) 


where f(A) is a known function (4) and c a constant which determines the amount 
of reddening. The function f(A) is chosen to be such that 10° is the factor by 
which the absolute Hf intensity must be multiplied in order to correct for 
extinction. Let [,(A) denote the calculated recombination intensities. Our 
procedure is to determine c on comparing J, and J, for the stronger Balmer lines 
and then to compare /, and /, for all of the H1, Het and Hert lines. Table II 
gives values of c and Table III gives values of J,, 7, and I, for NGC 7662 and 
2392. The values of J, and J, are made to agree for H1 A 4861, Het A5876 and 
He 11 A 4686; intensities for these lines are given in italic type. 


Tasce II 


Values of the reddening constant 

NGC c 

7662 0°46 

2392 0°32 

7027 1°22 

The agreement between theory and observation is seen to be good and it is 

doubtful if there are any differences greater than probable errors of measurement. 
To see if there is any suggestion of systematic errors, in Fig. 2 we plot log (/,/1,) 
against log J, and in Fig. 3 we plot log (/,/J,) against A, in both cases for NGC 7662. 
With the single exception of the Paschen line furthest in the red (A 8863) the 
values of log (/,//,) are never greater than 0-06 or less than — 0-18. Fig. 3 suggests 
that the values of log/, may be overestimated by about 0-15 for A< 4000 but 
in this region all lines considered have small value of J, and will be difficult to 
measure accurately. 
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Taste III 


Intensities for NGC 7662 and 2392 


NGC 7662 NGC 2392 
Series A z. I 
Hi 8863 . 
iu— 3 751 
S665 


H1 6563 y 380 294 257 

N—+2 4861 100 100 r00 
4340 441° 48°5 49°38 
4102 28 32°5 29°0 
3835 a g°0 
3798 5°6 69 «3 
3771 45 6-0 50 
3750 29 3°6 40 
3734 21 2°6 32 
3712 1°5 1°9 2‘1 
37°04 1'4 18 1°97 
36097 ' sid = 
3092 . . : —_ 
3687 : : : si _ 
3683 . : . oo 


es 


Naw UIE 
Neb DOs) Nv 


4471 ° . ° ° 2°9 
4026 q a ° 
3520 


6678 
4922 
4388 
4144 


4686 


5412 
4542 
4200 
3858 


tog I, 


Fic. 2.—Log (1,/1,) against log 1, for NGC 7662. 
(Hi: @, Heri: X,Heu: £.) 
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a = -_ : > 
4000 5000 6000 7000 


Fic. 3.—Log (1,/I,) against 4 for NGC 7662. 
(Hi: @, Her: x, Hen: @.) 

4. Relative intensities for NGC 7027.—We consider the intensities /, 
measured by Aller, Bowen and Minkowski (2) and ignore lines for which the 
intensities may be uncertain due to blending or are marked as uncertain for 
any other reason. ‘Table II gives the value of ¢ obtained from the stronger Balmer 
lines and ‘lable IV gives values of 7,, J, and J,. It is seen that there is fair agree- 
ment between /, and /, for the stronger lines but that the overall agreement 
between /, and /, is far from satisfactory. ‘The general trend of the discrepancies 
is such that the measured intensities are too strong for small 7, and too weak in 
the infra-red. 

We are faced with two possibilities: either (a) that the observed intensities 
are correct but that the recombination theory is insufficient or (6) that the dis- 
crepancies are due to errors in the observations. The first of these possibilities 
has been discussed, for the case of the H1 and Het lines, by Aller, Bowen and 
Minkowski (2) and by Aller and Minkowski (14). 

(a) Collisional effects.—If the recombination theory is suspect another method 
of determining reddening must be sought. Aller and Minkowski (14) considered 
the ratios of Paschen to Balmer intensities for lines having the same upper state, 
these ratios being calculated assuming the condition (1) to be satisfied. This 
method gives a much smaller reddening constant, c~ 0-6 compared with c= 1-22 
from the stronger Balmer lines assuming recombination theory. With c~o-6 
the corrected decrement for the stronger Balmer lines is much steeper than the 
recombination decrement. Aller and Minkowski suggested that this was due to 
collisional excitation from the ground state and showed that this requires N, ~ 4N, 
where N, is the density of H1s atoms. This very large concentration* of neutral 
hydrogen is difficult to reconcile with the high degree of ionization for other 
elements (14) and it is also difficult to see how the condition N,~ 3N, could be 
maintained; with N,~10'cm~%, linear dimensions of order 2 x 10!7 cm and an 
absorption cross section of order 10~' cm the optical depth would be of order 
10* for the Lyman continuum and ionizing radiation could not penetrate. One 
cannot overcome these difficulties by supposing (14) T, to be much greater 
than the accepted value (~ 1-5 x 10 K) since this would give a Balmer discon- 
tinuity much smaller than that observed. 

Whether or not one assumes some collisional excitation intensity, anomalies 
remain for the weak lines and for the continuum (1) It has been suggested 
that the anomalies for the weak lines may be due to collisional redistribution 
effects (1, 14). Approximate calculations by Aub and Seaton (to be published) 
suggest that, with V,~10'cm~* and 7,,~ 10' K, such effects will be important 


* Due to the use of Born collision cross sections the required value of N, will be underestimated. 





No. 4, 1960 Hi, Hei and Het intensities in planetary nebulae 333 


in H for n~ 100. It seems very improbable that they will be important for n ~ 15 
as required to explain the intensity anomalies. The idea of collisional redistri- 
bution may also be criticized on more empirical grounds. Such processes would 
tend to make the usual 5, factors (3) tend to unity, corresponding to conditions 
of thermodynamic equilibrium. With any reasonable 6, values for the lower 
levels the observed intensities for H1 and Heri give 6, which, for large n, are 
greater than unity and are increasing rapidly with n. 

(6) Possible errors in the observed intensities.—We assume that the recom- 
bination theory gives correct relative intensities and that the measured intensities 
are correct for the strong lines; we then have c=1-22. Assuming that there are 
calibration errors for the weak lines and in the infra-red we may deduce a correcting 
factor €=(J/,/J,). We consider the possibility of putting 


€=,(,)E9(d). | (20) 


In Fig. 4 we plot log (/,//,) against log J, for all lines of ‘Table IV having A < 8000 
(that is, we exclude the Paschen lines and He 11 A8237). For /,<5 we see that, 
as 1, decreases, (/,/1,) decreases systematically and shows an increasing amount of 
scatter. The fact that similar results are obtained from the various series of 
H1, Het and Het would be expected if calibration errors depend systematically 
on J, but would not be expected if the anomalies were due to the recombination 
theory being inapplicable. We take it that the dashed curve of Fig. 4 gives us 
the function log €,(/,). | Intensities corrected for systematic errors which depend 
only on J, will then be given by €,(/,)J.. 











1 
2 





tog |. 
Fic. 4.—Log (1,/1,) against log I, for NGC 7027 
(Hi: @, Her: xX, Hen: @,.) 


To examine the position in the infra-red, in Fig. 5 we plot log (/,/£,/,) against A. 
To reduce scatter we exclude the weakest lines (J,< 1-5). We include points 
obtained from continuum intensities. In the spectral region of Fig. 5 a good 
first approximation is to take the continuum intensities /,, per unit frequency, 
proportional to exp(—/v/kT.). We use the Mt Wilson continuum intensities 
of Minkowski and Aller (1). The intensities measured per unit wavelength 
are converted to intensities per unit frequency and corrected for reddening 
with c=1'22. For the continuum we take €,(/,)=1 and take (/,//,)=1 in the 
5000 A region. 
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TABLE IV 
Intensities for NGC 7027 


Series n 
Hi 10 
n —» 3 11 

12 


13 


3 
4 
5 
6 
8 
9 
10 
II 
12 
13 
15 


17 
18 
19 
20 
21 
22 
23 


—~— -— — - NN 


7 ¢ = 
ossd 


3 
4 
6 
rj 
8 
9 
° 


SN OM aw 
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TasB_e [V—(continued) 
Intensities for NGC 7027 


A n I, 
4686 4 45 


¥ 
} 


5412 
4542 
4200 
3924 
3858 
3834 
3814 
3782 
3769 
3749 
3740 
3733 
3715 


= NN 
Nu os 


NS 8 AO AW ~7 











Fic. 5.—Log (I,/£,1,) against A for NGC 7027 
(H1: @, Het: X, He: g, continuum: O). 


Fig. 5 suggests that, for A>6000A, all measured intensities are under- 
estimated by a factor which decreases with increasing A. We take it that the 
dashed curve of Fig. 5 gives us the function log &,(A)._ Aller and Minkowski (14) 
consider it possible that the Paschen intensities may have been underestimated 
due to the continuum having been drawn too high; we consider it more probable 
that all intensities have been underestimated in the infra-red. 

(c) Further discussion of the continuum intensities.—In Fig. 6 logarithms of 
continuum intensities, per unit frequency, are plotted against o=104/A. The 
full line curve is calculated (13) allowing for: H1 and Heil recombination and 
free-free continua using asymptotic expansions for the Gaunt factors; the Her 
recombination and free-free continuum, assumed hydrogenic; the H1 two- 
quantum continuum. The calculations are made for 7,=1-7x10'°K and 
N,=2:2 x 10'cm~’, which are the best estimates obtained (13) from forbidden 
line data and Hf brightness. The intensities measured by Minkowski and 
Aller (1) are corrected for reddening with ¢ = 1-22 and, in the infra-red, multiplied 


25 
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by the correction factor £,(A) of Fig. 5. The measurements are in arbitrary 
units and we therefore have one adjustable scale parameter when comparing 
theory and observation. Fig. 6 shows that this may be chosen so as to obtain 
very good agreement both for long wave-lengths and for the short wave side 
of the Balmer limit but that the measured intensities then show a depression 
to the long wave side of the Balmer limit. Referring to Table IV we see that 
it is in this region to the long wave side of the Balmer limit that there are the 





2.0 
o= 10/4 


Fic. 6.—Continuum intensities for NGC 7027. 
(Mount Wilson observations OQ, McDonald observations Xx). 
biggest discrepancies between J, and J, for the very weak lines, particularly those 
of Het. We suggest that, from about 4500A to the Balmer limit at 3647 A, 
the continuum may have been drawn too low and that in consequence the intensities 
of very weak lines may have been overestimated. 

(d) Conclusions on NGC 7027.—The intensity list of Aller, Bowen and 
Minkowski (2) for NGC 7027 contains many weak lines which have not been 
measured in other nebulae. NGC7027 is therefore of particular interest for 
theoretical work such as that involved in the determination of chemical composi- 
tions but in such work one cannot simply ignore the intensity anomalies for the 
weak lines of H1, Het and Heil. We cannot find any satisfactory explanation 
of these anomalies in terms of physical processes in the nebula but we have been 
able to show that the general trend is consistent with the suggestion of calibration 
errors depending on the observed quantities 7, and A. If this suggestion is correct 
the true relative emission intensities are given approximately by ¢/, where 
€=£,(1,)€2(A). 

This whole question can, of course, be settled only by further observational 
work, 

5. Helium abundances.—The condition (15) is satisfied for the three nebulae 
considered and we may therefore neglect He® in determining He/H abundance 
ratios. Results are given in Table V. 

TABLE V 
Helium abundances 
N(He*?) N(Het) 
N(H*) N(H*) 


o’o71 0°126 
0°079 0'083 
0°067 O°127 
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RADIO-SOURCE PROBLEMS 
F. Hoyle 


(Received 1959 July 29) 


Summary 


The present paper is based on the theory first proposed by Ginzburg and 
Shklovsky that all non-thermal, non-stellar, cosmic radio emission arises from 
synchrotron radiation by electrons. The electrons will be regarded as 
secondaries from nuclear collisions of cosmic rays with ambient gas. Evidence 
for this point of view can be adduced from the observed radio frequency 
spectra. ‘The consequences of cosmic rays existing (i) in galaxies, (ii) in 
clusters of galaxies, (iii) in space generally, are examined, on the assumption 
that the cosmic-ray energy density is always ~1 eV cm™. 

Apart from the introduction, the paper is divided into three parts : 

I. Weak emission from galaxies, clusters, and intergalactic space. 
Il. Strong emission from extragalactic sources. 

III. Strong emission from galactic sources. 

Halo emission from spiral galaxies—from M31 and the Galaxy in particular 
—-is considered in I. ‘The cases of Cygnus A and NGC 1275 are investigated 
in II—a new class of optically non-visible source being also proposed in this 
part. The cases of the Crab Nebula, Casseiopeia A, and IC 443 are discussed 
in IIT. 

So far as strong sources are concerned, the argument depends in an 
important way on the conditions that develop at the shock front of a cloud 
expanding at great speed into an ambient medium. ‘The present situation 
differs in a curious way from ordinary gas dynamics, in that the momentum 
and energy terms of the Rankine~Hugoniot relations are carried by quite 
different components of the assembly—the momentum by the material of 
the gas cloud and the ambient medium, and the energy by a reservoir of 
relativistic particles and by a magnetic field. ‘To overcome the difficulty that 
the detailed physical processes that take place within the shock front itself 
are not well understood, a rule is proposed that receives strong empirical 
support from a wide range of observed phenomena, including such diverse 
cases as the radio emission of Cygnus A and the optical emission of the Crab 
Nebula. 

When the magnetic field dominates in the energy terms, the rule yields 
results that differ markedly from those given by the de Hoffman-Teller 
equations. The field does not preserve a uniform structure at the front, as in 
the de Hoffman—Teller case, but develops a highly wrapped, contorted struc- 
ture—such as is suggested by the observed optical filaments of Casseiopeia A. 

If the rule were also applicable to cases where the energy and momentum 
were carried by different components, the momentum by a comparatively 
high-density, low-temperature gas and the energy by a low-density, high- 
temperature non-relativistic gas (the gas considered in the present paper is 
relativistic) a new process would be available for producing very high 
temperatures. 





Introduction.—The total energy e of an electron or positron may be written 
e=yme*, y=(1—v7/c*)—2, (1) 


where m is the mass and @ the velocity of the electron, and c is the velocity of light. 
If H,, is the magnitude of the component of the magnetic field perpendicular to v, 
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then for a highly relativistic positron or electron (y 5, 1) the rate of energy loss P(y) 
through synchrotron radiation is given by (1) 

P(y)= 1°58 x 10° H_*y% ergsec™, (2) 
where H , is in gauss. The energy is always radiated in a spectrum of the same 
form, provided a scale factor proportional to H_, y? is used in the frequency unit. 

In the following work, electron distributions of the form dy/y’ often arise. 
For such distributions the radio emission can be estimated to a satisfactory 
accuracy by calculating as if each electron emitted the whole of its radiation at a 
single frequency v given by 

v=4:19x 10°H , y?¢/s. (3) 
This simple procedure is particularly accurate when the exponent r is close to 3, 
as it will be in many of the cases examined below. 

It may be noted that for an electron distribution confined within a range 
¥1<¥Y <Y2 this method of calculation becomes inaccurate near its end points; i.e. 
for v< ~2v, Or v> ~0°5 vg, where v,, vy are given by substituting y,, v2 in (3). 
Likewise the method is inaccurate near the join of two distributions dy/y’, y <y3;3 
dy|y*, y>v4; i.e. for v< ~ 2, or v> ~0°5 vs, where v, is given by substituting ys 
in (3). 

The present point of view is that all non-thermal, non-stellar, radio emission 
arises from synchrotron radiation by electrons produced in the nuclear collisions 
of cosmic rays with ambient gas (2). For simplicity, both the cosmic rays and the 
ambient nuclei will be taken to be protons. ‘The number of cosmic rays per unit 
volume O(E) with kinetic energies between FE and E+ dE is given by 

S (4) 
where E is in BeV. Experimental evidence (3) suggests that (4) remains valid 
down to about 0-5 BeV. The mean kinetic energy over the whole distribution is 
easily shown to be ~3°5 BeV. At very high energies, E > 10° BeV, a distribution 
somewhat steeper than (4) must be used. There will be no occasion in the 
following work, however, to consider electrons derived from such very energetic 
protons. 

The multiplicity s of electrons and positrons appearing in a nuclear collision 
is somewhat uncertain, but theory and observation seem to be tolerably well 
represented by 


O(E)dE= 5x 10- 


s=2E™4, (5) 
The total energy of the electron-positron secondaries may be taken as a constant 
fraction ~ 0-1 of the proton kinetic energy (cf. 4), so that 
y = 107 E34, (6) 
with Fin BeV 
The probability of a nuclear collision of a cosmic ray per unit ttme ~cNo(E), 
where N is the number density of hydrogen atoms in the ambient medium and 
o(E)is the collision cross-section. In the following, a further simplification will be 
effected by taking o to be 2 x 10-**cm~? independent of E (5). (A slight error at 
this stage would influence later work by distorting the calculated radio spectra 
in some small degree.) 
The probability of a collision is thus 6 x 10-'® Nsec-!. With N~10~5 atom 
cm~* in intergalactic space, the probability of collision is very small over the 
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cosmological time scale ~ 10!’ sec. This is even true for heavy nuclei among the 
cosmic rays. Within the Galaxy and its halo, where N probably averages about 
10~* atom cm~*, the chance that a proton will survive for a time ~ 10" sec is still 
appreciable. 

The rate at which electrons and positrons are produced by cosmic rays between 
E and E+dE is 


6x10-*sNO(E)dE cm~-*sec—, 
a factor s being introduced to take account of the multiplicity. Substituting for 


O(E) and s from (4), (5), and expressing E, dE in terms of y, dy with the aid of 
(6), the electron—positron production rate per cm® becomes 


83 d 
~1§ x10 N r a 
Ser TG AP wi 
Since we shall not normally be concerned with values of y as low as 10%, it will 
usually be sufficient to approximate (7) by 
~1°5 x 1077! N dy/y*, (8) 
For any electron energy it is possible to define a radiation time scale t, by 
t, P(y) = yme?. (9) 
Evidently ¢, is the order of the time required for synchrotron radiation to produce 
an appreciable degradation in the electronenergy. Using (2) for P(y), and writing 


the total magnetic intensity H in place of H , (this will not affect order of magnitude 
estimates), the condition (g) for t, is readily found to be 





t= § x 10°y-? H-* sec. (10) 
Eliminating y between (3), (10) gives 
t, = 101? y-!? H-3?2 sec, (11) 


H_, being also written as H in (3). 


PART I 


WEAK EMISSION FROM GALAXIES, CLUSTERS, AND INTERGALACTIC SPACE 


1. The radio emission of the Galaxy.—The following considerations support 
the somewhat startling conclusion that the majority of cosmic ray protons have 
traversed 30 grams or more of hydrogen. ‘The implication of this conclusion in 
its relation to the cosmic ray abundances of lithium, beryllium, and boron is 
discussed in a separate paper. 


Throughout the following work it will be assumed that secondary electrons and 
positrons do not escape from the Galaxy, and that secondaries produced in the disk of 
the Galaxy can pass freely into the halo. 


It is reasonable to take the halo as a spheroidal region with polar diameter 
~30kpc and equatorial diameter ~ 50kpc, the volume being ~10%cm*. The 
present day total mass of neutral atomic hydrogen lying in the disk of the Galaxy 
~1'5 x 10°) 23x10"%gm. Such a mass distributed uniformly throughout the 
halo would yield a mean density of ~ 3 x 10-7? gmcm~*. Assuming the cosmic ray 
protons to move freely to and fro between the disk of the Galaxy and the halo, this 
gives a lower limit to the mean density of the matter through which the protons 





No. 4, 1960 Radio-source problems 341 


move; for it ignores both the contribution of molecular hydrogen in the disk and 
the contribution of ionized hydrogen in the halo. These might well double the 
mean density. Moreover in the past the total mass of gas in the Galaxy must almost 
certainly have been greater than the present day value, so that any calculation 
made over a time interval ~ 10!" sec must take account of the higher density that 
presumably was operative in the past. All in all, a mean density ~ 10-** gmcm~* 
is as good an estimate as it seems possible to make, averaging being both spatial 
and temporal. That is to say, the mean density N of the ambient medium should 
be set at ~ 10-* atom cm~* in applying (8) to the case of the Galaxy. 

A relativistic proton travels ~ 3 x 1077 cmintime ~ 10!’ sec. Atamean density 
~ 10-*6 gm cm_~ it therefore moves through ~ 30gm of material. Remembering 
that the cross-section for a nuclear collision ~ 2 x 10-**cm~, it is clear that an 
appreciable fraction of the cosmic ray protons within the Galaxy have experienced 
such a collision. 

To obtain the secondary electron distribution, consider three time scales 
t., t,, ta, Where t, ~ 10"” sec is a cosmological time scale of the order of (but less than) 
the age of the Galaxy, #, is the synchrotron life-time given by (10), and tf, is the 
lifetime of an electron against serious energy loss by atomic collisions. The 
latter time scale is given by (6) 


t,~2 x 10!4yN-"sec. (12) 


Then the electron distribution is determined to a satisfactory approximation by 
multiplying (8) by min (¢,, t,, t,), 


~1'5 x 107%! n& min (10!7, 2 x 10!#yN-1, 5 x 108H-*y-"). = (13) 
Y 


This simple procedure ignores the detailed effects of energy degradation. As 
an electron loses energy, either by atomic collisions or by synchrotron radiation, it 
augments the distribution at a lower value of y. This augmentation is not 
important, however, because of the steep dependence of the dy/y** factor in (8). 
This steep dependence means that the electron distribution cannot be appreciably 
changed by electrons falling from higher values of y. ; 

With N= 10-2 atomcm~*, H= 10~* gauss, the distribution is 


~3x 10% y-53 dy, 102? <y<5§ x 10%, (14) 
~1'5 x 10S y-*3 dy, 5 x 10?<y<5§ x 10%, (15) 
~7:5 x 1073 y-13 dy, y25 x 103, (16) 


The main contribution to the total electron energy density comes from (14) and 
(15), and amounts to ~10-ergcm~*. This is about o-1 of the energy of the 
protons, reflecting the fraction of the proton energy that passes into electron 
secondaries in a nuclear collision. 

According to (16), the total number of electrons with energies > 5 BeV is 
~6x1074cm~-%, The corresponding value for protons ~ 2 x 10~"! cm~* showing 
an electron concentration of about 1/3 per cent of the protons. At lower energies 
the proportion of electrons increases, rising to about 5 per cent at a kinetic energy of 
o's BeV. 

With H = 10-6 gauss, and v> 10’ c/s, (14) is not required, since the y values 
appropriate to such frequencies lie appreciably above 500. The emissivity is 
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given by multiplying by P(y) and by replacing y, dy by v, dv. Thus (15) gives 


dv 
~ 107% sa ergcem~*sec~!, v<10%c/s, (17) 


and (16) gives 


dv 
~ 10779 =a erg cm~*sec-!, v>10%c/s. (18) 
Vv 


The discontinuity of slope of the emissivity at the join of (17), (18) evidently 
arises from the method of calculation, which becomes inaccurate at the join of the 
distributions (15), (16). An accurate procedure would obviously give a smooth 
transition from the slope 5/6 to the slope 4/3. 

Measures by Baldwin (7) at »=8 x 10’c/s indicate an emissivity of 1-8 x 10° 
watts per pce*® per steradian per cycle per sec—i.e. 7-7 x 10-“ergcm~*sec™! 
(c/s). At v=8x10%c/s (17) gives ~2-6x 10-“ergcm~*sec~! (c/s). The 
agreement is even better than this comparison would suggest, for Baldwin’s 
value was derived on the assumption of a halo volume smaller than that used 
above by a factor ~5. The discrepancy largely vanishes when an adjustment is 
made for this difference. 


2. The steady radio emission from galaxies in general.—The radio emissivity 
might be considered to vary from galaxy to galaxy because H, N vary, or because 
secondary electrons might not be wholly trapped. 

To examine the effect of changing the magnetic intensity, consider a corres- 
ponding calculation, first for the case 

H= 1078? gauss, N=10-* atom cm-, 

In place of (15), (16) we have, 
~ 15 x 10°F dy/y93, 5 x 10?<y<5§ x 104, 
~75x 10% dy/y"8,  y>5 x 108. 
And in place of (17), (18) we now have 


34 dy a 


ys 


~10 , v<3x10°%C/s, 


30 pa 


dv 


~5 x10 v> 3x 10°%C/s, 


ian 
for the emissivity in ergcm~*sec™!. 

The effect of decreasing the magnetic intensity is thus to flatten the spectrum in 
the range 10° — 3 x 10°c/s, and to decrease the emissivity, particularly for v ~ 10% c/s 
where the reduction is by a factor of about 10. The reduction is less at higher 
frequencies. It is of interest that the emissivity at v=8x10’c/s is about 
3x10 ergcm~*sec™!(c/s)~!, in reasonable agreement with the value 
~ 10 ergcm~* sec™! (c/s)! given by Baldwin (7) for the Andromeda Nebula. 
The observed spectrum of M31 is in tolerable agreement with (21), but appears 
somewhat less steep—the observations are rather fragmentary, however. 

A stronger magnetic intensity increases the emissivity, and for frequencies less 
than ~10%c/s steepens the spectrum. Thus it can readily be shown that for 
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H=107"!? gauss, N= 10~? atom cm-*, the emissivity is given by 
p -2 dv —- a 
~ 15x 10°" —; erg cm sec (23) 


at all y>10’c/s. The emissivity at v= 8 x 10’ c/s in this case is ~ 5 x 10-* ergcm~* 
sec~! (c/s). 

The total emissivity taken over the frequency range 10’ <v< 10° is readily 
found to be ~1°5 x 107%! erg cm~* sec~? at H=10~"? gauss, ~8 x 10-* at 
H=10~ gauss, and ~ 10~*? at H=10~1%? gauss. For a halo volume of 10° cm® 
the total emission of a whole galaxy would thus fluctuate from ~ 10" ergsec™ to 
~ 10% erg sec, 

On the basis of these results it seems that the emissivity per unit volume should 
be about ten times greater in the disk of the Galaxy than it is in the halo.- Indeed 
(23) should be applicable to the disk. Because of the much smaller volume of the 
disk, 10°’ cm’, it is clear that the major contribution to the total emission must 
come from the halo, however—the disk should contribute only some ro per cent of 
the total. 

The next step is to examine the effect of reducing the amount of gas contained 
within a galaxy. In elliptical galaxies, for example, the value of N might well 
be as small as 10-4atomcm~*. Consider the case 

H=10~-* gauss, N=10-4 atom cm-*, 

The electron distribution again follows from (13), 
~ 15 x 10 8 dy/y43, 10?<y<5 x 10%, (24) 
~75x10%dy/y3, = y>5 x 10%. (25) 
The emissivity is 
dv 


~ 107% 


+t v<10%c/s, (26) 
v? 


dv 
~ 10731 = v>10%c/s, (27) 
i 


the distribution being essentially the same as in (17), (18), except that the 
coefficients are reduced by 10? because of the reductionin NV. ‘Thus the emissivity 
of a galaxy is reduced to the extent that the amount of gas within it is reduced. 


‘Two important conclusions follow : 


The steady emission from normal elliptical galaxies is very low because such 
galaxies contain comparatively little gas, and hence the ambient medium is 
ineffective in exposing the cosmic rays to nuclear collisions. 

The argument often used to overcome the lithium, beryllium, and boron 
paradox—that cosmic ray protons in the Galaxy have traversed only three or four 
grams of matter—cannot be supported (8). If this were so N would need to be 
~ 107% atom cm~’, and at this value there would be a tenfold discrepancy between 
the observed radio emission of the Galaxy and its calculated value. 

A further reason why galaxies may differ one from another was mentioned 
above—secondary electrons may not be trapped or may be only partially trapped. 
In such cases the radio emission will be reduced to the extent that secondaries 
escape from the system. It is possible that the comparatively low radio emission 
of the Magellanic Clouds is to be explained in this way. 
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3. The correlation between photographic and radio magnitudes.—A correlation 
m,—m, =constant, between the apparent radio magnitude m, and the apparent 
photographic magnitude m,, has been obtained for Sb galaxies and for Sc galaxies 
by Brown and Hazard (9, 10) and by Mills (11). It is possible that this result can 
be explained in a simple way. 

The major class of galaxy (masses ~ 10!'©)) all seem to be of about the same 
optical size (this statement is valid for ellipticals, although this fact will not be 
used here). Hence the halo systems may be expected to be of much the same size. 
The other factors promoting radio-emission—the magnetic intensity, the 
mean density N—may well also be fairly constant from one galaxy to another. 
That is to say, the absolute radio magnitudes of the galaxies probably have fairly 
small dispersion among themselves. Since the absolute photographic magnitudes 
also have small dispersion it follows immediately that variability in apparent 
magnitudes—radio or photographic—arises essentially only from variability in 
distance. This variability obviously is removed when the difference m,—m.,, is 
taken. 


4. The general emissivity within clusters of galaxies.—The emissivity within 
clusters of galaxies will be calculated on the assumptions 
(i) that the cosmic-ray spectrum is given by (4) throughout a cluster, 
(ii) that the ambient gas density ~ 3 x 10°" gmcm~3, i.e. N= 1-8 x 10-4 atom 
cm~3, 
(iii) that the magnetic intensity ~ 3 x 10~‘ gauss, 
(iv) that secondary electrons remain trapped over a time scale ~ 10" sec. 
All these assumptions tend to promote a high emissivity, so that the value calcu- 
lated below is in the nature of an upper limit. 
With N~1°8x 10-%atomcem~*, H~=3.10~7 gauss, t= 10" sec, the electron 
distribution given by (13) is 
~ 2:7 x 1077 dy/y*, 10?<y<5'5 x 108 
~1'5x10°%dy/y'3, -y25°5 x 104. (28) 
The emissivity is given once again by multiplying (28) by P(y) and by replacing 
y, dy by v, dv with the aid of (3). For H~=3 x 10-7’ gauss, the result for 
5 x 10°>v> 10’ c/s is 


ar dv P 
~2x10-% — ergem™*sec™', (29) 
v° 


which is less, at v= 10° c/s, than the emissivity in halo of the Galaxy or in M31 bya 
factor of ~ 107, 


5. The radio emission by the universe as a whole.—The present calculation will be 
based on the extreme assumption that cosmic rays with energy density ~ 1 eV cm~* 
exist everywhere throughout space. It will be shown that the magnetic intensity 
in space must then be less than 10~7 gauss in order that there shall be no conflict 
with observations. 

For convenience, H < 10~7 gauss will be assumed at the outset. Because the 
observed radio band is confined to v > 10’ c/s it follows immediately from (3) that 
we are concerned only with values of y greater than ~ 5 x 10°. Reference to (11) 
also shows that t, is greater than the cosmological time scale, ~ 10" sec, for all v up 
to ~10!'c/s. Hence for electrons yielding synchrotron radiation within the 
normal radio band there is no appreciable degradation of energy, even over the full 
cosmological time scale, t, say. 
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The spatial density of electrons and positrons therefore builds up over time f,, 
and is given simply by multiplying (8) by ¢,, 
~15 x 107% Nt, dy/y*? cm-. (30) 
Build up is of course limited to the time scale t,, since the expansion of the universe 
prevents any continued accumulation of the electron-positron reservoir on a time 
scale longer than this. ‘The general emissivity of space is now obtained by 
multiplying (30) by P(y), and by replacing y and dy by v, dv from (3), viz. 
—37 r dv wl m 
~10-87¢t.N HM — ergem "eee", (31) 
Now the general emission in space is only lost through the expansion of the 
universe. The radio energy density therefore also builds up for atimet,. * ‘Thus to 
obtain the energy density, (31) must be multiplied again by ¢,. With t,~10" sec, 
N~ 10-5 atom cm~, the energy density is therefore 


: dv 
~10~* H16 —__ ergcm~. 
pls 


By equating this expression to the Planck distribution 

8rv?c kT dv, hv <kT, 
an effective temperature 7(v) can be defined as a function of the frequency. ‘The 
observed condition that T does not exceed ~ 300 deg. K at v= 10*%c/s requires H 
to be not greater than ~ 10~’ gauss, as stated at the outset. 


PART II 


STRONG SOURCES OF EXTRAGALACTIC ORIGIN. 


1. Preliminary remarks.—It is known, or suspected, that in all intense sources 
of radio waves considerable masses of gas are in rapid dynamical motion. Ina 
series of papers Burbidge has demonstrated a remarkable result (12, 13, 14). 
Even setting the dynamical energy very high, the energy that must be present as 
cosmic ray protons is at least as great as the dynamical energy—if secondary 
electrons are to be produced in situ at a rate adequate to explain the observed radio 
emission. 

For the Cygnus A source, Burbidge finds the necessary proton energy to be 
~10* erg. For the dynamical energy to be as great as this, say for a motion of 
1000km sec~!, the mass of moving gas would have to be ~10"%©. An even 
more extreme situation seems to exist in the Crab Nebula. The mass of gas 
in this case can hardly be much greater than the Sun, and the speed of motion is 
known to be ~ 1000kmsec~!. ‘The dynamical energy cannot therefore be much 
in excess of 10ergsec~!, Yet Burbidge (4) finds the necessary proton energy to 
be appreciably greater than 10” ergs. 

Attempts have been made to argue that in some, if not all, of these intense 
sources the electrons might be primary. Because no proton concentration would 
then be necessary the energy requirement could be reduced by a factor ~ 10?, 
thereby bringing the energy of the reservoir of relativistic particles down to a more 
plausible level. But this is a council of desperation. Except perhaps in the case of 
the Crab, the frequency spectra of the intense sources are of just the form to be 
expected from electrons that are secondary to the normal cosmic ray distributions. 
To throw away this most striking agreement is to lose the one solid foundation on 
which a theory of the radio emission should be built. 
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If the cosmic ray protons are produced in situ, the difficulty is to understand 
how the dynamical energy of mass motion of the gas comes to be converted quickly 
into relativistic protons with almost 100 per cent efficiency. The problem is 
made worse by the requirements that must be placed on the magnetic field. To 
keep the total proton energy at all within bounds it is necessary that the proton 
energy and the magnetic energy shall be comparable. ‘This implies a field of 
10-*— 10-5 gauss throughout the large extended Cygnus A source and a field of 
10-?— 10~* gauss in the Crab Nebula. These values are very much above normal, 
so that the dynamical motion must not only generate cosmic ray pistons with great 
efficiency but must also be capable of building magnetic fields of large intensity. 
No such wonder-process is known. The very nature of the gas motion—radial 
expansion in the case of the Crab Nebula and probably also in the other intense 
sources—would seem to preclude the building of a magnetic field by turbulence. 
Polarization measures (15) of the Crab Nebula show indeed that the magnetic 
field does not possess the detailed, irregular structure that it would have to possess 
if the lines of force had been lengthened by turbulence. 

The attempt by Burbridge (4) to evade these difficulties in the case of the Crab 
Nebula by supposing that the proton distribution is not derived directly from the 
gas at all, but was supplied by the initial explosion that gave rise to the expanding 
gas, fails for two reasons: first, the high magnetic intensity remains unexplained ; 
second, the proton energies would long ago have been degraded. (If there was no 
gas to hold them back, the protons would stream rapidly into space without loss 
of energy. But with gas to hold them back, through the agency of a magnetic 
field, the protons simply lose energy through mechanical work done on the 
expanding gas.) 

The only possible escape from the dilemma along these lines would be to suppose 
that the relativistic protons are produced currently by a supernova remnant, as 
indeed was suggested by Oort and Walraven (15) for electrons. This suggestion 
is still unsatisfactory, however, because it leaves the magnetic field problem 
unsolved. It is also a suggestion that could scarcely be applied to the intense 
extra-galactic sources, since there are no mysterious remnants to fall back on in 
such cases. 

The present point of view is that all relativistic electrons are produced as 
secondaries from the normal cosmic ray protons, and that there is no special con- 
centration of this distribution in any of the intense sources. The energy problems 
and difficulties outlined above do not then arise at all. 

An essential step is that electrons are not produced concurrently with the 
episode of intense emission. The electron distribution is built up on a long-time 
scale, in precisely the manner discussed in the foregoing part of the paper. An 
intense source comprises a physical condition that allows the electrons of a 
pre-existing reservoir to be used quickly. The physical condition in question is 
one that operates at the surface of any rapidly expanding cloud of gas. Its 
nature will now be discussed. 


2. The shock-front condition.—A discontinuity must exist at the surface of any 
mass of gas that expands rapidly into an ambient medium. Over a limited portion 
of the surface the effects of curvature can be neglected. The ambient medium 
being of essentially uniform density, the front can be taken as moving with the 
speed V of the expanding gas. According to recent investigations (16, 17) the 
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thickness of the front is determined by the Debye length calculated for a thermal 
speed V, not by the mean free path for atomic collisions or by the Larmor radii of 
particles in a magnetic field. Oscillating electric fields within the front prevent the 
high-speed penetration of one gas into the other. Since the Debye length is 
always very small compared with the radius of the expanding surface, the front 
approximates closely to the concept of a discontinuity. 

Material from the ambient medium flows through the front, say from right 
to left. On the right of the shock the speed of flow relative to the front is V and 
is supersonic. On the left the flow speed is markedly less than V and is sub- 
sonic. A pressure gradient causes the subsonic flow to tend ultimately to zero 
on the left—i.e. the ambient material, having flowed through the front, ultimately 
joins the expanding gas cloud. i 

Write p, u, p, U for the mass density, the speed of flow relative to the front, 
the pressure, and the internal energy per unit mass, and denote values immediately 
to the right and left of the front by the subscripts 1,2. "The Rankine-Hugoniot 
relations, expressing the conservation of energy, momentum, and mass for a 
normal shock, are 

Pill = polls (mass) 
Pit pitts = pot polls” (momentum) (32) 
U, +P, /p; + 4uy? = Uz+ po/po+ dus® (energy). 


In all cases with which we shall be concerned the terms in p,, U, are small. Thus 
with u, = V we have 
pPiV = pote, 
piV? =p, + pats”, (33) 
V2=2U,+ 2po/p.+u,’, 
which solve to give 


a 2 V2 
ug = yr— >» pe= «ll ’ (34) 


i 
Pe A.’ ’ +t yt+1 


when the ratio of specific heats y, defined by 


y=1+p/pU, (35) 
takes the same value at all points behind the front, so that 


(y—1)p2U2=pr. (36) 
These equations determine the condition of material behind the shock front— 
y, V and also the density p, of the ambient medium being given. 

The above equations are very well known. What is new here is that the 
equations will be used in cases where a reservoir of relativistic particles, or a 
magnetic field, makes the main contribution to U, not the usual thermal motions 
of the ordinary ambient medium (this gas still determines p,, p, of course). 

As they pass through the front the lines of force of the magnetic field become 
plastered onto the expanding gas cloud. It will be assumed that the field suffers 
sufficient distortion for the relativistic particles to remain entirely trapped as they 
pass through the front. The following important assumption will also be made: 
whichever component—ordinary ambient gas, relativistic particles, or 
magnetic field—makes the chief contribution to U, will also make the chief 
contribution to U3. 
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The changes across the front of the minor components, whichever they may 
be, are to be calculated in accordance with an ordinary straightforward com- 
pression determined by the ratio p,/p,. 

Consider as an example the case p,=10-%gmcm~-*, H=10~* gauss, 
V=10%cmsec™!. Suppose the contribution of the ordinary thermal energy of 
the ambient gas to U, happens to be less than the contribution of the cosmic rays, 
viz ~1-6x10-"ergcm~*. Since H?/87 is also small compared with the latter 
value, it is clear that the cosmic rays then make the main contribution to Uj. 
The above rule therefore requires the relativistic particles to make the main 
contribution to U,. 

For isotropic cosmic rays, y= 4/3. Equations (34) accordingly give 


I 6 
P2=7P1» ¥e= 7 V, p= i 


Also by (36), 
8.» 
p2U = 3p2= = pl ’. 


Inserting values of p,, V we obtain p,U, ~ 2-6 x 10- erg cm~. 
The increase of energy density from 
~1:6x 10-Mergem™= to ~2°6x 10-" ergcm-* 
is due in part to the sevenfold rise of the density of the relativistic particles that 
accompanies the general rise of density from p, to p,. But a further factor ~ 23 
must come from a rise in the energy of each particle. 

The increase of the magnetic intensity is calculated from the sevenfold com- 
pression through the front. Since this compression is one-dimensional, the 
magnetic intensity increases proportionately to the density—i.e. sevenfold to 
7 x 10~* gauss. 

Gas dynamics is of no help toward understanding what physical process 
could promote such an increase in the energies of individual cosmic rays. The 
above rule implies that conditions within the front are controlled in such cases 
by the relativistic plasma rather than by the ambient gas. 

The considerations of the following sections provide strong support for the 
rule, provided we add the further requirement: 

The energy gain of relativistic particles follows a Fermi statistical dependence in 
which the proportionate gain of energy is the same for all particles. 

Two additional issues can be dealt with at the present stage, at any rate partially. 
First, the propriety arises, when the magnetic vector is not parallel to the expanding 
surface, of using the Rankine—Hugoniot relations for a normal shock. So long 
as p,, U, can be neglected from the equations, and so long as the contribution of 
the magnetic field to p,, U, is small—as it is in the above example—the obliquity 
of the shock is small. Equations (32) then represent a satisfactory approxi- 
mation. 

Even if the magnetic field makes the main contribution to the energy density, 
the above treatment is still probably substantially correct in cases where the 
magnetic energy density rises by a very large factor through the front (provided 
py U;, are still negligible). Such cases are not required in the present part of 
the paper but will arise in Part III. The very large increase of the energy density 
demands a high degree of wrapping of the lines of force, the wrapping being 
given by a process occurring within the front, again probably of a statistical type. 
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With a random wrapping of the field there is no important deviation from the 
normal condition implied by equations (32). 

A second question arises concerning the value 4/3 chosen for y in the above 
example. ‘This is correct only for relativistic particles with isotropy. There is 
no guarantee, however, that the particles gain speed equally in all directions. 
Rather would a greater gain normal to the front be expected than in transverse 
directions. But a reconsideration of the above example shows that no important 
change occurs even on the extreme hypothesis that the whole energy gain takes 
place in the normal direction. The magnetic field, plastered on the expanding 
surface, serves in this case to randomise two components of velocity, yielding 
the value y= 3/2 appropriate for a two-dimensional relativistic gas. The value 
y=3/2 must accordingly be used in (34) and (36). The rise of the energy 
density through the front is thereby reduced by 28/45 to p,.U,=8p,V?/5. Inthe 
above example the rise would therefore be by a factor ~ 100 instead of by ~ 160. 
Evidently the process is not appreciably affected, although it is slightly reduced 
in effectiveness. 

Throughout the following work the isotropic case y = 4/3 will be used. 


3. The Cygnus A case.—It will be supposed that as an outcome of a collision 
between two galaxies a mass of gas M has been set in expansive motion with 
speed V. The condition that V shall not have been much reduced by an ambient 
medium of density p, is 


4” R%p, < M, (38) 


where for convenience the cloud is assumed to be a sphere of radius R. 

We suppose further that the collision has occurred in a cluster of galaxies 
where the density p,}~ 10-**gmcm~*. This is a high value, but since we are 
now concerned with a radio source of very exceptional intensity it is permissible, 
and even desirable, to take rather abnormal values for all the quantities that affect 
the determination of the radio emissivity. An important practical advantage of 
choosing this value is that the particle density becomes close to the value used 
in Part I, Section 1, viz. N~10~*atomcm~*, Since moreover the cosmological 
time scale for building the reservoir of relativistic electrons is ~ 10! sec, again 
the same as that used in Part I, Section 1, the electron distributions (14), (15) 
and (16) can immediately be taken over to the present case. This means that 
the electron energy density in the ambient medium is ~10~ergcm-3—i.e. 
about o-1 of the proton energy density. 

Each relativistic particle experiences a rise of energy at the shock front (there 
should be no difference between electrons and protons in this respect). ‘Thus 
after passing through the front the electron energy contained in unit mass of material 
rises to 0-1 U;, where U, is given by (37). Write f for the fraction of the electron 
energy lost by synchroton radiation in the time scale R/V determined by the 
expansion of the gas cloud. Then the total radio emission rate over all frequencies 
is given simply as a fraction f of the rate at which electron energy is swept up by 


the expanding gas, viz. 
o-1fU2(47p,R°V) = hp, fR°V®, (39) 
using (37) for U,. 


These considerations can be subject to rather stringent tests, on the basis 
that V, R, f be chosen so that (39) takes its maximum value in the Cygnus A case. 
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Galaxies in clusters can collide at relative speeds up to about 3000 km sec". 
Write M,, M, for the colliding masses of gas. If no energy were lost by radiation, 
the combined gas cloud, M = M, + M,, would possess sufficient energy to expand 
outwards from its mass centre at a speed V less than the relative speed at collision 
by the factor (M,M,)"?/M. It follows that V must always be set less than half 
the relative speed of collision. Allowing for some energy loss by radiation, a 
value V ~ 1000 km sec~! would seem an entirely reasonable choice. 

The emission is at its greatest when R takes the largest value permitted by 
(38), (3M/47p,)'*. The emission that can be achieved is therefore proportional 
to M**. In considering a wholly exceptional case M should thus be set at 
~10'© corresponding to the gaseous content of a newly formed galaxy, rather 
than to an old galaxy such as M31 where M is only ~10°Q. With M=10"0, 
p,;=10-*% gm cm~%, R= (3M/47p,)'3=1-5 x 10% cm. At the distance of Cygnus 
A—7 x 10% cm, using 4 x 10! sec for the reciprocal of the Hubble constant—a 
source of this radius has an angular diameter of ~ 100”. The observed diameter is 
about 80”. 

The time scale of the expansion R/V is thus ~ 1-5 x 10! sec~ § x 10 years. 
Inserting this time scale for t, in (11), and putting H = 2 x 10~° gauss, it follows that 
f=1 for v> ~5 x 10’ c/s—i.e. over essentially the whole radio band. A magnetic 
intensity ~ 2 x 10~° gauss after compression requires ~ 3 x 10~® gauss before com- 
pression. Once again this is rather higher than might normally be expected, but 
once again this is an advantage, since all factors should take favourable values in 
the special case of Cygnus A. 

With f=1, p,=10-* gmcm~-3, R=1°5 x 10%cm, V = 10° cmsec™, (39) gives 
a total emission of ~10“ergsec-!. The observed radio flux corresponds 
to a total emission greater than this by a factor ~ 3, if Cygnus A radiates iso- 
tropically. It is certain, however, that radiation will not take place isotropically 
for an object of the present type. The initial angle made by the lines of magnetic 
force with the expanding surface must affect the detailed emissivity in a con- 
siderable measure, and it is by no means unlikely that a variation of emission with 
line of sight by as much as a factor 3 takes place. If the Earth lies in a direction 
of preferential emission (39) gives closely the correct intensity, but in any case 
(39) is correct in order of magnitude. 

It remains to consider the form of the radio frequency spectrum. The 
electron distribution in the ambient medium before compression is given by 
(14), (15) and (16). Compression across the shock front lifts the energy of each 
electron by ~20. Hence the join of (14), (15) is lifted to y~1o*. The 
emitted frequency at y=10', H=2x10~5gauss is about 10!°/c/s. The radio 
spectrum for v<10!c/s is thus given essentially by the lifted form of (14). It 
is important in this connection to notice that, since (14) contains an appreciable 
fraction of the total electron energy density, the radiated energy lies mainly 
below 10!’ c/s. 

For f=1, y< ~10*, the energy radiated by electrons between y, y+ dy is 
proportional to dy/y** (in this case the electron distribution is not multiplied by 
P(y) but by ymc* to obtain the emission). Replacing y, dy by v, dv it is thus 
easily seen that the frequency spectrum is dv/v**, in excellent agreement with 
the observed spectrum for the range 4 x 10°< v<10!c/s. The observed spectrum 
becomes less steep at v < 4 x 10°, however, and this is exactly to be expected, since 
frequencies < 4 x 10° are radiated by electrons at the very base of the distribution. 
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The problem of calculating the spectrum is then quantitatively awkward, because 
the flatter distribution (7) must be used in place of (8), and because energy 
degradation effects cannot be ignored at the base. 

The lifted form of (15) must be used for »y>10". ‘The electron distribution 
then has an extra factor y~', and this steepens the spectrum to dv/v*, 


4. The NGC 1275 case.—It is of interest to compare the Cygnus A case, in 
which p,, H, M were favourably chosen, with a case in which these quantities 
take more average values, e.g. p,=3x10°% gmcm-, H=10~ gauss, 
M=3x10°©. The value chosen for H, referring to the magnetic intensity before 
compression, is comparable with the halo fields within galaxies. It is hence a 
reasonable value for the field, say at distance 30 kpc from a galaxy within a cluster. 

The expansion velocity V will again be chosen as 1o*cmsec™', since oebser- 
vation suggests that a value of this order is attained in NGC 1275. In this 
connection it may be noted that, whereas p,, M, and possibly H also, decline as a 
cluster ages, there is no appreciable decline in the dynamical velocities with 
which galaxies may collide within the cluster. 

With N=1-8x 10°*atomem~* and H=10~* gauss, (13) gives the electron 
distribution before passage through the front, 


~ 2:7 x 107‘*dy/y*3, 1? <y<5 x 10%, 
~ 1°3 x 10-4dy/y"8, y25x 10%. (40) 


The next step is to calculate the electron distribution after passage through 
the front. From (37), p,U,=7-8x10-"ergem™ in this case. The ratio 
py U,/p,U,, with p,U, ~ 1-6 x 10-"’ ergcem~, is thus ~ 49. Again allowing for a seven- 
fold increase of density, the energy of each relativistic particle is lifted by a factor 
~7. After passage through the front the electron distribution therefore becomes 


~ 6:9 x 10~-Fdy/y, 7.107 < y < 3°5 x 104, (41) 
~ 24x 107Ndy/y"3, y>3°5 x 104. (42) 


The magnetic intensity after compression is 7 x 10~* gauss. 

The maximum value of R permitted by (38) is ~ 8 x 1072 cm, and the volume of 
the expanding cloud at the stage where R takes on this value is ~ 2 x 10% cm’. 
The time scale R/V is thus ~8 x 10'sec. With this time scale in (10), and with 
H=7 x 10 * gauss, we see that the electron energies are not seriously degraded by 
synchroton radiation up to y = 1-2 x 10‘, i.e. for frequencies up to ~5 x 10%c/s. 
‘Thus to obtain the emissivity at v <5 x 10%c/s, (41) must be multiplied by P(y) 
and y, dy must be replaced by »v, dv, 


dv , 
~ 15 x 10°91 — ergem 3sec—}, (43) 


At higher frequencies the electron distribution must be multiplied by yme? 
instead of by P(y). This steepens the spectrum to ~ 10~**dv/v*3, v> 5 x 10% ¢/s. 

The total emissivity is ~ 5 x 10°-**ergcem~*. Since the volume of the system 
~2 x 10° cm3, the total emission is therefore ~ 10" ergsec—'. It follows that the 
changes made in p,, H, M reduce the emission by ~ 10%. The reduction is 
contributed roughly equally by each of these parameters. 

The radio source NGC 1275 is an example of the present case, the total 
emission of NGC 1275 being less than Cygnus A by a factor ~ 10°. 


26 
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5. Sources not associated with galaxies.—We now take an example differing 
from Cygnus A, not in p, H, M, but in V. Let p,=10-*gmem-, 
H=3 x 10~* gauss, M = 10"'©, as before, but take V = 3 x 10’cmsec™'. Formula 
(39) with f=1 is again applicable to this case. Thus with R~ 1-5 x 107*cm as 
the maximum value permitted by (38), the maximum total emission 

~3x10"ergsec™. 

The electron distribution in the ambient medium before compression is 
again given by (14), (15) and (16). Compression across the front is less effective 
than before by a factor ~ 10, however. Thus instead of each individual electron 
energy being lifted by ~ 20, the increase is now only by a factor ~2. Hence the 
join of (14), (15) is lifted only to y ~ 10%. Hence the emitted frequency at y = 10%, 
H = 2 x 10~* gauss is about 10° c/s, the radiation spectrum now changes at v = 10* c/s 
instead of at y= 10! c/s, being of the form dv/v** for v< ~ 10° c/s and of the form 
dv/43 for v>10%c/s. 

In a former paper (18) it was shown that galaxies probably form by a cooling 
process in which dynamical motions of order 300 km sec™! are developed in the 
final stages. For a blob of cooled gas to remain compact, gravitation must be 
able to control the dynamical motions, and this requires the mass to exceed 
~10'©. If the mass < ~10''© the material simply flies apart again (except 
that an occasional blob may be able to hold together if the dynamical motions 
happen to be small in a particular locality). Here then we have a process in 
which the conditions assumed above may be satisfied fairly closely. 

In many cases the value of p, probably falls near 3 x 10-** gmcm~? (this was 
the value obtained in (18)) rather than at 10-**gmcm~*. The effect of lowering 
p, is to reduce the total emission to ~ 10 ergcm~* and to extend the spectrum 
dv/v'* down to the base of the radio spectrum. 

The present work, taken with that of the previous paper therefore leads to 
a prediction of the existence of a new class of radio source. ‘The main properties 
are: 

(i) M is high and V is low, the opposite to the case of NGC 1275. 

(ii) The frequency spectrum is of the form dv/v**. 

(iii) ‘There is no association with a visible galaxy. 

(iv) The total radio emission has an upper limit of about 10” erg sec™!. 
The majority of the unidentified radio sources could well be of this type. It is 
already of great interest that Whitfield (19) finds that unidentified sources tend 
to possess steep frequency spectra with indices of about — 1-2. 


PART III 
STRONG Sources oF GALactic ORIGIN 
1. Preliminary remarks.—Exactly the same considerations as those used in 
Part II will be employed in the investigation of non-thermal galactic sources. 
The numerical results are altered, however, in an important way, because p,, 
H, M are all changed very markedly: 
Extragalactic value Galactic value 
Pr ~3 x 10-" gmcm-4 ~3x10-%gmecm~* 
H (before compression) 10~* — 10~7 gauss ~ 10~° gauss 


M 10°— 10" © O-10°O 





No. 4, 1960 Radio-source problems 353 


The galactic values all refer to the disk of the Galaxy rather than to the halo, 
for the reason that the strong galactic sources all lie near the plane of the disk. 

‘The situation is somewhat simplified in the present case because the electron 
distribution has already been obtained in Part I, viz (14), (15) and (16). These 
formulae will be used throughout the following work. 

2. The optical emission of the Crab Nebula.—With p, = 3 x 10-**gmcm~? and 
with V given by the observed expansion velocity of the filamentary system of the 
Crab, 1-1 x 10° cmsec™', (37) gives p,U,~9°3 x 10 “ergem~*. Again taking the 
main contribution as arising from the reservoir of relativistic particles with 
p,U,= 16x 10 "ergcem~, it follows that the energy of each individual particle 
is lifted by 


: U,U, = 83 x 10%, {44) 


the factor 1/7 being due to the straightforward compression from p, to p). * 

Now the main electron energy, before passage through the shock front at 
surface of the expanding Crab, comes from vaiues of y immediately below the 
join of (14) and (15)—i.e. from values of y somewhat less than 5 x 107. After 
passage through the front, the main contribution therefore comes from y 
values of order 3 x 10°. 

It is a necessary condition (according to our rule) that the initial magnetic 
energy per unit volume shall be less than the energy of the relativistic particles. 
Taking 1°6x10-"ergcm * for the latter, this requires H< ~6x 10~* gauss 
before compression. Passage through the front increases H by the factor p/p, = 7, 
i.e. to a maximum of ~4x10~°gauss. Reference to (3) now shows that the 
emitted frequency at y=3 10°, H=4~x 10-gauss is v>1-5 x 10!%c/s. The 
synchrotron emission is thus lifted into the optical range, the bulk of the emission 
taking place at wave-lengths longer than 2000 A. 

The time scale of the expansion of the Crab is known to be ~ 3 x 10! sec, 
since this is the time that has elapsed since the supernova of A.D. 1054. Equation 
(10) gives a time scale t,~10"sec for the degradation through synchrotron 
radiation at y= 3 x 10°, H=4x 10~° gauss. Thus only a moderate fraction of the 
electron energy has been degraded by synchrotron radiation. The total energy 
of relativistic particles in the Crab (protons and electrons) is therefore given 
simply by multiplying the volume of the system by p,U/,/7. ‘The rate of emission 
is therefore 


o-1(p,U/,/7)(Volume)t,~!, (45) 


where a factor 0-1 must be included because only some ro per cent of the energy 
density is contributed by electrons. 

It is possible to proceed in two ways. First, entirely conventional values will 
be used. On the basis that the Crab has always expanded at the present rate, 
the radius R must be ~3x10!%cm, the volume must be ~10%cm', and 
the distance 1 kpe. Accepting this distance, the observed apparent magnitude 
demands a white light emission rate of ~2x10%ergcm~*. Inserting 
p,U,=9°3 x 10 ergem~*, t,= 10" sec, in (45) yields a calculated emission rate 
of 1°3 x 10% ergsec™!. 

The second method of calculating leads to very interesting results, and is 
perhaps more likely to be correct. In the Cygnus A case, where t,< R/V, the 
emission declines as soon as condition (38) is no longer satisfied—i.e. as soon 

26* 
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as the resistance of the ambient medium produces a deceleration of the expanding 
cloud. ‘There is no such decline when t,> R/V. Indeed the epoch of maximum 
emission occurs during the deceleration of the cloud, since it is not until decelera- 
tion is in progress that the energy storage in the system takes on its maximum 
value. ‘The Crab Nebula being a very outstanding object, it is possible, and even 
likely, that this phase of deceleration has now been reached. Writing V; for 
the initial expansion rate, we then require 


4 MV? (Z R) (o, U,) , (46) 


With R~3 x10" V;, and U, given by 18V,?2/49 in accordance with (37), this 
condition becomes 
M~8 x 10%p, V3. (47) 

The attractive feature of these considerations is that V; can now be set equal 
to the normal expansion rate of Type I supernovae, V;~2000kmsec™!. It has 
always been a remarkable feature of the Crab that its present expansion rate 
falls markedly below the values observed for supernovae. It is worth 
remarking in this connection that former discussions (4, 15, 20) do not provide 
this opportunity of resolving this discrepancy, since according to these discussions 
the gas cloud should be accelerating. 

With p, ~ 10-"* gmem~*, V;~ 2 x 10° cmsec™, (47) gives M~@©/3, an entirely 
plausible value. It will be noticed that p, has been set lower than the value 
3x 10°%*gmecm~* used previously. ‘This ensures that U, remains substantially 
unaltered, thereby keeping the same y values for the electron energies. ‘This is 


necessary, since the radiated energy must continue to lie mainly in the optical 
band. 

The radius Ris now ~6 x 10!%cm, the volume ~ 10** cm?, and the distance 
~2kpe. Accepting this increase of distance, the observed apparent magnitude 
demands a white light emission rate of ~ 10" ergcm-*. Inserting in (45) a 
volume 10% cm’, and maintaining p,U, ~ 9°3 x 10-* ergcm~*, for the reason stated 
in the previous paragraph, the calculated emission rate is now 


~ 1°3 x 10° ergsec™}, 

The cumulative effect of these results would seem to leave little doubt that 
the phenomenon of the Crab Nebulae is to be sought in terms of a compressive 
effect exerted on the electron reservoir of the Galaxy. ‘The energy dilemma is 
removed, the observed intensity is explained, and no difficulty arises concerning 
the magnetic field with an intensity no greater than 4 =x 10-°gauss. A simple 
explanation can also be offered for the difference between the Crab and the Tycho 
and Kepler supernovae. A similar process would be expected to operate in 
these cases also, but the degree to which the electron energies are lifted will 
be different, because p, is different in the regions where the supernovae occurred. 
A larger value of p, would lift the emission largely into the ultraviolet, while a 
smaller value would lower it into the infra-red. Moreover the Kepler and ‘Tycho 
supernovae have only about 40 per cent of the age of the Crab, so that the volumes 
of the expanding clouds are less by a factor ~ 15. On this score alone the emission 
should be reduced by a comparable factor. 

Polarization measures of the Crab show that the configuration of the magnetic 
field may be consistent with a twisted tube of force that has become plastered onto 
an expanding cloud, the field initially having much the form of a force-free filament 





No. 4, 1960 Radto-source problems 355 


(21). ‘That the field is derived from an interstellar source was already proposed 
by Shklovsky and Ginzburg. ‘This suggestion was rejected by Oort and Walraven 
(15) on the ground that compression would destroy the directivity of the field to 
an extent that would contradict the observed high degree of local polarization. 
This does not seem to be so for the systematic compression envisaged here— 
although it is a valid objection to the idea that the field is enhanced by small 
scale turbulence. It is true that some local field distortion is necessary to contain 
the relativistic particles, but a general modulation of the field would be sufficient 
for this purpose. Photographs by Baade (15, page 289 for example) show the 
white light zone to be dappled in a way that might be expected to arise from an 
undulation of the field along the envelope. 

The filamentary structure of the Crab presumably arises from the compression 
of the normal interstellar gas at the expanding surface, and hence must be taken 
as marking the outer boundary of the system. It may be noted that the detailed 
structure of the white light zone must vary according to the direction from which 
it is seen, since there is never any emission along the line of sight whenev cr r the 
magnetic vector points towards the observer. 

3. The magnetic case.—The example of the Crab Nebula shows that radio 
sources in the Galaxy do not arise in exactly the same fashion as in extragalactic 
space. ‘The enhancement of the y-values of the electrons is so great, and the 
initial magnetic intensity is so high, that (3) yields frequencies far above the radio 
wave-band. 

‘There remains the second possibility, not so far considered, that the main 
contribution to U, is made by the magnetic field. Since the cosmic rays always 
contribute ~ 1-6 107" ergem~%, the condition for the magnetic case to be 
applicable (in accordance with our rule) is 

H? , ; ‘“ 
—< s 16x 10°" ergem-*, H> ~6x 10~* gauss. (48) 

‘The mean magnetic intensity in the Galaxy is of about this value. It is hence 
clear that both the case considered previously, and the magnetic case, must take 
place within the Galaxy. Both cases may even occur on different portions of one 
and wd same ne as me, sje hoop on ad Crab Nebula is Sats a radio 


at different beasdinien on the same precisa env a. In intergalactic space, on 
the other hand, the condition (48) is presumably never satisfied, so that the 
magnetic case does not arise there. 

Condition (37) continues to hold, but with the magnetic field now making the 
main contribution to U,. Thus the magnetic intensity after passage through the 
shock front of an expanding cloud is given by 


H = 12 (7p,/7)'? V. (49) 


Putting p, = 3 x 10-* gmcm~, V = then gop 1 gives H~1-4x 10% gauss. ‘The 
magnetic intensity accordingly exceeds 10~* gauss over the parts of the surface 
of the Crab where the magnetic case is applicable (a closely similar value for H is 
obviously obtained if the alternative values p,~10 249m cm-3, V=2x 10" 
cm sec~! are used in (49)). 

In the case of the radio source Cass A a portion of a rapidly expanding cloud 
with V perhaps as high as 3000 km sec~! seems to be in collision with a more or less 
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stationary cloud (22). In this case p, is probably as large as 3 x 10 ™gmem™. 
Substitution in (49) yields a remarkably high magnetic intensity of ~ 10 * gauss 

Plainly these high magnetic intensities cannot be achieved through a one- 
dimensional compression from ambient density p, to py. Such a compression 
would raise an initial field ~ 10~° gauss to less than 10 ‘gauss. Hence in order 
that the equations (32) be satisfied in accordance with our rule, the field must 
become highly wrapped. Because of this wrapping the value of y to be used in. (34) 
can still reasonably be taken as 4/3 (if the compression had been one-dimensional 
the appropriate value would have been 2). Also because of the large tield dis 
tortion the energies of relativistic particles probably increase only by the factor 
(p,/p,)'* —rather than by p, p,, as they would for a one-dimensional gas (25) 
‘Thus the energies of the relativistic particles probably rise by no more than a factor 
~2inthis case. ‘The clectron distribution per unit volume ts therefore lifted by 
passage through the shock front from (14), (15) and (16) to 

~35 x10 "y 53 dy, y<~ 10°, (50) 

~35 x10 %y "3 dy, ~10°<y< 104, (51) 

~ 35 > 10 ly 3 dy, Y > 10', (52) 
a factor 7 being alsu introduced to allow for the increase of the electron number 
density. ‘The total electron energy density is now ~1°5 ~ 10 “ergem 4, the 
main contribution coming from y values somewhat below 10*—1.c. somewhat 
below the join of (50) and (51). 

Where conditions are not uniform over an expanding surface, in particular 
where the field is highly wrapped over separate localities on the surface, some 
unwinding of the lines of force is to be expected. Such a process, occurring 
near the boundaries of the radio and optical zones, may well be responsible for 
the light ripples that have been reported in the Crab (15). More certainly, the 
writhing filaments found near Cass A would seem to be due to just such an un- 
wrapping, taking place near the boundary of the collision zone between a mass of 
rapidly expanding gas and an interstellar cloud. 

4. The case of Cassetopeta A.V the radio emission per unit area were always 
the same over the surfaces of all expanding gas clouds, the radio surface brightness 
would always be the same, so that sources would have apparent radio magnitudes 
proportional to the square of their angular diameters. 

Now the squares of the angular diameters of Cygnus A and Cass A are in the 
ratio ~1:10. ‘The observed radio fluxes per unit frequency are, however, in 
the ratio ~ 7:10 at v= 10%c/s and in the ratio ~ 1:2 at v= 10""c/s. These two 
sources do not therefore have the same surface brightness, Cygnus A being 
brighter by ~7:1 at v= 10%c/s, and by ~5:1 atv <10"c/s, The reason for this 
difference is readily understood. 

According to the work of Part I], Section 3, the energy of cach relativistic 
particle is lifted by a factor ~23 at the expanding surface of Cygnus A. For 
galactic sources, however, the energies are lifted only by a factor ~2. ‘The ratio 


of the rates at which electron energy is swept per unit surface area is accordingly 
greater in the Cygnus A case by the factor 


~10(Vp,l eye! YP U dees (53) 


where p,l/,, 18 the energy density of the electron reservoir outside the expanding 
gas cloud. According to Part I], Section 3, (p,U 4) 


yy #8 Closely the same as p, U4 
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for the Galaxy —i.c. for Cass A, while the expansion velocity V is probably twice as 
great for Cass A. Hence (53) is a factor ~ 5, which is just the order by which the 
surface brightness of Cygnus A exceeds that of Cass A. 

It is implicit in this agreement that the time scale for synchrotron radiation ¢, 
shall be less than the expansion time scale R/V. It was already seen in Part I, 
Section 3, that this is so for Cygnus A. ‘The condition that it shall also be so in 
Cass A, even down to frequencies as low as 10% c/s, places a requirement on the 
magnetic intensity. ‘Taking R/V ~3 =< 10sec for Cass A, and using (11) at 
v~ 10" c/s, shows that if t,< R/V then H> ~2« 10 gauss. ‘This was precisely 
the order of the ficld obtained in the previous section for the case of Cass A. 

The radio spectrum obtained in Part II, Section 3, for Cygnus A was of the 
formdy/v**, ‘This should apply also for Cass A, as indeed it does. 

Finally, a direct computation of the total radio emission rate can readily be 
made. ‘The surface area of Cass A is somewhat uncertain, because the surface ts 
probably not a simple sphere, but a value ~ 5 x 10% cm* must be correct in order of 
magnitude, At an expansion velocity ~ 3000kmsec ! the volume swept out 1s 
thus ~ 10" cm* per sec. ‘The electron energy ~ 3 x 10°“ ergem ~ (this includes 
the factor ~ 2 by which the individual clectron energies are lifted), so that the rate 
of sweeping of electron energy ~ 3 « 10“ erg sec™', and for t, < R/V this gives the 
total emission rate —in excellent agreement with observations. 

5. The radio emission.of the Crab Nebula. —Vhe Crab Nebula and Cass A 
possess approximately the same angular diameter. If conditions at their surfaces 
were the same they would therefore appear as sources of equal intensity, ‘The 
observed radio fluxes are in fact in the ratio ~ 1:1 at v= 10" c¢/s, but whereas the 


spectrum of Cass A is of the form dv/v"* for vy < 10" c's the spectrum of the Crab 
isdviv"™) ‘Thus at vy 10" c/s the ratio of the observed fluxes is ~ 1: 10. 

The clear implication is that, while in Cass A the magnetic intensity 1s large 
enough for the synchrotron time scale ¢, to be < R/V for frequencies v > ~ 10% c's, 
in the Crab ¢, < R/V only for v> ~ 10’ c/s. ‘This requires the magnetic intensity 
to be greater in Cass A by about 5, a result already obtained in Section 3 from a 


consideration of equation (49). 

For vy < 10! c's, the frequency spectrum of the Crab is therefore given by first 
multiplying (50) by (7) and then by the usual procedure of replacing y, dy by v, dv 
in accordance with (3). “The result is a spectrum dv/v'* in noteworthy agreement 
with observation. 

Accepting the conventional value of 1 kpe for the distance, the surface area of 
the Crab is ~ 10“ cm*—about five times less than the surface of Cass A. Since 
the velocity of expansion of the Crab ts also less than that of Cass A by a factor ~ 2, 
the rate of sweeping of electron energy is thus greater in Cass A by ~10. A 
further factor of two arises in the total emission rate because of the less efficient 
radiation in the Crab—on account of the lower magnetic intensity. (The energy 
contained in the spectrum (10'/v)'? dv, 10° < v < 10!’ c/s, is less than the energy 
contained in (10! /v)**® dv, 10° < v < 10! c/s, by about 2.) ‘Thus the total emission 
rate for the Crab should be less than the rate ~ 3 » 10“ erg see“! calculated above 
for Cass A by a factor ~ 20-—1.c. ~1°5 x 10“ ergsec”', ‘This result is too small by 
about 5, a difference of importance. 

A number of suggestions can be made for resolving this discrepancy. ‘The 
simplest is to say that the Crab is anon-isotropic emitter (this is almost certainly the 
case) and that the Earth lies in a strongly preferential direction. ‘This possibility 
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is made plausible by the circumstance that the line of sight is almost perpendicular 
to the long axis of the Crab, which in turn points nearly along the spiral arm. 

6. The 1C 443 Case.—The effect of reducing the velocity of expansion V 1s to 
lower the magnetic intensity. ‘hus for V ~ 1ookmsec™', p, =~ 10° gmcm~* the 
value of H given by (49) is ~ 10 ‘gauss. Large, slowly expanding clouds, such 
as 1C 443 and the Cygnus Loop, may be expected to possess fields of this order. 

‘The emissivity in such cases is given by multiplying (50), (51) by P(y) and by 
writing in terms of v, dv 

dv 


~5 x10 ¥—, v< ~4x 10%C/s, 
all 


o, dv 
~ 10-8 —_ , v2 ~4 10%C/s, (54) 


ps 

the units being erg cm~*sec~!. ‘lhe emission rate per unit volume for the fre- 
quency range 4 x 10°<v<4x 10°%isthus ~10~*ergcm *sec~!. Foran expand- 
ing spherical shell of radius 15 pe and thickness 1-5 pe the volume is ~ 10°" cm’, 
yielding a total emission (for the same frequency range) of ~ 10” erg sec~'. 
Such a cloud because of its greater volume is therefore capable of emitting at a 
rate that approaches the Crab Nebula. This result is consistent with the observed 
intensities from the non-thermal galactic sources of large angular diameter—on 
the basis that brighter sources are distant about 1-2 kpc. 

7. Remarks on the shock front condition.— he above results have been obtained 
from an interesting, and perhaps unconventional, use of the Rankine-Hugoniot 
relations. Reference to (37) shows that the energy density behind the shock 
depends on the speed V of the expanding cloud and on the density p, of the 
ambient medium into which the cloud is moving. By making p,, V large this 
energy density can therefore be made essentially as large as we please. Under 
suitable conditions the energy storage behind the front takes place, not in ordinary 
material (density p,), but in an assembly of relativistic particles. ‘he energies of 
the latter particles can then be lifted essentially in any degree simply by increasing 
Py» V. 

In the absence of a thorough analysis of the physical processes taking place 
within the shock front, the present treatment must remain unproved. Yet the 
considerable successes achieved in explaining a wide range of phenomena give a 
strong reason for accepting the implications of this process. One cannot help 
wondering whether a similar effect might not arise even where the high 
speed particles are non-relativistic—in the case for example where the 
material ahead of the shock front consists of a high-density gas at temperature 
~10'deg. K mixed with a low-density gas at temperature ~ 10®deg.K. With 
the low temperature component making the main momentum contribution, and 
the high temperature component the main energy contribution, a situation might 
well exist whereby the temperature of the hot component was considerably 
enhanced, perhaps even to values ~ 10° deg. K that would be of interest in the 
thermonuclear field. A lack of collisional equilibrium between the two com- 
ponents of the assembly is of course implicit in these remarks. 

Lastly, we may mention that the present treatment of the case where the 
magnetic field supplies the dominant energy term differs considerably from that 
given by de Hoffman and ‘Teller (24, see also Helfer, 25). ‘The difference lies 
in the admission of detailed irregularities. Thus the ‘leller-de Hoffman 
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treatment assumes that uniform conditions are maintained over the shock front. 
‘The examples of Cass A and of the filamentary structure of the Crab Nebula would 
seem to argue against uniformity. 


St. John’s College, 
Cambridge : 
1959 July. 
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ON THE DIMENSIONS, SPATIAL DENSITY, AND POWER 
OF THE EXTRAGALACTIC RADIO SOURCES 


H. E. C. Tunmer 
(Communicated by F. Graham Smith) 
(Received 1959 August 11) 


Summary 


Although the majority of the discrete radio sources are unidentified, it is 
believed that they are very distant extragalactic objects. ‘They may be 
galaxies with few visible stars. ‘This paper examines physical conditions in 
large gas clouds, and shows that if these constitute the major class of radio 
sources a lower limit may be placed on their radio powers. If a large pro- 
portion of the mass is condensed into stars the conditions for gravitational 
stability, which form a part of the present argument, will not apply; how- 
ever it appears that even under these circumstances the radio power must 
exceed 10% w. ster~! (c/s)~! at 80 Mc/s. 





Introduction.—It has been shown (Ryle 1958) that the majority of radio stars, 
apart from those concentrated towards the galactic plane, must in all probability 
lic at great distances from the Galaxy. ‘This result has the most important 


implications in cosmological theory, and further consideration should be given 
to the limits which may be placed both on the average spatial density and the 
radio luminosity of these extragalactic objects. 

Basing his argument on the number of radio sources falling within given 
ranges of flux density, on their isotropic distribution, and on the observed limit 
of their maximum integrated emission, Ryle was able to show that most of the 
extragalactic objects had a radio luminosity P greater than 10% w. (c/s)~! ster! 
at 80 Mc’s and a space density p less than 10-**pe~*. An extension of his 
arguments was based on measured angular diameters of the brighter sources, 
and on the failure to identify most of them with objects; this argument suggested 
that P was greater than 3 x 10% w. (c/s)~! ster“! and p was less than 10-*° pe~*. 
This further argument depended, however, on the assumption that extra- 
galactic radio sources were likely to have physical dimensions or an optical 
brightness as great as those of normal galaxies, and although this is justified 
to some extent by the most recent identifications (Dewhirst 1959), it is desirable 
to examine this assumption further. 

The radio emission from Cygnus A apparently comes from two regions on 
either side of the visible galaxy associated with it (Jennison and Das Gupta 1956). 
‘These regions themselves are optically dark. It is possible that the association 
of radio sources with galaxies may be fortuitous and that the interaction of colliding 
gas clouds of smaller dimensions could account for a major class of radio sources 
with a high spatial density and a low luminosity. Such clouds would not 
necessarily contain any stars and might belong to a completely new class of 
astronomical objects. In this paper, therefore, we shall consider radio sources 
which are simply bodies of gas. The question of supplying sufficient energy 
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to account for the radio emission provides a basis for this discussion, and this 
factor, taken in conjunction with the requirement for gravitational stability, 
will be found to reinforce the previous arguments. 

The relation between the size and luminosity of an extragalactic radio source. 
For a typical extragalactic radio source we will use the relation (Ryle 1958) 
between P in watts (c/s)~!ster~! at 80 Mc/s and p in pe: 


\pP3? = K x 10! 


where K is found observationally to lie between limits of 2 and 4. Now it ts 
observed that very few of these sources have angular diameters greater than 
2’ arc, so that if we assign a radius R (in cm) to a source we can obtain the 
relation 

P"?R-1> 36x 107". (1) 
We may now represent possible types of source on a graph of log P against: log R 
as in Fig. 1; the above condition shows that combinations of P and R representing 
real sources must lie in the area below and to the left of the line A in this graph. 








Ikpe 


Log. R 


Fic. 1.—The limits on the power P (in watts (c/s)~* ster~!) and size R (in cms) which may be placed 
on radio sources, 


The energy contained in extragalactic sources.—It is generally accepted that 
the radio emission is the synchrotron radiation from relativistic electrons in a 
magnetic field. ‘The energy of the electrons in any given source may be calculated 
provided the field is known, but the total energy will include the energy of the 
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field itself and the energy of the other fast particles which must be present. It 
will appear that the total energy must be so high that the balance between magnetic 
energy and particle energy must be such as to give a minimum total energy; 
otherwise the supply of sufficient total energy will be impossible. ‘The precise 
balance depends on the spectrum of the electron energy, but for observed spectra 
it is near equipartition. 

In integrating the total energy of a source, it is necessary to consider the lower 
limit ¢, of energy of the electrons. ‘This is taken throughout at 10’eV which is 
just sufficient for the electrons to be counted as relativistic and therefore con- 
tributing to the radiation. By excluding electrons of energy less than e, we are 
excluding those which do not contribute to the synchrotron radiation and are 
estimating a minimum energy in the electrons. 

Since the radio emission varies on the average as y~! (Whitfield 1957) we 
infer that the number of electrons in unit energy interval varies with energy 
e as e *, and we obtain for the total electron energy 

Ke~*e de= K/e,, 

where Ke “de is the number of electrons of energy between e and e+de. ‘The 
energy of all fast particles will lie between this value and about one hundred times 
larger according to the relation of the electrons to the total cosmic-ray flux. 
We will denote this factor by QO. ‘That O is of the order of 20 is suggested by a 
consideration of the detailed relationship between cosmic rays and electrons 
(Tunmer 1959). For weak sources like the Galaxy a factor of 30 fits the 
observations. Values much higher than these involve very severe difficulties in 
the supply of energy, as will be seen later. 

Synchrotron emission from an electron flux Ae~* gives a power P at 80 Mc/s 
where 

a7P = 1-43 x 10° KH? watts (c/s)-!. 
Hence the total particle energy OK/e, = 6°25 x 10'OK is 
5°47 x 10° POH ~ ergs. 

‘The total magnetic energy is H*R*/6. If we set this equal to the total particle 
energy we have the following relations between H, P, R and the total energy FE: 

H = 1°34 x 10°O'MP14R-S4 gauss, (2) 
E=6 x 10% Q!2Pl2 R32 ergs, (3) 

‘These relations have been plotted as contour lines in Fig. 2, a graph of log P 
against log & for various values of E and H. 

An upper limit to H is provided by consideration of the observed smooth 
continuation of the radio spectrum down to 20 Mc/s or lower. ‘The electrons 
responsible for this must have energies of at least 10’eV, and their critical 
frequency v, must be below 20/3 Mc/s. 

Since v,= 4:25 x 10°H(e/mc*)? it may be deduced that’ H<4x10~% gauss. 
We can use the expression (2) above for H to show that possible combinations 
of P and R lie to the right of the lines B, depending on the value of O (Fig. 1). 

The source of the energy. —¥ rom the line B on Fig. 1 and the energy contour 
lines on Fig. 2, we can see that extragalactic sources, which according to Ryle 
must have P?> 10** w. (c/s) 'ster~!, require energy of at least 10° ergs. It is 
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apparent from ‘lable I that this could only be supplied by the translational 
energy of gas clouds with masses approaching those of galaxies. 


‘TaBe | 
The magnitudes of various possible sources of energy 
10** ergs: the gravitational energy of the Sun in the galaxy (10!! A7O) 
10°" ergs: the nuclear energy released by 1 per cent A/© of hydrogen 
turning into helium, e.g. in a supernova. 

10°" ergs: the kinetic energy of a gas cloud of 104/0© moving at 300 
km/sec, the escape velocity at the Sun’s orbit. (This is the 
same as its gravitational potential.) 

10° ergs: the kinetic energy of the gas in a galaxy (10°. MO) moving 
at 1000 km/sec. 

10°" ergs: the gravitational energy of the galaxy as a whole. 


2. $ , - s 5! 
w “q* wes | 10°°@ * 10 @* 








28 8 


£ inergs . Log wk 
H in gauss 


Fic. 2.—The minimum energy E and the magnetic field H which could obtain in sources of power 
P (in watts (c/s)~' ster”') and size R (in cms). Q is the factor by which the total energy in 
relativistic particles exceeds the energy in the relativistic electrons alone. 


The kinetic energy may be converted into a magnetic field which will grow 
both by its compression between the two bodies and by its stretching by the 
irregular motions of the gas. Subsequently particles may be accelerated by the 
irregularities in the field. On this picture it is possible that Parker’s Mach 1 
principle (Parker 1958) is applicable, in which case the energy in the particles 
approaches that in the magnetic field and the equilibrium is maintained. The 
condition for minimum total energy is then satisfied. 
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Gravitational stability.—The total energy E in the particles and magnetic 
field acts as a disruptive pressure which must be contained by the gravitational 
potential of the whole body. This is obviously true in the case of an isolated 
body of gas. In the case of colliding bodies, their kinetic energy is converted 
into fast particles and magnetic fields of E only so long as the collision continues. 
But the gas clouds will start to disperse as soon as the value of E approaches that 
of the gravitational potential. ‘The total energy EF therefore cannot exceed the 
gravitational potential of the whole mass; consequently the mass must exceed 
the mass M, where: 

aa 
= +e 
where y is the gravitational constant. From the equation for E (3) we can show 
that 
M, =3°'0 x 10°(OPR’)"*. (4) 


Spectral index.—An upper limit for the mass of the clouds forming a typical 
radio source is determined by consideration of the spectral index. It is so steep 
that the electrons must be losing energy more quickly by synchrotron radiation 
than by bremstrahlung. For this to be true the number density of hydrogen 
atoms may not exceed 10°H? (‘Tunmer 1959). ‘The total mass of gas must 
therefore be less than M,, where 


M,= = R°m10°H? 


or, using equation (2), 
My =2°1 x 107", E= 12:6 x 10 3Q'2P!2R, 5) 


This upper limit is not too low for the release of kinetic energy in the collision 
of two gas clouds to be the source of EF provided that the velocity of collision V 
is large enough. It is given by 


\M,V2q=E 


where gq is the efficiency of conversion of the kinetic energy to magnetic fields 
and fast particles. Consequently the velucity of collision must exceed 
100 g~'* km/sec, which is certainly possible unless q is very small. 

‘The upper limit /,, and the lower limit M, are only consistent if M, > M,. 
This is only true for values of P and R below the lines C on Fig. 1 along which 
M,,=M, for the extreme values of Q=1 and Q= 100. 

The optical visibility of the radio sources.—Since very few radio sources have 
been identified optically, their optical luminosities must be small. ‘They are at 
the same time quite massive. In this section, we shall estimate the free—free 
emission to be expected from the minimum mass, M_, of hydrogen and show 
under what conditions it would escape detection in photographic surveys. 

The energy of the collision is sufficient to ionize the hydrogen so that its 
temperature is at least 10 K. On the other hand, the power of the collision is 
not in general large enough to maintain the total optical radiation if the temperature 
much exceeds 10°°K. Since the power emitted in unit frequency interval at 
A=5500A is almost independent of temperature between these two values, we 
may assume a temperature of 1of K. The power emitted by a cloud of gas of 
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density N and radius R at 5500A is 1°68 x 10 N?R$ ergs/sec™! (c/s) ster. 
We invoke the gravitational stability condition (4) which may be rewritten as 
N?R® > 1-86 x 108 P!2Ol2R-12 (6) 
and hence shows that the power emitted is greater than 
3°13 x 107° P!2O!2R-!2 ergs sec—! (c/s)! ster~}. 
Using a distance determined by the radio luminosity P and the flux S, we 
can put a lower limit on the optical flux: it is 
, O'S 
° 25 e scac-l (e/s)-i -2 
3 13 x 10 Przpie TBS sec (c/s) _ s 
This can be converted to photovisual magnitude in the manner used by Woltjer 
(1958). The Sun’s m,,, is — 26-94, its area is 2°89 x 10° square seconds of arc, 
and its surface emits 9°36 x 10~°ergem~*sec~!(c/s)“! at 5500A. ‘The Sun’s 
flux at the Earth is therefore 2-02 x 10~° ergsec!cm~? (c/s)~!. 
The radio source must therefore be brighter than. m,,, where 
3°13 x roAOUZS ; oe 
= —102'4—$logS +P log— - (7) 





m,,= ~ 26-94 $ log 


2°02 x 10~P! 2R! 2 
Few radio sources having a flux density S equal to or less than 10~*4 have 
been identified with optical objects. We may conclude then that for these sources, 


m,,> +180. ‘This condition corresponds to the requirement that 


PR 


Points satisfying it lie below the lines D on Fig. 1. 

Discussion.—We have seen that various arguments limit the range of values 
of radio power P and size R which can represent typical radio sources. 

(i) The observed number-intensity relationship and the observed limits of 
angular diameter place radio sources below and to the left of line A. 

(ii) The condition for minimum total energy led to expressions for that 
energy and the magnetic field as plotted on Fig. 2. ‘The observed smoothness 
of the spectrum even at low frequencies led to an upper limit for the magnetic 
field. Radio sources therefore must lie to right of the lines By or Byyo, depending 
on the distribution of energy between the relativistic electrons and all other fast 
particles. 

(iii) Gravitational stability, together with the condition for efficient rate of 
loss of energy by the electrons, places the radio sources below the lines C. 

(iv) ‘The requirement of low optical visibility places radio sources below the 
lines D. 

Radio sources must therefore lie in the area of the P-R diagram between the 
lines A and B and below D on Fig. 1. ‘This area is shaded in Fig. 2. We can 
see that the minimum power within this area at the intersection of D and A is 
2 x 10*° watts ster !(c/s)~! or approximately 1 per cent of the power of Cygnus A. 
The minimum energy at the intersection of B and D is, from Fig. 2, 4x 10°O"? ergs. 

If more accurate measurements show that the angular diameter limit is less 
than 2’ arc, the line A will be shifted to the left and its intersection with D will be 
at even higher values of P. If the fast particle energy is not mostly in the electrons, 
i.e. O> 1, the limiting line D which is applicable is lower. Again this possibility 
increases the minimum value of P. At each stage in our discussion the extreme 
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values have been taken, i.e. the minimum stored energy and the minimum mass 
which could maintain gravitational stability. It is therefore most probable that 
the lower limit of P which we have derived is a very conservative estimate. 

We have assumed throughout this paper that it is the mass of gas alone which 
maintains gravitational stability, and so it is possible to put a lower limit on it 
(M,). ‘This argument fails however if a mass of condensed material, for example 
stars greater than M_, is allowed. Then the criterion of gravitational stability 
is that 


g? 





YMaas x Meonaensed matter >E (9) 
R 


and it is not possible to put a lower limit on the mass of gas. On the other hand, 
if the mass of condensed matter is in stars with normal luminosities, the sources 
are simply galaxies. ‘Then, from expression (g), the known ratios of the mass 
of gas to stars and the mass-luminosity ratios for different classes of galaxies, the 
optical luminosity may be estimated. ‘This gives a minimum power of 
5 x 10” watts ster! (c/s). 

It is possible that sources are short-lived and their temperatures are very 
high indeed, e.g. 10° K. The free-free emission which above a million degrees 
varies at 7~'? at A=5500A, will be less and the limits Dy and Djgy must be 
raised towards lower values of P. The limits Cy and Cj) will remain, however. 

Conclusion.—The power of radio sources must be at least of the order of 
1 per cent of the power of Cygnus A. This is true whether they are purely 
gaseous or normal galaxies. ‘They could be less powerful only if they were 
massive associations of under-luminous stars. Even with this exception in 
mind, it is reasonable to assume that P> 2 x 10*° and hence that the observed 
number-intensity relationship should be sensitive to cosmological effects. 


Mullard Radio Astronomy Observatory, 
Cavendish Laboratory, 
Cambridge : 
1959 August. 
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GRAVITATION AND THE PRINCIPLE OF STATIONARY ACTION 
C. Gilbert 
(Received 1959 September 16)* 


Summary 


A quadratic action principle is used to derive equations for pure gravita- 
tional fields. Matter is described by solutions of the field equations which 
are free from singularities and have a non-vanishing action density. A 
cosmological model is determined by this type of “‘ matter-field’’, but the 
discrete nature of matter in the universe is taken into account by replacing 
parts of the matter-field by mass particles of equal action and by free-space 
solutions of the field equations. Inertial systems are defined by the matter- 
field and the inertial mass of a particle is defined from the form of its action. 
The inertial mass is not constant, but varies as (epoch)! ®. 

The properties of the model are shown to agree with results found by 
Dirac for a cosmological model. Also the equations of motion of free particles 
in a local field are shown to be approximately the same as the Newtonian 
equations of motion under an inverse square law of attraction with the 
gravitational power of matter proportional to (epoch)~!*. For small angular 
velocities there are paths which are equiangular spirals, and it is shown how 
these paths can give rise to spiral arms in a galaxy. 

A brief discussion is given of the relation of the model to the Einstein- 
de Sitter universe, and of certain resemblances to E. A. Milne’s kinematic 
model. 





1. Introduction.—The gravitational field equations adopted by Einstein for, 
free-space were 


0, (1) 
where G;; is the Ricci tensor of a Riemannian space-time. It is well-known 
that these equations may be derived from a principle of stationary action, with 
an action density GV —g, where G is the scalar curvature and g is the determinant 
|g,;| formed from the components of the fundamental tensor g,;.__ In consequence 
of the equations (1) the action density is identically zero, so that these equations 
cannot be used to describe the field in a region containing matter, unless the 
material particles generating the field are situated at singular events. When 
matter can occur at non-singular events, we shall call the field a matter-field. In 
order to describe such a field in Einstein’s theory one has to introduce additional 
terms on the right-hand side of equation (1), which include the energy-momentum 
tensor. Since matter is composed of charged particles this tensor will in general 
include the energy of the electromagnetic field. The field equations then imply 
the principles of the conservation of mass and the conservation of charge. 

If instead of deriving field equations from the linear action principle, one 
derives them from the quadratic action principle 


5 [rv —gaV*=o, (2) 


where dV‘ is the element of coordinate 4-volume, one finds that matter-fields can 
be determined from the equations without the use of an energy-momentum 


* Received in original form 1959 January 5 
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tensor. ‘This is because, instead of the condition G=o, of the linear action 
principle, the field equations then give the condition 


AG = Oo, (3) 


where J is the operator g’/(_),,;, the comma denoting covariant differentiation, 
and solutions of the equations (3) exist for which G#o. ‘The action does not 
therefore necessarily vanish, and the field of matter can be determined from 
solutions of the field equations which are free from singularities. 

The field equations used in this paper are derived from the variational principle 
(2). They determine pure gravitational fields, and the theory is not linked with 
any specific theory of electricity, although it is possible that it may be associated 
with either the unified field theory of H. Weyl (1), or that of C. Lanczos (2). The 
matter field of a cosmological model is determined from solutions of these field 
equations. ‘The matter is assumed to exist in the form of mass particles, without 
random motions, and having a uniform average density throughout space. ‘The 
particles are assumed to be of finite size and the action of a particle in any time- 
interval is assumed to be equal to the action of that part of the matter-field which 
it replaces in the same time-interval. Each particle generates a local gravitational 
field which joins onto the matter-field at a 3-space where certain boundary 
conditions are satisfied. 

We shall make use of the postulates and principles stated below. 

A. ‘The geometry of space-time is Riemannian. 

B. The field equations are derived from the principle of stationary action (2). 

C. The cosmological principle. 

D. The principle of equivalence in a restricted form. 

E. Mach’s principle. 

F. On the small-scale the action is discontinuous, and an elementary particle 
is a world-tube in free-space, with a non-singular action density in the interior 
of the tube. 


G. The action density is invariant under transformations to conformal space- 
times. 


(The transformations are equivalent to the gauge transformations in Weyl’s 
generalized geometry (1). In the following work frequent use will be made of 
the gauge in which the action in any space-time volume is proportional to the 
invariant volume. In the case of the large-scale, or the average field, this will be 
called the natural gauge, and in the case of a local centrally symmetric field it will 
be called the local gauge. When the average field is defined by the fundamental 
tensor g;; alone, proper intervals will be said to be measured in the geodetic gauge.) 

H. When the invariant interval in the natural gauge is separated into a cosmic 
time interval and an orthogonal space interval, the world-tube of an elementary 
particle has a time-section of constant radius. 

I. The proper period of monochromatic light emitted by a particle, which is 
at rest relative to the mean motion of the matter in its immediate neighbourhood, 
is constant in the natural gauge. 

J. The path of a free particle is a geodesic of the quadratic form for the geodetic 
gauge. 

In the following discussion it will be assumed that the line element for the 
model is written in the form of equation (30), where x* are spatial coordinates of 
a co-moving coordinate system, and ¢ is cosmic time. 
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A solution of the field equations giving a matter-field for a cosmological model 
is obtained in Section 3. On account of the observed discrete nature of matter 
in the universe this solution is assumed to give the average field for a model in 
which matter has been smoothed out uniformly throughout space, and on the 
small-scale the action is assumed to be discontinuous (postulate F). The average 
field is therefore assumed to describe the large-scale properties of a model in 
which matter exists in the form of world-tubes separated by free-space solutions 
of the field equations. We consider one world-tube and join the associated 
free-space field onto the average field. The free-space field is defined by a 
Schwarzschild line element with the cosmological constant A=o, but more 
information is obtained from the boundary conditions by assuming this field to 
be the limiting case of fields with A#o. 

Boundary conditions play an essential part in the theory and the gauge trans- 
formations were first introduced in order to make certain solutions satisfy these 
conditions. From Weyl’s generalization of Riemannian geometry it is known 
that these transformations leave the field equations unchanged. The meaning 
of the boundary conditions, which are derived in Section 2, may be understood 
by considering the effect of replacing a world-tube of matter in the average field, 
bounded by the world-lines x*=constant, by the world-tube of an elementary 
particle and its associated centrally symmetric field. In so far as the principle of 
stationary action is concerned the impact of the change on the exterior average 
field comes through the 3-dimensional surface integral for small arbitrary variations 
of the field variables on the surface. ‘The question arises as to whether the 
substitution of the centrally symmetric field can be made, without affecting the 
surface integral, if certain conditions are satisfied at the surface? It is found 
that this is indeed the case, when the average field has zero spatial curvature, 
provided the additional variables given by the gauge transformations are 
introduced. 

The significance of the gauge transformations is also apparent from a different 
point of view. Onaccount of the discontinuity of the action at a 3-space separating 
different regions of space-time, some, or all, of the components of the fundamental 
tensors and their first derivatives may be discontinuous at the 3-space. There is 
therefore no a priori reason why the Riemannian geometry of these different 
regions should be based on the same standards of interval measurement. (In the 
case when the fundamental tensors and their first derivatives are continuous 
it is assumed that standards may be compared by parallel transport). By intro- 
ducing the gauge transformations the geometry of the whole of space-time may 
be described in terms of one of these standards of interval measurement. From 
this point of view postulate G is a natural consequence of postulate F. 

Postulate J may be considered to be a consequence of both D and E. Since 
mass particles do not occur as singularities of the field, the usual method of 
derivation of the path of a free particle cannot be used (3). It is possible that the 
equations of the free paths could be derived from the field equations by a suitable 
method, but without such a method at hand we assume the form of the equations 
and justify the assumption by use of the principles D and E. At any space-time 
event in the average field a transformation of the coordinates can be made such 
that the first derivatives of the metric tensor are zero in the neighbourhood of the 
event. Then by an appeal to the principle of equivalence, the general equations 
of the paths in the field may be derived as a natural generalization of these equations 


28 
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in the special theory of relativity. The principle of equivalence is not used in 
this manner to determine the equations of the free paths in the local field (this is 
the restriction on the use of the principle), because then the influence of all the 
other masses in the universe on the motion of a free particle would be neglected. 
Instead Mach’s principle is taken into account by making the paths continuous 
at the boundary with the paths in the average field. ‘This is achieved by assuming 
that the paths are geodesics of the quadratic form determined by the geodetic 
gauge of the average field. 

In Section 5 inertial systems are defined and the concept of inertial mass is 
derived from the form of the action of a mass particle. The inertial mass of a 
particle is found to be proportional to 7'*, where 7’ is the measure of the epoch 
in the inertial system, and a free particle is shown to have constant momentum. 
Also in Section 7 it is shown that the equations of motion in the neighbourhood 
of an attracting mass have approximately the Newtonian form, with the gravita- 
tional power of matter proportional to T~-'?. The paths given by this law of 
force were calculated long ago by J. H. Jeans (4), who found that certain paths 
were equiangular spirals. ‘These paths are assumed to be described by matter 
which is continually detached from the rim of a rotating nucleus of a galaxy, so 
as to form spiral arms, in the manner suggested by E. A. Milne (5, 6). ‘The equa- 
tion of a spiral arm is obtained in very nearly the same form as that found by 
Milne. 

Postulates H and I are based on the assumption that the natural gauge gives a 
fundamental system of measurement in atomic physics. It appeared necessary 
to make this assumption in order to get agreement with Dirac’s principle (7) for 
measurements made in an inertial system. In the models of relativistic cosmology 
the unit of measurement of the distance coordinate is undefined, but in the present 
case, on account of H, we are able to define the unit in terms of the radius of an 
elementary particle. ‘lhe constant ¢, which appears throughout the work makes 
the measures of proper intervals and of intervals of coordinate time in the different 
gauges, equal to the values they would have in the natural gauge at epoch t,. The 
description of the model in terms of atomic units is, however, independent of ty, 
the constant being completely absorbed in the measures of the atomic quantities. 

The main conclusions are: 


(1) A cosmological model is determined by solutions of the field equations 
derived from the quadratic action principle, which are free from singularities. 


(2) Inertial systems are determined by all the masses in the universe, in 
accordance with Mach’s principle, and in these systems the gravitational power 
of matter is proportional to T~*?. 

(3) The value of Hubble’s constant is }7~!. 

(4) The properties of the model are in agreement with results derived from 
Dirac’s principle. 

(5) The paths of free particles in a local gravitational field are geodesics of a 
metric which is conformal to the Schwarzschild metric. These paths do not 
differ appreciably from the geodesic paths in the latter field for the cases of 
astronomical interest in the solar system but, in the galaxy, matter with small 
angular velocity will spiral away from the nucleus. 


(6) The theory gives support to Milne’s gravitational theory of the formation 
of spiral arms in a galaxy. 
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Notation.—Latin suffixes take the values 1, 2, 3, 4 and Greek suffixes the values 
1, 2, 3, with summation for a repeated suffix in each case, except that when 
explicitly stated Greek suffixes have the values 2, 3, 4. The time coordinate is 
x* and spatial coordinates are x*. When spatial polar coordinates are used x! 
is the distance coordinate. 
2. The field equations and boundary conditions.—In Weyl’s theory (1) the 
metrical structure of space-time is characterized by a quadratic form 
ds? = g;,dx'dx! (4) 
and a linear form 
dl/l= x,dx', (5) 
where di is the change of length of a displacement which is carried by parallel 
transport from x’ to an infinitely near event x‘+dx'. The g,, and «, are general 
functions of the coordinates, and the forms (4) and (5) are said to define a gauge- 
system, which is invariant under transformations of the coordinates. The 
metrical structure of a domain is also invariant under gauge transformations which 
are made by replacing g,;; and x; by the quantities y*g;;, «;+ (1/x)ex/ex', where 
x is a scalar function of the coordinates. 
When the «; have the form of a gradient 4-vector 
K;= — 0a/dx'= —9;, (6) 
where o is a scalar function of the coordinates, the Weyl space is Riemannian, 
It is then possible, by a change of gauge, to make the linear form vanish everywhere 
and to transform the quadratic form into 
ds? = g,, dx'dx’, (7) 
where 
Big = Big. (8) 
The gauge will then be said to be geodetic. 
We write G,,, G and G,,, G for the Ricci tensors and curvatures of (4) and 
(7) respectively, and we assume that G is proportional to e~”. Then 
G =(G + 60; + 60%, te = nat (9) 


where A is a constant, o! = gia, and o;'= 


2 (ot = 9). 


\ —g = 


The action . of the gravitational field, in a region # of the space-time (7), will 
be defined by 


A = sa | OV=8 gdV'4. (10) 


The integral is taken throughout &, C is a constant giving the dimensions of 

action, g=|g,,| and dV* = dx! dx* dx*dx*. We find from (8), (9) and (10) that 

f is proportional to the invariant volume when the quadratic form is (4). The 

system of measurement will then be called the natural gauge, or the local gauge. 
The principle of stationary action 


e — 
Sof = ax? | OV —EdV*=o, (11) 


for small arbitrary variations 5g,,, 5(dg,,;/dx*), 50,, 5a will be used to derive field 
equations. When the field in a region # bounded by a V, is considered without 
reference to external conditions, it will be assumed that the small variations are 


28% 





372 C. Gilbert Vol. 120 


zeroon V,. On the other hand, when the fields in different regions #,, #, have 
a common boundary V3, this assumption will not be made and the boundary 
conditions will be deduced from the integrals over V3. 

The variation of the action 5. may be written as the sum of an integral J, 
throughout & and an integral J, over V3. These integrals may be readily deduced 
from expressions given by Eddington (8). If V has the equation x'=constant, 
we find 


h=C [{{— (GI bg'G) +6(—o'0! + Aetayot) — Age, 
- 120,50] V —gdV4, (12) 


a4 a. a 
I,=C | [ ev —gs (=) + 68(o'V —g)+ 1201V =380 | dS,, (13) 
1 
where 
~Y i rs -s | 

dS,=dx*dxtdx', 9,5= si (giv —g) 
and ¥ is the generalized Lagrangian formed from the Riemann-Christoffel 
symbols {ij, k}, 


L = giv —g[{il, m}{ jm, }} — {i, B{lm, m}). (14) 
When use is made of the symmetric properties of 5g;; and 5(dg;;/0x*) with respect 
to interchange of the suffixes 7 and j, we find from (13) and (14) 


1=C {| (e%e'*— dete — bere 9( 2) 


+ {t(2gingitg + agingligin — agingigin 
gotn — 
— gimgilgin — gimgilgin) =e 
— 304g! — za'g! + 301g") dg 
+ 6g°80;+ 120160 | Vv =gd5,. (15) 


Since 4g,;,50 are arbitrary, their coefficients in the integrand of (12) must be 
zero. We thus obtain the field equations, which may be written in the form 


o'=0. (17) 
These equations can also be derived from results obtained by H. A. Buchdahl 
(9). 

We now consider the boundary conditions at V, when the fields in #, and 
&, have fundamental tensors g,;,", 2;,° and 4-vectors —o,, —o. We assume 
that the fields are centrally symmetric, that o and o® are functions of x* only, 
and that V, is the hypersurface x!=a, where a is a constant and spatial polar 
coordinates are used. Then g,,%=g,,%=0, (i#4%,a=2,3), o,%=0,%=0, 
(8 =1, 2,3), and in the variational principle the variations of these components, 
and their derivatives may be set equal to zero. Let us assume for the moment 
that suitable boundary conditions have been found which enable the fields to 
be joined at x'=a in such a manner that the principle of stationary action is 
fulfilled for arbitrary small variations of the field quantities. At x!=a relations 
then exist between the field quantities for #, and those for #, and in consequence 
some, or all, of the arbitrary small variations of the field quantities for #, can 
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be expressed linearly in terms of arbitrary small variations of the field quantities 
in #,. Conversely in order to determine the boundary conditions we make an 
assumption concerning the linear dependence of the arbitrary small variations 
of the field quantities. We investigate the simple assumption that the arbitrary 
small variations of the non-vanishing quantities for the two fields satisfy the 
relations 
S09 =80, 80, =80,?, (18) 
dg = Sg, 8(0gj/dx*) = 8(0g,j/Ox*), (i=J), (19) 
bea = Seu, 5(2gq/x*) =3(8gq1/2x*). (20) 
From the principle of stationary action (11) we find that equations of the forms 
(16) and (17) must be satisfied by each of the fields, and that in addition 
I,=1,, (21) 
where J,, J, are integrals of the form (13) for each of the fields. All the small 
variations in (21) can be expressed in terms of arbitrary variations 
8¢,,, 8(ag,/axk), 80, 80,0 
by use of (18), (19) and (20), and the coefficient of each of these variations in the 
integrand must therefore vanish. This gives the boundary conditions at V3, 
which may be stated in the form that the following expressions must be continuous 
at x'=a: 


V-gg" and V—gol, (22) 
V —2(giig*! — hgikgit — } gikpit), (23) 
iV -g a| (2ehemet + agg might — agligmign! — gimgnigil — gimgnigi1)_im Bim 


ox” 


— 120/g'l — 120'g/1 + 120 | (24) 


with i=j, or 1=4, j=1 in (23) and (24). 
In order to satisfy (22) we take g,,5%=g,,% =0 at x'=a. Then we find from 


(23) that the values of V —gg'g*, («=2, 3, 4, not summed) are continuous at 


x'=a, Hence we have 

[gis] = Lele, (@=J) (25) 
where w is in general a function of x*, x3, x* and the bracket signifies that the 
values of the quantities are calculated at x'=a. Also from (24) we find when 
i=j, and i=4, j=1 respectively that 


Be a “Ese , («=2, 3,4, not summed) (26) 


and gi ey +g% 32 ee +g oe SEs _ ste + 1204, (27) 


are continuous at Pe 
If we write 
vig = Og, (28) 
where 4 is an arbitrary constant, the conditions (25), (26) and (27) are satisfied, 
provided that 


Yip» ms , Yet (i =j,a = B=2, 3,4) (29) 


are continuous at x'=a. 
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When o is discontinuous at x' =a, the metrics for #, and #, with fundamental 
tensors g,;, g,® give measures of proper intervals between any two events in 
the hypersurface x! =a, which are in the ratio exp{3(o — o®) — 3(b® — b™)} to 1. 
In order to describe the geometry of the whole of V, with one system of measure- 
ment, the gauge must be changed in at least one of the regions #, and #,. 

3. The solution for an expanding universe.—We shall now find solutions of 
the field equations (16) and (17) when (4) has the form of the line element for a 
cosmological model. It has been shown by H. P. Robertson (10) and A. G. 
Walker (11), from general considerations based on the hypotheses of homogeneity 
and isotropy, that this line element has the form 


ds? = c? dt? — Rth,, dx* dx*, (30) 


where c is a constant, R is a function of t only, and h,, is the fundamental tensor 
of a 3-space of constant curvature k= —1,0,or +1. Then the lines of parameter 
t are geodesics giving the world-lines of fundamental observers, who are at rest 
relative to the mean motion of matter in their neighbourhood, and light paths 
are the null geodesics. 
From the equations (16) we find 

5 ee nan (31) 

R R? : 


38 + 667 = Ac?, (32) 


where the dot denotes differentiation with respect to t. From (31) we find 
R? = — ke? + A®/R*+ JAR*2, (33) 
where A is a constant of integration, and from (32) and (33) we find 
6= + A/R*. (34) 
This value of 6 satisfies equation (17), and the solution of the field equations 
may therefore be found by integrating (33) and (34) to obtain R and o as 
functions of t. It will, however, be more convenient at present to use a dimen- 
sionless coordinate o instead of the time coordinate t. Then writing R(t) =R’(c), 
the metric (30) becomes 
ds* = cy? R’*do? — R'*h,,, dx*dx*, (35) 


where c,=c/A is a constant having the dimensions of a length. 
Also from (33) and (34) 


where 
a? = |A|co?/3, (37) 
and the upper and lower signs in the last term of (36) correspond to positive 
and negative values of A respectively. 
In the case k=o, we find the following solutions: 


(i) A>o, R’=asinh3(d—o), (38) 
(ii) A<o, R’=acosh3(d—c), (39) 
(ii) A=o, R’=e-%, (40) 
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where d,d’ are arbitrary constants. Case (iii) may be considered to be the 
limiting case of (i), or (ii), when ao and do so that 
Limit (e*4a/2)=e8%’. (41) 
a=0,d=a@ 
If we write 
eVv=ty, e%=(t/ty)'*, (42) 
the metric (35) becomes 
ds? = c2dt® — t?3k,, dxdx*, (43) 
where ¢=c,/3 and k,, is the fundamental tensor of the 3-space withk=o. Then 
writing 


T= jtto 8, Ty= ito, (44) 
the metric for the geodetic gauge is 


7) kg dxedx*. (45) 
Ty ; 

4. The mass particle in an expanding universe.—When the quadratic form (4) 
has a centrally symmetric form with g;; independent of x‘, it is found that equa- 
tions (16) and (17) have solutions with o,=constant. We consider only the 
case o,=0 which gives the Schwarzschild exterior solution of Einstein’s equations 
for free-space. ‘Taking spatial polar coordinates x!=p, x?=0, x*=¢, and a 
time coordinate x*= 7, the metric can be written 


ds’ = 2 dr? —e %dp? — p*d6® — p®sin? 0 d¢?, (46) 


where 


(47) 


In this solution A’ has been written for A in equation (16), and I’, m are constants. 
In Einstein’s theory [ is the gravitational constant and m the gravitational 
mass, but we shall use these terms to describe different quantities later on. 

We shall now show that, when k=o, the field of the mass particle can be 
joined onto the field of the expanding universe so as to satisfy the boundary 
conditions (29). Let the 3-space with fundamental tensor k,, be described in 
terms of polar coordinates (r,4,¢), and let b=o in (28). Then y,; can be taken 
to be the fundamental tensor of a quadratic form 


ds’? = c2dt'? — S(t’ )?(dr? + r2d6? + r? sin 76 dd*), (48) 

where 6, e*°R'(c)'do=cdt’, e*R’'(c)= S(t’). (49) 
From (49) we find for both the cases (i) and (ii) given by (38) and (39), that 

2 ' 3d 
a (SF 9) = 30 + ae-*4S(2")?, (50) 
On the other hand y;; for the local field is the fundamental tensor of (46), since 
for this field o=constant. 
The conditions (29) which must be satisfied by y,; are included in the 


conditions derived by S. O’Brien and J. L. Synge (12). It has been shown 
(13) that the latter conditions can be satisfied at a 3-space r=a, by taking 


t=t(t',r), p=S(t')r, (51) 


dS(t’)\?_ 2m | NAG ng 
(Fr) = a5) > (52) 





provided that 
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Comparing (50) and (52) we find that the conditions (29) are satisfied at r=a, 
provided that 

Do? .. .« Aq? 
aa’ ated a =" 

A’ is therefore positive, and from (37) and (53) we deduce 


ulti (53) 


wa lAl Pm _ _ gel 


KN? VGA) - 


In the limiting case A=o, A’=o we find from (41), (42) and (54) that 


Limit -S4l 
A=0, A’=0 V(3A’) 
ad 
a ty 


2 
a (55) 
(56) 


In this case we also find from (40), (42) and (49) that we can write 
t’ = t*/2t, (57) 


and hence the quadratic form (48) becomes 


2/3 
ds’? = c2dt'2 — (*) t'*3k, dx*dx*. (58) 


This metric gives measures of intervals in the local gauge, and is similar to the 
metric of an Einstein-de Sitter universe (14) but the description in the geodetic 
gauge will be shown to be more significant, and this will lead to different properties 
of the model. 

5. Inertial systems.—E. Mach held the view that the inertia of matter is due 
to the presence of other matter, and that the inertial frames are determined by 
all the masses in the universe (15). These ideas of Mach are incorporated ina 
cosmological model by the following steps: 

(i) The matter-field solution of Section 3 is assumed to give the average 
field of a cosmological model consisting of mass particles, and an inertial system 
is defined to be a system of reference at a given event for which spatial and 
temporal coordinates are natural coordinates for the average field in the geodetic 
gauge. 

(ii) Mass particles and their associated Schwarzschild fields replace the 
average field in regions of space-time, in the manner described in Section 4, 
and the action of a particle is equal to the action of that part of the average field 
which it replaces. 

(iii) The inertial mass of a particle is defined by making the action of a 
particle in a local inertial system similar to action in mechanics. 

We consider the limiting case A=o0, A’ =o for which the average field in the 
geodetic gauge is defined by (45). Then at an event having coordinates (7, x*) 
we can define local coordinates 

T\} 


p-uge(Z) 0) 
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such that the metric locally has the form 
ds? = c*dT? — (dé"" + dé* + dé"). (60) 
Then the metric (60) will define the local inertial system at the event*. 

In the general formula (10) for the action ./, we can choose any value of C 
which will give « the form of action in mechanics, but the same value of C 
must be used in different regions #,, #, in order that the boundary conditions 
shall be applicable. We find it convenient to take the value 

aA 
C= eT A" (61) 
Then in the limiting case A =0, A’ =o we find that there is a non-vanishing action 
density for the average field. Calculating the action in the region between the 
hypersurfaces x! =a, T= 7,, T=T,, we find in this limiting case from (9), (10), 


(43) and (44) 
— A2 T. 7 1/2 
of =49ca%ty Limit (rr) | (z-) dT. (62) 
A=0, A’=0 PA T; T 


Using (55), (56) and remembering that c= c)/3, we find 
T™ /T\12 
oP | - (x) cd. (63) 
tr, \To 


Since dT is the interval of proper cotime in the local inertial system, the equation 

(63) shows that of has the form of action in mechanics if the inertial mass of 
the matter in the 3-space is 

m,=m(T|T,)*2. (64) 

When the average field in the region is replaced by a mass particle and its 

associated Schwarzschild field, the action of the latter field is zero. This is 

because the action is equal to the invariant volume occupied by the field multi- 

plied by the factor limit (2CA’), which vanishes on account of (61). Hence 
A=0, A’=0 

all the action (63) is equal to the action of the particle and m; is identified with the 
inertial mass of the particle. 

By the postulate J the equations of the free paths in the field are given by 


5 [di=o, (65) 


From the metric (45) we find that these equations give 

de Tw to??T , dx* dxP\12 

7 eT | ° ae oir aT) : 
where v* are arbitrary constants. 

From (59), (64) and (66) we can write 
mV 
va-Vie) (7) 

where p is a constant 3-vector and V is the velocity 3-vector at the event x*=o0, 
having components dé*/dT and magnitude V. 


(66) 


* The coordinates 7, £* are not natural coordinates according to the usual definition (8, p. 80), 
but the difference is unimportant for the following work. A transformation of the time co- 
ordinate T=7,+ T’|c—(€"* + €* + &")/4c°T, is required in addition to (59), in order to make the 
first derivatives of Z,, zero at the event T=7,, £*=0, and also to make the velocity of light unity. 
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Equation (67) expresses the conservation of linear momentum according 
to the special theory of relativity, but on account of the variation of m, with T 
it is easily seen that the velocity of a free particle decreases as T increases. On 
the other hand, the energy E of the particle given by the special relativity formula 


mc — 
V(1—V8c8) 
is found to increase with 7' on account of this variation. 

In the neighbourhood of a fundamental particle of inertial mass m,, a free 
particle will be acted on by a gravitational force in an inertial system. The 
equations of the paths in the local field are derived in Section 7 where it is shown 
(equation (105)) that, for sufficiently small values of €, there is a Newtonian 
attractive force on the particle of intensity 


in! 6 
PAyy 


In Newtonian theory it is assumed that gravitational mass m, and inertial 
mass are equal. In order to preserve the equality m,=m, we shall therefore 
write the intensity of the force in the form ym,/&, where the gravitational power 


of matter is 
a 2a 
nt ao bet oo. 6 
Y ( ?) — (69) 


y rather than I’ is the measure at any epoch of the Newtonian constant of 
gravitation. 

Also since inertial mass m; is enclosed in a spherical surface of radius 
t,'3(T/T,)'*a at epoch T, we find that the mean density of matter in the universe 


is 
= 3% _(Te\" _ _9F_ 
as oral 7) 167a° T° (70) 


6. Dirac’s principle.—The cosmological theory of Dirac (7) was based on the 
principle that ‘‘ any two of the very large dimensionless numbers occurring in 
nature are connected by a simple mathematical relation in which the coefficients 
are of the order of magnitude unity’’. Since the values of the gravitational 
constant, the mean density of matter in the universe and Hubble’s constant, 
when made dimensionless by the use of constants of atomic theory, were found 
to give very large numbers of the same order as the reciprocal of the very large 
number ¢ formed from the age of the universe and an atomic unit of time, 
Dirac concluded that the first three quantities were approximately equal to t~!. 
In order to describe a cosmological model by the equations of the general theory 
of relativity he then changed the units of time and distance so that the gravita- 
tional ‘‘constant’’’ was an absolute constant and the velocity of light was 
unchanged. It was later pointed out by H. Bondi (16), that the latter procedure 
was not justified, because the change of units would make the value of the 
electronic charge vary with the time, and this would be contrary to the principle 
of the conservation of charge which is inherent in the general theory of relativity. 
Nevertheless it would appear necessary to take into account the results derived 
from Dirac’s principle in constructing a cosmological model. Theories in 
which the gravitational power of matter varies inversely as the epoch have also 
been proposed by P. Jordan (17) and R. H. Dicke (18). 


E, (68) 
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We shall show that the results derived from Dirac’s principle are obtained 
from measurements in an inertial system, in consequence of postulate H, when 
appropriate units of the time and distance coordinates are chosen. Let these 
units be chosen so that a is the radius in the natural gauge of an elementary 
particle of inertial mass m,, and 


c= 4/3.4. (71) 
Then, in the geodetic gauge of the inertial system, we can define atomic units 
of length and time by 


T 14 a T 14 
s=a(7-) = Sz) ; (72) 


Also from (44), (71) and (72) a dimensionless measure of time in the inertial 
system is given by 
T/e=cla. T**T,"* =t. (73) 
Then from (64), (69), (70), (71), (72) and (73) we find : 


(74) 


(75) 


On account of postulate I we calculate Hubble’s constant from the metric 
(43) for the natural gauge. We find in the usual manner (19) the value }¢-'. 
Since coordinate time in the natural gauge is by (73) equal to the dimensionless 
quantity 7/e, the value of Hubble’s constant H in the inertial system is 


H= \(et)™. (76) 


At epoch t=1 proper distance in the natural gauge is equal to coordinate 
distance, and consequently the average field of the model is a complete matter- 
field, in the sense that there are no free-space regions. The description of the 
model in terms of discrete masses will therefore only hold for t>1. We then 
find that the Schwarzschild field is free from singularities. The condition for 
this is that 2I'm/c2p, <1, where p, is the radius of an elementary particle in the 
local gauge. Since py=(7/T,)"*6, we find from (56) and (44) that 

i as (77) 
C“py t 

Hence the free space regions are free from singularities. 
If we write y*, w* H* for the values of y,w, H when they are made dimen- 


sionless by the use of the quantities m;,5,« we find from (74), (75) and (76) 


*— a *— - e ¢— s ~ 8 
eo ee " 3t (78) 
Since tf is the epoch T divided by the unit of time «, these results agree with the 
deductions from Dirac’s principle (7), if m,, 5 and « are assumed to be quantities 
occurring in atomic theory. 
From (72), (74) and (76) we find 


(79) 
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and since the quantity on the right-hand side of (79) has the dimensions of length, 
we can use the values of the quantities in c.g.s. units to calculate the value of 5. 
Assuming that m; is the neutron mass, and that H = 75 km/sec/megaparsec, we find 
that 6 is approximately 10~'? cms, which is of the same order of magnitude as the 
nuclear radius. This result justifies the assumption that m,, 5, « are the measures 
of quantities occurring in atomic theory. 

7. Spiral orbits in a local gravitational field.—The local field of a mass M 
has the form 

dt=2 (: * ez 
c*p 

This field will be joined onto the average field in the local gauge defined by (58), 
at a 3-space r=b, by the use of boundary conditions (29). These boundary 
conditions are satisfied by the method described in Section (4). 

We make the transformation 


2rM\",, 
jan (:- =) dp*—ptd0*—p*sin*@dg*. (80) 
'p 


r=7(t',r), p= (4) "rm, (81) 


to 


a(t a) Or _ op op (82) 


c2p } ot’ or at’ ar’ 
so that the metric (80) has the form 
ds’? = c?e*dt'? — ebidr? — (4/t,)**t'**r?(dé* + sin? Ad¢*), (83) 


where ¥,,\, are functions of t’,r. When (82) is solved for 7, subject to the 
boundary conditions at r=5, the functions ¥,,%, can be calculated from 
expressions involving 07/ét', dr/dr. 
We find from the boundary conditions that 

[rM_ 8 

aaa OE ates 8 

ra (84) 
which corresponds to the condition (56). Also when #, and ¢, are expressed 
in terms of the variables 7, r and the constants 7), t,,b, we find, after a lengthy 
but straightforward calculation, that 


r 2 y 2 
on1e(Z) {r+ ¥(5 ‘ aD} (85) 
0 
V ae 


(86) 


(87) 


T is defined by (44) and terms of higher order in N/T are neglected. Since 
N/T =20'M/c*p,, where p, is the radius of the boundary in the local gauge, the 
neglected terms will be <1. 

We shall now determine the free paths in the local field. The metric in the 
local gauge is given by (44), (57), (83), (85) and (86), but the free paths are 
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geodesics in the geodetic gauge, and the metric is then given by the conformal 
space-time which is continuous with (45) atr=b. We find 

dg? = e’"dT? — e*dr* — e*(d@ + sin® 6 d¢*), (88) 


e = 4Tt-¥9 fi + 7(; ie ae (89) 
e = Tt,- "7, (90) 


rf 


eae a —}>- (91) 
The Christoffel symbols for the metric (88) were calculated from general 
formulae (19), and the equations of the geodesics were then found to be 


d’r dr\? dr dT d¢\? 
+4 Ce. he 29(° 
a 4)’ (5) +X5 a bet 4’ sin (=) 


~4e'w' () +e 0 (F) =0, (92) 


do ~—s, dr d d¢ dé dT 
a ui GF -sindcos0 (FE) + ji bea 


ddr dé 


where 


=0, (93) 


e2cor9 0% , ab aT 


mt aa Za’ "aa (94) 


(3) 


-G) {ay -eme(S)- 


Considering paths which satisfy the initial conditions 6=2/2, d6/d§ =o, we find 
from (93) that d*6/ds*=0. These paths will therefore always lie in the plane 
6=7/2. In this case we find from (94) 


Tr + =ha*!c, (96) 


where A is a constant of integration. Then from (g2) and (95) we find 


d% A dr dr eae dT > oe 
on) goes CE —Saee om 


(2) (se SRE) eC) 


We obtain an approximate form of the equations by writing 


pee. epee (99) 
dT dq’ " a’ ” 
and expressing the terms in powers of T-'. Equating to zero the coefficients 


of the lowest powers of 7-1 in (96) and (97), we find 


no =f, (100) 


d™, 


oa - +45 - G0 (101) 


n= M/m. (102) 





382 C. Gilbert Vol. 120 


In deriving these results use has been made of equations (56), (84), (87), (96), 
(97) and (98). When » is small compared with the other terms in (101), these 
are the Newtonian equations of motion under a central force —n/2n*, and the 
paths are therefore conics. The solution in the more general case has been 
discussed by Jeans (4). 

If we use proper cotime T and proper distance 


*\L2 
g=1 (7) = Td, (103) 


0 


where d= 2at,~"*/,/3, the equations (100) and (101) transform into 
) V3 q 


dp _ 
fis =H, (104) 


7 M 
dT? i ne 
The constant H=hd? and # is given by 


. 1/2 
an ay u (=)  s (106) 


(105) 


where .V, is the inertial mass of the attracting body and y is given by (69). For 
small values of € the term —£/47? in (105) is negligibly small and then the 
equations of motion are the same as the Newtonian equations of motion under a 
central attraction. #/&, where. #@ccT-**. The nature of the orbits in this case 
have also been discussed by Jeans (4). In particular there are orbits given by 
equations (100) and (101) for which 7 = constant, and these give the equiangular 
spirals 
E=F, en, (107) 
where &,€, are constants and 
I 


ph dT ~ 2¢T 


The angle «, measures the inclination of the spiral at any point to the circle 
through the point, whose centre is at the attracting mass. From (108) we see 
that this angle will only be appreciable when ¢ is extremely small, of order 
T-', or less. 

When r<bé terms of the next higher order in N/T have to be retained in 
forming the equations (100) and (101). ‘These terms come from the approxi- 
mate expressions for e* and e” given by (89) and (g1), the contributions from 
the new terms of the higher order approximations for e* and e” being negligibly 
small in comparison. From the equations the calculated values of the cen- 
tennial advance of the perihelion of Mercury and the gravitational deflection of 
light by the Sun agree with the results of Einstein’s theory. Also the calculated 
value of the gravitational shift of the spectral lines, for light emitted by atoms 
at the surface of a star, is found, from (88) and (g1), in conjunction with postulate 
I, to be in agreement with the result for the latter theory. 

The theory of the structure of a spiral nebula can be developed from the 
equations (100) and (101) along the lines proposed by E. A. Milne (5, 6). We 
assume that the nucleus of the nebula has a constant radius x, in the natural 
gauge. This assumption is made because of its similarity to the postulate H 


(108) 
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for an elementary particle, although a closer resemblance to Milne’s results 
would be obtained by assuming the radius to be constant in the local gauge (cf. 
Section 8). At epoch 7, let »=7, for events on the boundary of the nucleus, 
so that using (43) and (44), we find 


Xi = ToT nd. (109) 


Then at epoch T the radius of the nucleus in the geodetic gauge of the inertial 
system is 


€,= TMT “yd. (110) 


We assume that the nucleus is rotating with constant angular velocity, and 
that matter is continually detached from a point (£,, 7/2,¢,), with fixed angular _ 
coordinates, on the rim of the nucleus. We assume that each portion of matter 
detached describes a circular orbit y=», where the value of 7, depends on the * 
epoch 7, at which it became detached. At epoch 7 this matter has coordinates 
(€,7/2,¢), where 

£€= T™*y,d. (111) 


Also from the equations of motion (100) and (101), we find 


é~&s “log (7); (ve2) 


h* n 
Mis ia (113) 


Eliminating h,»,, 7, from equations (110)-(113) and using (106), we find the 
equation of the spiral arm at epoch T in the form 


T2 12 
b-6=4( = :) log (=). (114) 


c 
$ $1 


This may be compared with the corresponding equation found by Milne, 
(5,6) which in our notation becomes 


t-te (eH) "i 9 


In the kinematic treatment, however, y is proportional to T and M, is a constant. 
At distances from the nucleus which are small compared with the distance 
&=T1*ba-'d, which is the radius of the sphere of influence of the local field, 
we find yM,77€-*> 4. Consequently Milne’s general conclusions regarding his 
equation for a spiral arm will apply also to the equation (114). 

8. Discussion.—We have seen in Section 4 that in the local gauge the average 
field has the same metric as an Einstein-de Sitter universe. Also from (44), 
(57) and (63) we find that if proper intervals are defined by the metric (58), the 
action .& will have the same form as action in mechanics if m is defined to be the 
inertial mass of a particle. The question therefore arises as to whether the system 
of measurement which corresponds to measurements in the physical universe 
should be defined by the local gauge rather than the geodetic gauge. This also 
raises the question of the validity of postulate J, because then from the principle 
of equivalence we should be led to make the assumption that free paths are 
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geodesics of the metric (58). There would then be no essential difference 
between the properties of this model and the Einstein-de Sitter universe. There 
are at least two reasons against making this interpretation. One is that it does 
not appear to admit of a satisfactory derivation of Dirac’s results, and the other 
that it does not incorporate Mach’s principle into the theory. 

We consider the second reason from the point of view of Newtonian cos- 
mology (16). It was shown by E. A. Milne (20) that there is an exact corres- 
pondence between the Einstein-de Sitter universe and a Newtonian universe 
described in Euclidean space, with distauce coordinate p and Newtonian time 
t’ defined by (51) and (57) respectively. The position vector of a point in this 


space may be written 
1/3 
Pp == = t'? 3p, (1 16) 
ty 


where r is the position vector of the corresponding point in the subspace 
du* =k,, dxedx® of the space-time with metric (58). ‘The Newtonian model 
satisfies the principle of relativity for transformations to any of the fundamental 
observers r=constant. In Milne’s description with Newtonian time ?’, these 
fundamental observers were relatively accelerated, but if Newtonian time is 
defined differently by 


T = 3¢,490’8?/242 (117) 
they have uniform relative velocities 


V= (118) 


4 
T 


Now consider the motion of a free particle relative to a local fundamental 
observer, the relative velocity being small compared with the velocity of light. 
If we use postulate J and the integrals (66) of the geodesic equations, we find 
that the local velocity of the particle is dp/d7~v where v is the constant vector 
with components v*. In the given frame of reference the particle accelerates 
from one fundamental observer to another, but it always has approximately the 
velocity v relative to the frame of reference of a fundamental observer in its 
immediate neighbourhood. We can picture the motion as consisting of a 
uniform velocity relative to a uniformly expanding space. The frames of 
reference are thus a generalization, for an expanding space, of the Newtonian 
concept of inertial frames. On the other hand, if we postulate that the paths 
are geodesics of the metric (58) we find that the local velocity of the particle is 
de/dT ~vT,'"/T**. The particle then comes to rest relative to the expanding 
space as 7’ tends to infinity. ‘The former postulate rather than the latter is 
therefore more in accordance with Mach’s principle. 

The description of the model in terms of measurements made in the geodetic 
gauge therefore seems preferable to the other alternative. The fact that the 
inertial mass of a body increases with the epoch, does not imply the creation of 
particles, since the masses of elementary particles also vary with the epoch in 
the same manner. 

The description of the model in Euclidean space with distance coordinate 
p and Newtonian time T brings out certain resemblances to Milne’s kinematic 
model. The velocity-distance relation (118) is the same as that found by Milne 
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for his system of fundamental observers. Also from the equations (100) and 
(101) we readily deduce 


» db — 
27 —f 
par IT, (119) 


(hy ee ee pe 
dT ° (7 ae (7): " 
where H’ =a*t,**h/T,?. These equations resemble the equations of motion 
obtained by Milne for a particle in the neighbourhood of an attracting 
mass (§). A slight difference occurs in the first term on the right-hand side of 
equation (120), where we have —p,27*, whereas Milne’s equation has —p/T?. . 
‘The difference can be traced to the occurrence of the term — }y in equation (101), 
and arises from an essential difference in treatment of the problem. In the ° 
present theory the gravitational field of the mass .V merges into the background 
field of the universe, whereas in Milne’s theory the gravitational field of the 
condensation is superposed on the field of his substratum. _ It is on account of the 
close resemblance of equations (119) and (120) to Milne’s equations that the 
equation of a spiral arm can be obtained in very nearly the form (115) if 
p=constant at the boundary of the nucleus. 

The nature of the field in the interior of an elementary particle, or a nebular 
nucleus, has been left undetermined. ‘This problem may require a more complex 
field theory than has been used here, but the following property of the fields 
should be mentioned, as indicating that the matter-fields of Section 3 may be 
significant for the description of fields in the interior of matter. Using the 
notation of equations (19) and (20), we find from (30), (33) and (46) by the 
method of Section 4, that 

Ag, = A'g, (121) 


A'\ES /2P'm\"3 
r -(<) (Gr) R. (122) 


at an expanding boundary 


If in the equation (10) defining ./ we had omitted the factor A~', we should have 
arrived at the boundary conditions (29) by taking Adg,;Y = A’ég; ete. in 
equations (19) and (20). It seems possible that these conditions on the small 
variations may therefore be related to the continuity relations (121) at the 
boundary of the particle, whereas the derived conditions (29) express the 
continuity relations with the average field. 


Kings College, 
Newcastle -upon- Tyne : 
1959 September 15. 
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